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ABSTRACT 

Complementary  binary  sequences  were  invented  by  Golay  in  an  in- 
vestigation of  infra-red  multi-slit  spectrometry.      This  dissertation 
formalizes  the  basic   results  obtained  by  Golay  and  develops  new  concepts 
and  techniques  for  examining  the  characteristics  of  these  special  binary 
codes.     This  work  has  developed  new  understanding  of  the  structure  and 
methods  for  the  decomposition  of  complementary  sequences. 

Complementary  sequences  have  the  property  of  an  infinite  correla- 
tion peak  to  ambiguity  ratio  when  detected  with  a  matched  filter.      These 
binary  sequences  should  find  much  application  as  pseudo- random  noise 
modulation  signals  for  both  radar  and  communications  systems. 

A  discussion  of  the  need  for  such  sequences  is  included  in  the  intro- 
duction and  is  followed  by  a  state  of  the  art  description.      An  operations 
group  on  the   sequences  is  formulated  and  the  proofs  of  several  theorems 
concerning  the  operations  group  are  given  in  a  rigorous  manner.      One 
reason  for  developing  the  operations  group  is  the  application  to  elimina- 
tion of  redundance  in  the  computer  search  for  new  codes. 

Several  invariant  properties  of  complementary  codes  are  proved 
through  the  use  of  the  Hamming  distance  concept.     Many  more  invariant 
properties  of  the   sequences  are  demonstrated  through  the  introduction  of 
a  Hamming  vector.      The  concept  of  a  Hamming  vector  is  extremely  use- 
ful as  a  complementary  code  decomposition  tool.     A  large  number  of 
theorems  are  proved  to  enhance  its  use  in  this  field. 

Several  computer  searches  for  complementary  sequences  are 
described  and  the  actual  computer  programs  for  the  CDC   1604  are  in- 
cluded in  the  Appendix. 
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CHAPTER    I 
INTRODUCTION 

Search  radars  of  today  are  required  to  obtain  smaller  and  smaller 
targets  at  ever  increasing  ranges.      This  has  led  to  many  different  tech- 
niques of  summing  radio  frequency  energy  to  obtain  these  long   ranges. 
Complementary  sequences  are  well  adapted  for  coding  the  RF  energy 
pulses  to  give  a  summing  of  the  returning  RF  energy  in  the  echo.      The 
present  work  is  more  concerned  with  the  characteristics  of  complemen- 
tary codes  such  as  their  operation  groups,    Hamming  distances  and 
Hamming  vectors,   than  it  is  with  their  application  as  modulating  signals 
for  RF  carriers.     However,   the  introduction  will  develop  the  radar 
problem  as  a  vehicle  to  emphasize  the  use  of  these  codes  in  future  appli- 
cations in  radar.      Their  application  will  probably  be  equally  important 
in  the  communications  field.      The  initial  work  in  complementary  sequen- 
ces was  directed  toward  an  infra-red  application,   but  it  is  felt  that  the 
eventual  application  of  these  sequences  will  spread  to  many  different 
fields.1  '    Z 

The  usual  approach  to  the  problem  of  increasing  radar  range  is  to 
look  for  ways  of  varying   some  of  the  parameters  in  the  radar  range 
equation  which  will  lead  to  a  greater  maximum  range. 

The  radar  range  equation  is 
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where     P    =  peak  transmitted  power  in  watts, 

P  =  minimum  peak  detectable  signal  in  watts, 

mm  r  b 

S  =  scattering  cross  section  of  target  in  units 
consistent  with  range, 

G  =  gain  of  transmitting  antenna, 

A     =  effective  area  of  receiving  antenna  in 
consistent  units, 


D  =  a  composite  loss  factor  for  transmission  lines, 
atmospheric  losses,   etc. 

The  parameters  under  the  control  of  the  design  engineer  are  A    ,G, 

P     .     and  Pi#     A       and    G    are  fixed  by  various  considerations   such  as 
mm  t  e  ' 

antenna  beamwidth,    side  lobe   suppression,   physical  size  of  the  antenna, 

etc.  ,   and  are  not  truly  available  as  adjustable  parameters  beyond  small 

variations.      Therefore  the  two  parameters  for  the  designers  to  optimize 

are     P  and    P.. 

min  t 

P     .       is  basically  limited  by  three  factors:     Johnson  noise  in  the 
min  ' 

input  circuits,    shot  effect  and  other  noises  in  the  first  tube,   and  cosmic 

noise  picked  up  by  the  antenna.      P  approaches     KBT     in  a  perfect 

mm  a 

receiver,   where     K    is  Boltzmann's  constant,    B    is  the  bandwidth  and  T 

a 

is  the  absolute   space  temperature.      P  is  currently  being  reduced 

toward  the  above  value  through  the  utilization  of  parametric  amplifiers, 

masers  and  other  low  noise  devices.      Another  method  of  reducing     P 

°        mm 

is  through  pulse  integration,    since  noise  voltage  averaged  over  a  period 
of  time  has  a  value  approaching  zero  while  the   signal  values  are  additive. 
However,   in  most  cases  this  integration  of  pulses  can  also  be  considered 
as  an  increase  in  transmitted  energy.      It  will  be  considered  from  that 
viewpoint  in  this  paper.      There  are  two  basic  types  of  integration,   co- 
herent and  noncoherent.      Coherent  integration  gives  a  gain  of  approxi- 
mately   N,    where     N    is  the  number  of  pulses  integrated  and  noncoherent 

integration  gives  a  gain  of  approximately  */N.      Both  of  these  values  are 

3 

compared  to  a  single  pulse    .      The  values  are  respectively  the  upper 

bound  for  coherent  integration  and  the  lower  bound  for  noncoherent 
integration. 

There  are  basically  three  different  integration  time  bases  possible 
in  a  radar  system.      One  is  the  antenna  scan  to  antenna  scan  where  the 


operator  or  a  computer,   decides  if  the  echoes  from  previous  scans  were 
signals  or  noise.      The  second  method  of  integration  is  the  interpulse 
method  where  the  cathode  ray  tube,    a  delay  line,   or  a  computer  sums  up 
the  hits  within  one  antenna  scan  of  the  target.      The  third  method  is  intra- 
pulse  integration,    or  pulse  compression,   where  different  portions  of 
one   radar  transmitted  pulse  are  summed  up.     Complementary  codes  are 
used  in  one  scheme  of  intrapulse  integration. 

An  examination  of  how  the  transmitted  power  or  energy  since 
E     =  P     A   T,     can  be  varied  is  now  in  order.      From  this   relation,    it  is 
seen  that  there  are  two  basic  methods  for  increasing  the  transmitted 
energy:    either  increase  peak  pulse  power,    or  increase  the  pulse  width. 
A  common  method  of  increasing  radiated  energy  has  been  to  increase 
the  peak  power.      This  method  has  proven  both  expensive  and  wasteful  as 
many  radar  transmitter  tubes  are  peak  power  limited  while  only  dissipa- 
ting a  very  small  percentage  of  their  average  power  capability,   and 
radar  modulators  are  taxed  with  a  very  similar  problem.      The  alterna- 
tive did  not  look  very  promising  at  first,    since  to  increase  pulse  duration 
normally  reduces  the  high  frequency  content  of  a  pulse,   which  in  turn 
reduces   range  resolution  capabilities.      However,    intrapulse  modulation 
retains  the  high  frequency  content  of  a  narrow  pulse,   while  increasing 
the  energy  content  to  that  of  a  pulse  of  long  duration.      Intrapulse  modula- 
tion allows  the   summing,    required  in  intrapulse  integration,   for 
detection  in  a  radar  receiver. 

There  are  two  methods  of  intrapulse  modulation,    the  digital  method 
and  the  analog  method.     Although  the  emphasis  in  this  paper  is  on 
digital  modulation  schemes   -  the  analog  system  will  now  be  briefly 

mentioned  since  it  is  a  part  of  the  overall  picture.      Figure     1.  1  contains 

4 
the  graphs  used  in  the  explanation  of  the  chirp  radar  scheme    .      In 
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Frequency  modulation  of  a  chirp  radar, 
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Signal  addition  through  the  use  of  a  delay  line. 
FIGURE     1.2 


this  analog  method  of  intrapulse  modulation  the  frequency  of  the  outgoing 

pulse  is  swept  linearly  upward  with  time  as  is  shown  in  Figure    1.  1-b. 
This  gives  a  large  high  frequency  content  to  the  pulse  while  at  the   same 
time  maintaining  a  long  pulse  duration.      Upon  reception  the   signal  is 
passed  through  a  reversed  time  delay,    delaying  the  low  frequencies  more 
than  the  highs,    as  is   shown  in  Figure     1.  l-d„     This  allows  all  the  fre- 
quency components  of  energy  to  arrive  at  the  output  simultaneously, 
giving  a  large  and  narrow  output  pulse.      The  type  of  filter  used  in  the 
delaying  process  is  known  as  a  matched  filter  and  is  discussed  in 
Appendix  I.      The  procedure  here  is   seen  to  be  quite  simple:    first  we 
generate  energy  components  with  various  frequencies  at  different  times, 
and  then  we  delay  all  components   so  that  they  arrive  at  the  output 
simultaneously. 

As  a  preliminary  to  the  digital  type  intrapulse  modulation  consider 
a  transmitted  waveform  such  as  in  Figure     1.  2.      This  waveform  con- 
sists of  three  pulses  of  unit  length  with  spacing  between  them  of  one  and 
three  time  units  respectively.      This  waveform  is  then  fed  into  a  delay 
line  in  the  receiver  which  has  the  delays   shown  in  Figure     1.2-b.      The 
output  is  either  a  one  or  a  zero  except  at  the  exact  match  of  signal  to 
filter;     the  output  is  then  three  as  is   shown  in  Figure     1.  2-c.      The  three 
pulses  have  therefore  been  summed  as  though  they  were  one,    at  the 
time  between     6    and     7    units.      These  pulses   still  contain  all  the  high 
frequency  components  of  the  short  pulse  of  unit  length  but  have  effec- 
tively three  times  the  energy  transmitted  in  a.  single  pulse  at  the  time 
of  exact  match.      There  are  two  features  to  be  noted  about  this  output, 
first  it  has  been  delayed  in  time  and  second  there  are   secondary  pulses 
at  various  undesired  times  which  might  be  classified  as  coding  noise. 
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Fig.   1.3-  Matched  filter  detection  of  a  random  code. 

Fig.    1.4  -  Matched  filter  detection  of  an  optimal  noise  code, 


I  will  consider  any  undesired  signal  which  might  be  mistaken  for  a  true 
signal  as  noise.      These  two  features  will  be  discussed  in  some  detail 
in  the  true  digital  intrapulse  modulation  schemes  to  be  discussed  later. 

The  scheme  just  considered  is  somewhat  inefficient  since   seven 
units  of  time  are   required  to  radiate  three  units  of  energy.      If  in  the 
transmitted  signal  all  the  spaces  with  zero  energy  output  in  the  pulse 
envelope  had  also  had  an  output,    seven  units  of  energy  would  have  been 
radiated  in  seven  units  of  time.     One  might  use   some   sort  of  frequency 
shift  or  phase  modulation  to  accomplish  this  purpose,    still  utilizing  of 
course  the  same  output  tube.     Consider  the  same  waveform  as  before 
with  the  units  of  energy  coded  as  plus  and  minus  by  the  modulator  by 
means  of  phase   shift  modulation  of  either  zero  or     180    .      The  waveform 
is  as  shown  in  Figure     1.  3-a.      Upon  reception  this  waveform  is  now  fed 
into  a  delay  line  in  the  receiver  as   shown  in  Figure     1  „  3-b.      "A"     repre- 
sents a  straight-through  amplifier,,      Figure     1.  3-c     shows  the  output  of 
the  delay  line  summer  with  respect  to  time.      The  output  energy  is  now 
seven  times  that  of  a  single  pulse  and  has  the  frequency  content,   and 
therefore  the  range  resolution,   of  one  of  the   single  pulses.      A  time  delay 
in  the  output  is  again  apparent,   as  are  a  number  of  noisy  sub- peaks. 
This  was  a  code  picked  at  random  and  the   results  were  good0      There  are 
certain  classes  of  codes  which  are  optimal,   having  even  smaller  noisy 
sub-peaks.      These  codes  are  of  length  2-1     and  have  been  studied 
extensively.     '     '        They  are  known  as     L    codes  or  pseudo-random  noise 
codes.     An  example  of  one  of  these  pseudo- random  noise  codes  is 
-  +  -   -  +  +  +,     where  the  minuses  and  pluses  indicate  the  phase  of  the  RF 
signal.      When  this  particular  noise  code  is  impressed  upon  its  matched 
filter,    Figure     1.  4-b,   the  output  of  this  summed  delay  line  is  as   shown 
in  Figure     1.  4-c.      It  is  to  be  noted  that  the  maximum  height  of  the 
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Matched  filter  detection  of  a  complementary 
code  of  length    8. 


FIGURE     1.5 


"hash"  (sub-peaks)  is  one,    with  a  main  peak  of  seven,   whereas  in  the 
previous  example  picked  at  random  the  hash  had  a  height  of  two. 
Elaspas  has  made  a  study  and  design  of  a   radar  utilizing  this  form  of 

o 

modulation  and  matched  filter  detection,,        He  infers  that  codes  of 
reasonable  length  should  be  used.      For  example,   a  code  of  length     127 
would  be  a  reasonable  length  to  obtain  the  benefits  of  a  pseudo-random 
matched  filter  radar.     An  optimum  code  of  this  length  has  a  maximum 
subpeak  of     13,    which  gives  a  good  coding   signal-to-noise  ratio. 

There  is  a  way,    however,   to  completely  eliminate  the  hash,   al- 
though this  required  the  transmission  of  two  separate  sequences.      Both 
of  these  sequences  could  be  transmitted  through  one  output  tube  with  a 

moderate  band  width.     These  codes  which  eliminate  the  hash  are  known 

9 
as  complementary  codes.        An  example  of  a  complementary  code  pair 

being  detected  by  a  matched  filter  pair  is   shown  in  Figure     1.  5.      The 
two  codes  being  transmitted  will  be  called    A    and    B    for  convenience, 
with  their  respective  matched  filters  also  being  labeled  in  a  similar 
fashion.      The  receiver  output,   after  the  complementary  codes  are  detec- 
ted by  the  matched  filter  and  their  outputs   summed,    is   shown  in  Figure 
1.5-g.      This  demonstrates,   as  was  mentioned  before,   that  the  hash  level 
is  zero.      This  lack  of  hash  is  indicative  of  an  infinite  coding  signal-to- 
noise  ratio.      This  particular  feature  is  the  prime  advantage  of  complemen- 
tary sequence  pairs  over  other  forms  of  noise  modulation.,      Figure   1.  5-g 

shows  the  same  amount  of  delay  as  did  all  the  other  forms  of  matched 
filter  detectors,   and  shows  that  the   sum  of  all  the  energy  sub-pulses  is 
incorporated  into  the  main  pulse  at  exact  match. 

Each  of  the  systems  described  in  this  introduction  has  had  the  char- 
acteristic that  the  transmitted  signal  is  spread  over  a  wide  frequency 
spectrum.     This  has  in  itself  two  advantages:     One  is  that  enemy  detection 
of  the   signal  becomes  more  difficult,   and  the  other  is  that  even  after 
detection  the  jamming  of  a  wide   spectrum  noise- like  signal  becomes 

10 


extremely  difficult.     Since  the  characteristics  just  described  are  de- 
sirable for  a  secure  communication  system  as  well  as  for  a   radar 
system,   all  of  the  advantages  that  accrue  for  the  complementary  se- 
quence modulation  in  radar  would  apply  in  a  communications  scheme. 

The  initial  major  objectives  of  this  investigation  were  two  in 
number. 

1.  The  formulation  of  a  set  of  operations  on  complementary 

sequences  which  form  a  group,   along  with  a.  theoretical 
study  of  the  invariant  properties  of  such  a  group. 
2C      The  search  for  new  complementary  codes.     An  exhaustive 
computer  search  for  a  new  code  of  length     26    was  made 
as  well  as  a  partial  search  for  codes  of  length    34. 
During  the  investigation  two  more  major  objectives  were  added 
to  the  work. 
These  were : 

3.  The  application  of  Hamming  distance  to  both  the  comple- 
mentary sequences  and  their  group  formulation,, 

4.  The  application  of  Hamming  vectors  to  complementary 

sequence  pairs  and  also  to  their  group  formulation,, 
The  original  objectives  were  oriented  toward  the  generation  of 
new  codes,   while  the  last  two  objectives  were  oriented  toward  the 
decomposition  of  codes0     All  of  these  objectives  were  attained  and 
consititute  a  large  portion  of  the  rest  of  this  workc      The  next  chapter 
is  a  study  of  complementary  sequences  designed  to  bring  the  reader 

to  the   state  of  the  art,   and  contains   some  proofs  which  have  not  ap- 

1    2    9    11    12 
peared  elsewhere.     '     '     '  Chapter     3     is  the  formulation  of  the 

operations  group  on  complementary  sequences  and,   along  with  several 

theorems  and  proofs  it  also  contains  the  identification  of  this  operations 

11 


group  with  a  known  group.      Chapter  3    as  written  requires  only  a  min- 
imum of  group  theory  knowledge.     Chapters    4    and     5    have  to  do  with 
the  Hamming  distance  and  Hamming  vectors  of  complementary  sequen- 
ces both  in  standard  form  and  in  group  form.     Chapter     6    is  a   short 
chapter  and  is  concerned  with  the  proof  of  two  theorems  used  in  the 
searching  for  new  codes.      Chapters     7    a.nd     8     describe  the  computer 
search  for  new  kernels  of  length     26    and     34     respectively,    the  actual 
computer  programs  used  being  listed  in  the  appendicie  s„      The  last 
chapter,    Number  9,   contains  the   summary,   conclusions,   and  suggestions 
for  future  research  in  this  interesting  field0 
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CHAPTER    II 
COMPLEMENTARY    SEQUENCES 

This  chapter  is  presented  to  provide  a  background  in  complemen- 
tary sequences  to  enable  the  reader  to  better  understand  the  material 
in  later  chapters.     Most  of  the  material  in  this  chapter  is  an  adaptation 

of  Golay's     "Complementary  Series"    with  expansions  and  deletions  to 

9 
fit  the  needs  of  this  paper.        The  formalization  of  some  proofs  and  the 

addition  of  new  proofs  given  here  have  not  appeared  in  print  elsewhere; 
however,  it  is  obvious  from  the  tone  of  Golay's  work  that  he  was  aware 
of  these  proofs.  This  modification  of  Golay's  work  is  offered  for  com- 
pleteness of  the  dissertation  and  not  as  new  work. 

In  the  introduction  a  pair  of  codes  of  length    n    were  discussed  which 
had  the  characteristic  that  when  detected  by  matched  filters,   the  sum  of 
the  two  filter  outputs  was  everywhere  zero,   except  where  the  signal 
patterns  had  zero  time  delay  with  reference  to  the  filter  patterns.     At 
the  zero  pattern  delay  time  the  output  was     2n,   leading  to  an  infinite 
coding   signal-to-noise  ratio.     A  study  of  the  necessary  and  sufficient 
conditions  for  this  to  be  true  will  now  be  considered. 

Assume  a  complementary  code  pair  (A,  B)  and  suppose  that  one  of 

the  pair,    i.  e.    A,   is  longer  than  the  other.     If  this  be  true  then  there  will 

be  an  output  from  the    A    filter  with  no  cancelling  output  from  the     B 

filter.     This  is  not  allowed;     therefore  the    A    and    B    codes  must  be  of 

equal  length.      Therefore  assume  a  code     A,   n  bits  in  length,   with 

elements    a,   through    a    ,     where  each    a      is  either    +1     or     -1.      Also, 
1  to         n  l 

assume  a  second  code     B,   n  bits  in  length,   with  elements    b.   similarly 
defined.     Detect  each  of  these  codes  with  a  matched  filter  (the  time 
inverse  of  the  code)  and  then  sum  the  outputs  of  the    A    filter  and  the 
B    filter.     At  time     1   (see  Figure   2.  1)  the   signals  a^    and    b-,    have 

entered  the  filter  segments     a       and    b       respectively  0     Since  at  this 

°  n  n  c  J 
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Time  1         a   a0a0a, 
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i  i 
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b3b2bl 


signal 


n       n-1        n-2.  »  .  „  „  0  b       filter 


Time  2 
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n     n-1      n-2 
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an+l-r    an+l-(r+l) 


Matched  filter  detection  of  a  complementary 
pair  of  length    n0 

FIGURE     2„  1 
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time  the  delay  between  signal  pattern  and  filter  pattern  is  not  zero,   the 

summed  output  voltage  must  be  zero.      This  gives  the  condition    a?a     + 

b,b     =  0.     At  time     2  (Figure  2.  1)  the  signals  are  one   segment  farther 

into  the  filter  and  the  output  must  again  be  zero,   giving    a,  a       ,     + 

a0a     +  b.  b      .   +  bob     =  0o 
L    n  1     n-  1  £   n 

At  time     3:    a,  a       ,+a-,a      ,+a0a   +b,b      T+b-.b       ,+b0b     =  0 
1     n-2      2n-l       3   n      ln-2      2n-l       3n 

At  time     r     (r  4  n) :    a,  a    .  .       +a->a    .  .     .      .  .  ,  , ,     ,  . 

x  '         1    n+l-r      2    n+l-(r-l)    ,  .  . .  +a   a   +b,b    ,,       +„  „  . 

r   n      1    n+l-r 

+  b   b     =0 
r   n 

At  time    n:    a,a,    +a->a0+....,)<)+a   a     +b!b,   +Oo  =  .o.+  b   b     =2n 
1122  nnll  nn 

for  at  time     n    the  pattern  of  the  signal  must  exactly  match  the  pattern 

of  the  filter.     The  explicit  statement  of  the  conditions  given  above  for  all 

times  is 

i=n-j  i=n-j 

=  0  j  4  0 


i=n-j  i=n-j 


J 


i=l         i    i+J  T     i=l        "^i+J 


=  2n        j  s  0 


(2.1) 


An  examination  of  the  equation  at  time     1  (j=  n  -    1)  shows  the  following 
eight  possible  solutions: 


1 


n 


1 


n 


1 
1 

-  1 

-  1 

-  1 
1 

-  1 
1 


These  possibilities  show  that  if  the    A    pair  are  alike  3   the     B    pair  must 
be  unlike,   and  if  the    A    pair  are  unlike,    the     B    pair  must  be  alike.      The 
equation  at  time     2(j=n-2)  indicates  that  if  the  number  of  likes  (like  pairs, 
i.e.    a    and    a       ,  ,)    of  the    A    code  and  the  number  of  unlikes  (unlike 
pairs)  of  the    B    code  are  equal,    and  if  the  number  of  unlikes  of    A 
equals  the  number  of  likes  of    B,    the  equation  is  satisfied.     Extending 
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this  reasoning  to  time     r     shows  again  that  if  the  likes  of    A    equal  the 

unlike s  of    B    and  if  the  unlike s  of    A    equal  the  likes  of    B    the  equation 

is  again  satisfied., 

It  will  now  be  proved  that  the  last  half  of  the  previous  statement  is 

not  necessary,    since  the  first  half  of  the   statement  contains  the  necessary 

and  sufficient  requirement. 

Let    U     =  number  of  unlike  pairs  in    A    for  the   spacing  specified, 
(j    in  equation    2.  1) 

Let     L     =  number  of  like  pairs  in    A    for  the  spacing   specified,, 
(j    in  equation    2.  1) 

Let    U,    =  number  of  unlike  pairs  in    B    for  the   spacing   specified, 
(j    in  equation    2„  1) 

Let     L      =  number  of  like  pairs  in     B    for  the   spacing  specified, 
(j    in  equation    2"„  1) 

Assume  that     L     =  U,  ,     we  want  to  prove  that    U       =     L, 
a  b'  ^  a  T)" 

Since  the   sum  of  the  number  of  like  pairs  at  each  spacing  with  the  num- 
ber of  unlike  pairs  at  the  same  spacing  must  equal  the  total  number  of 
possibilities , 

U     +  L     =  n-(n-r)     for     r  4  n    where     n    is  the  length  of  the  code 

and    (n-r)    is  the  number  of  units  of  time  delay 
from  exact  match,, 

For  the  same  reason 

Ub  +  H)  ~  n-(n"r)» 

But  since     U,    =   L    , 
b  a 

Ua  +  Ub  =  n-(n-r). 

Therefore     U     =  L,  . 
a  b 

The  definition  of  a  complementary  pair  of  sequences  will  now  be 
given.     A  pair  of  binary  sequences  of  equal  length  with  the  number  of 
like  pairs  of  one  sequence  equal  to  the  number  of  unlike  pairs  of  the 

other  sequence  for  each  possible   spacing  is   said  to  be  a  complementary 

code. 
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Two  schemes  of  representation  will  be  used  in  this  pape  r„      These 
schemes  will  be  used  interchangeably  at  the  convenience  of  the  author. 
The  first  representation  is  the  one  already  presented  where  each    a. 
and    b.     is  either    +1     or     -1„     The  operation  used  in  this  case  is 
ordinary  multiplication.     The  second  representation  will  use  the  opera- 
tion of  modulo  two  addition  denoted  by    ®,     where  the  elements    a.     and 
b.     are  either    0     or     1„      That,  these  two  groups  (1,  -  1;«  )  and  (0,1;  9  ) 
are  insomorphic  is  shown  in  the  following  tables. 


• 

1 

-1 

1 

1 

-1 

-1 

-1 

1 

9 

0 

1 

0 

0 

1 

1 

1 

0 

whe  re     1       -»        0 
-1       -        1 

The  first  form  will  normally  be  used  when  correlating  of  the  codes 
with  their  matched  filter  is  under  discussion  while  the  second  form  will 
normally  be  used  when  discussing  the  intrinsic  properties  of  the  indi- 
vidual codes. 

The  necessary  and  sufficient  condition  for  a  pair  of  sequences  to  be 

complementary,    is  in  terms  of  the  modulo  two  representation 
i=n-j  i=n-j 

F.=  E        (a.«a.,.)    =       £        (b  «b      .91)    for  all  j,    l<j<n-l       (2.2) 

J  i  _.  j  i+J  i  =  i       x      x_rJ 

or  an  equivalent  formulation  by  change  of  summation  index 

i=j  i=j 

V  y 

f .  =        .    ,    (a  9  a        L.)    =    .    .(b.«b  .0  1)  (2.3) 

j  i  =  l   *   i        n-j+i'  i=l    x    l        n-j+i  '  \    •-'/ 


The  above  £     is  standard  summation,      If  modulo  two  summation  is 

desired  the  symbol    2    will  be  utilized.      (The  equivalent  conditions  are 

given  since  they  are  both  used  in  formal  proofs  in  a  later  portion  of 

the  paper.     Due  to  the  possibility  of  confusion  it  seems  more  appropriate 
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to  introduce  the  second  form  now  rather  than  later  when  needed.) 

U     =  L,      has  been  used  to  represent  these   same  equations,    and  will  be 
a  b 

used  elsewhere  in  the  paper  when  it  is  more  convenient  to  use  it. 

An  interesting  property  of  any  system  of  codes  is  the  possible  in- 
variances  it  may  have  under  various  types  of  transformations.      It 
therefore   seems  appropriate  to  consider  the  possible  transformations 
on  a  complementary  sequence  pair  which  leave  the  pair  complementary 
and  of  the  same  length  after  the  transformation.     A  necessary  tool  for 
this  study  is  the  property  that    a. fa.    .  =  a.  0  a- ,  •    where  the   symbol 

a.     means  the  complement  of    a.,    0  =  1»    1=0. 

1^1 

Theorem  2.  1 

The  modulo  two  sum  of  a  pair  of  binary  numbers  is  equal  to 

the  modulo  two  sum  of  their  complements. 

The  truth  of  the  theorem  is  obvious  from  the  following  exhaustive 

table  of  four  possibilities: 

a .  w  a  . .  • 
i         l+j 

101  =  000  =  0 

10.0  =  0*1  =  1 

001  =  10      0  =  1 

000  =  101  =  0 

Therefore     a.  ta.  ,  .  =  a,  $  a.  ,  .  » 

i       i+j         i         i+j 

*  *  * 

In  each  of  the  following  proofs  it  will  be  considered  that    A    and    B 

are  a  complementary  sequence  pair  of  length    n.      Complementing  each 

element  in  a  code  is  called  complementing  the  code. 

Theorem  2.  2 

Complementing  the     A    code,   or  the     B     code,    or  both  the 

A    and    B     codes,    results  in  a  pair  of  complementary  codes.      (The 

proof  is  given  just  for  the    A    code  but  is  identical  for  the     B     code  and 
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a.  0  a. .. 
i         i+j 

1      0       1 

= 

0 

1       0       0 

= 

1 

0      0       1 

= 

i 
i 

0      0      0 

- 

0 

from  both  of  these  for  the    A    and    B    codes.  ) 

By  hypothesis 

i=n-j  i=n-j 

F.=       S      (a.tai+j)     =      S  (b.  %  b.+j  t  1). 

i=  1  i=l 

Complementing  the     A    code  transforms 

i         i 

but    by  Theorem    2„  1 

a.  $  a  ,  =  a.  •  a.  .  .  . 

i         i+j  i         i+j 

Upon  substitution 

i=n-j  i=n-j 

F.=       T.       (a    ©i^.)    =     S         (b.  •b.,.«l). 
J        i=1  i         i+J  i=1  *         1+J 

Therefore  complementing  the     A    code  has  no  effect  on  complementarity, 

and  similarly  for  the     B    code  or  for  both  codes. 

*  *  * 

Interchanging  the  first  and  last  bits  of  a  code,  the  second  and  next 
to  last  bits,  the  third  and  third  from  last  bits,  and  so  on  is  called  time 
inversing  or  time   reversing  a  code. 

Theorem  2„  3 

Time  inversing  the     A    code,   or  the    B     code,    or  both  the    A 

and    B     codes*    results  in  a  complementary  pair„      (This  proof  is  again 

just  for  the    A    code,    but  also  applies  to  the    B     code  and  therefore  to 

both  the     A    and    B     codes,,  ) 

By  hypothesis 

i=n-j  i=n-j 

f.  =       S       (a.  $  a       )      =        I        (b    «  bA  .  §  1). 
J         i=i         *         1+J  i=i  x         1+J 

The  time  inverse  of    A    causes  each 

a.    -    a    ,  ,     .  . 
l  n+1-3 

Applying  this  time  reversal  just  to  the     A    code  shows  that 
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i=n-j  i=n-j 

F.    =     S       (a.  ©  a,  .)    -     S        (a    , .    .  «  a    .  .    .....  ). 

j  j         i         i+j#  .  n+l-i         n+l-(i+j) 

Expanding  both  expressions  for  a  few  terms  and  remembering  that 

modulo     2    addition  is  commutative  as  is  regular  addition  we  see  that 

(on  the  left)  (on  the  right) 

for    j=n-  1 ,     F      .   =  a,   $  a  F      ,    =  a     ®  a, 

J  n-lln  n-lnl 

for  J=n-2,    F      ->  =  a,  0  a      ,    +  a->§  a       F      -,  =  a    ®  a ->  +  a      ,0a,. 
n-  l  1        n-1  2        n        n- 2         n        L         n- 1  1 

These  expressions  indicate  that  by  a  change  in  summation  index  the  two 
forms  of    F.    are  equal,    as  is  the  case.      Therefore  time   reversing  the 
A    code  has  no  effect  on  the  complementary  property,    and  similarly  for 
the     B     code  and  for  both  the    A    and    B     codes. 


The  altering  of  a  code  is  a  transformation  where  every  other  bit  of 
both  codes  is  complemented.     In  a  different  portion  of  the  paper  a  dis- 
tinction is  made  between  altering  odd  bits  and  altering  even  bits„      In  the 
following  proof,    however,   just  the  even  bits  are  altered,    although  a 
slightly  modified  proof  would  hold  for  altering  the  odd  bits0 

Theorem  2„  4 

The  result  of  altering  a  complementary  pair  of  sequences  is 

again  a  complementary  sequence  pair. 

Since  the  codes  are  a  complementary  pair, 
i=n-j  i=n-j 

F.  =       T.        (a.  0  a.x.)        =       £         (b.  ©  b...  ©  1). 
J         i=1  i         i+j'  i=1         i         i+j 

Divide  the     F.    into  two  groups,    F,      for    j=     an  even  number  and    F 
j  k  r 

for    j=    an  odd  number,,      The  portion  of  the  above  equations  for    j=  an 

even  number  is 

i=n- k  i=n-k 

Fk=    2,  <ai»ai+k>   =    £,    <bi»bi+k,1>- 

i=l  i=l 
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For  the  theorem  to  be  true,   when    j    is  an  even  number  the  above   set 
of  equations     F,      must  be  satisfied  for  all    k    when  the  pair  is  altered,, 
(The  proof  is  given  only  in  terms  of    A,   but  duplicate   steps  must  also 
be  applied  to    B.  ) 

When    i=an    even  number,   then     i+k=an    even  number.      In  the  alter- 
ing process   since  all  even  bits  are  complemented    a.  ®  a.  .,    -»  a.#a.  , ,    , 
°  c  1  l+k         1      l+k 

but  by  Theorem    2.  1     these  are  equal.      When    i=an    odd  number,   then 
i+k  =  an  odd  number  and  these  are  not  affected  by  the  altering  process, 
since  only  even  bits  are  complemented.      Therefore     F,      is  invariant 
under  altering,    or  symbolically    F,    =  F,      ,,. 

Examination  of    F       where     r    is  an  odd  value  of    j,    shows  that  of 

the  pair    a.,    a.  .        one  of  the  bits  is  complemented  and  the  other  is  un- 

r  i        l+r  r 

changed. 


a.  la.,  (a.®  a.  .         or    a.©  a.  .    )    =    a.  <&  a.  . 

i         l+r  -*       l        l+r  l        l+r  l         l+r 

Therefore     F        ,.     =  (n-l)-F       for  both  the    A    code  and  the     B     code. 

r  alt         l         '       x 

This  signifies  the  likes  of  the     A    code  for  odd  spacings  are  changed  into 
unlike s,   and  the  unlikes  of  the    A    code  for  odd  spacings  are  changed  into 
likes;     however,   for  each  change  in    A    there  is  an  opposite  change  in 
B     so  the  total  likes  in    A    still  equals  the  total  unlikes  in    B    for  each 
spacing.      It  is  to  be  noted  that  altering  both  the  even  and  the  odd  bits  is 
the   same  as  complementing  both  codes. 

*  *  * 

The  operation  of  interchanging  the     A    and    B     codes,   although 
trivial,    is  also  a  transformation.     There  are  a  total  of  sixty  four  pos- 
sible transformations.      These  can  all  be  generated  by  combining  the 
above  listed  transformations.,      They  will  be  discussed  much  more  thor- 
oughly in  Chapter  3. 
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Whenever  a  new  facet  of  a  pair  of  sequences  such  as  the  comple- 
mentary property  is  discovered,    it  is  natural  to  wonder  just  how  general 
this  property  might  be.      One  question  of  interest  would  be  the  possible 
length  limitations  of  the  code  pairs  that  might  be  forced  by  the  comple- 
mentary property.      The  following  proof  due  to  Golay  shows  that 

2      2 
complementary  sequences  must  have  a  length    n    of  the  form  n=x  +y    , 

where     n    is  an  even  integer  and    x    and    y    are  integers.      The  first 

theorem  to  be  proved  demonstrates  that    n    must  be  an  even  integer 

2      2 
and  the  second  theorem  shows  the     x  +y       form. 

Theorem    2.  5 

A  necessary  condition  for  a  sequence  pair  to  be  complementary 

is  that  their  length  be  an  even  number. 

By  hypothesis 

i=n-j  i=n-j 

F.  =        S     a. 0a ...       =         S       b  0b ...  0  1    . 
J         i=1       i      i+J  i=1  i      i+J 

Let    G.  =  F.  modulo     2,    so  that 

J  J       i=n-j  i=n-j 

G.  =       B      a.®  a. .        =        E        b.  0b.  ,.  0  1  . 

j        i=1        i       i+j  i=l         i       i+j 

An  expansion  of    G    for  a  few  different  values  of    j    will  suggest  the 
proof.     Start  with    j=n-l„ 

n-1  In  In 

G      ,  =a,0a      ,0a, 0a     =  b .  0b      ,010b.,0b    01, 
n-2  1      n-1       2      n  1       n- 1  2      n 

G      ,  =a,la.      ,0a       ,9a      ,0a, 0a    =b,0b      ,01   0b      .010b, 0b     9  1  „ 
n-3        1      n-2      n-1      n-1       3      n      1       n-2  n-1  3      n 

Change  the  order  of  the   sums  in    G      2     an<i    G       ,    to  achieve  the  following: 

2 


G      ,=a,  0a, 0a       ,0a  =b,0b,0b       ,0b     9     g    1, 
n-2       1        2      n-1      n      1.      2      n-1      n 


I 


G        =a,0a,0a,0a      ,0a      ,0a    =b ,0b, 0b ,0b       ,0b      ,0b   0      g  I  . 
n-3       1      2      3      n-2      n-1      n      1       d,      J      n- c      n~  1       n 

In  general  for    j=n-r 

G        =a,0a,0.  ..0a   0a      -     ,nte      /  ,4.?\*- •  •  *a„    ,0a   =b, 0b?0. .  . 
n-r       12  r      n-fr+1)      n-(r+Z)  n-1       n      i       l 
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®b   ©b      IA.,9b      ,   J_9V®.  „.«b      ,©b  ©     E     I  , 
r      n-(4  +  l       n-(r+2)  n- 1       n         \  ' 


and  for     j=n-(r  +  l) 


Gn-(r  +  l)=al§a2§-^ar§ar  +  l®an-r§an-(r+l)^-§an-iean=blfV--- 

r+1 

.  .»b    ©b    ,,©b        ©b      .        .©b      ,    .,v*„...»      S      1   . 
r      r+1      n-r      n-(r+z)      n-(r+2)  f 

Forming  the   sum    G         §    G      .    ,  . .    will  leave  only  five  terms  in  the 

n-r  n-(r+l)  7 

equation  since    G  and    G       ,     ,  . .     differ  only  in  a         ,    a    ,  ,  ,    b    ,  ,  , 

n-r  n-(r+l)  '  n-r        r+1'      r+1' 

and    a,   bcause 

r  r+1 

Sli        E     1  =  1  . 
1  1 

G  ©G      /    .    »  =  a    ,,»a         =  b    ,  ,  ©b        ©    1  . 

n-r         n-(r+a)  r+1      n-r  r+1      n-r 

Adding    b   +  1©b  modulo     2     to  both  sides  of  the  equation  gives 

a    ,  ,0a        ©b    ,  ,  ©b        =  *1.  (2.4) 

r+1      n-r      r+1       n-r  ' 

Let    n=2s-l     and    r=s-l;     then  to  satisfy  the  equation  just  derived 

a   ©a,       .        ..   ©b     ©b7      ,        ,  ,    =  a   ©a   ©b   ©b   =  1  ,  which 
s       Zs-l-s+1  S         £S-1-S+1  s        s        s       s 

is  obviously  false. 

Therefore     n    cannot  have  the  form    2s- 1     and  must  be  an  even  number. 

*  #  # 

Equation    2,,  4    which  was  a  step  in  the  previous  proof,   is  a  very 
important  necessary  condition  for  a  pair  of  codes  to  be  complementary 

and  will  be  referred  to  in  the  future  as  the  parity  check.      The  next 

2      2 
theorem  which  shows  the     n=x   +y     form,    starts  with  the  basic  assump- 
tion that    A    and    B     are  a  complementary  pair  of  length    n. 

Theorem    2C  6 

A  necessary  condition  for  a  pair  of  codes  to  be  complementary 

is  that  their  length  be  the  sum  of  the  squares  of  two  integers., 

By  hypothesis 

i=n-j  i=n-j 

F    =      2    (a    ©a       ©1)    =        S       (b  ©b       )  . 

J  i=1  1  l+J  i=1  1        ""J 

Assume  that  the     A    code  satisfies     F.     and  that  the    A    sequence  contains 

J 
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p    ones  and    n-p    zeros.      The     A    code  has  a  weight  of  p.      Similarly 
assume  that  the    B     code  satisfies     F.    and  that  its  sequence  contains 
q    ones  and    n-q    zeros.      The     B    code  has  a  weight  of    q„      For  the  total 
likes  in  the    A    code,    each  bit  is  matched  against  each  other  bit  exactly 
once.      For  the  total  unlike  s  in  the     B     code  each  bit  is  again  matched 
against  each  other  bit  once.      There  are  two  possibilities  of  likes  in  the 
A    code,      1     matched  against    1,    and    0     matched  against    0o     The  total 
number  of  like  pairs  for  the     A    code  is  therefore  the  number  of  com- 
binations of    p    ones    taken    2     at  a  time  plus  the  number  of  combinations 
of  n-p    zeros  taken    2    at  a  time,   for  the  unlike  pairs  in  the     B     code, 
q    ones     are  matched  against    n-q  zeros. 
Thi  s  gives  the  equation 

PlP-J!     +    (n-£)(n-p-l)    =    q(n_q) 

expanding  and  simplifying 

2  2         2 

n=n   -2qn+2q   +  2p    -2np, 

adding  and  subtracting  2pq  gives 

2        2    2 
n=n  -2qn+2q  +2p  - 2np+2qp- 2qp, 

collecting  terms  gives 

n=  n    +q    +p     +2qp-2np-2nq  +  p     -2pq  +q 
combining  terms  gives 

n=  (n-p-q)      +  (p-q)  (2. 5) 

Therefore  complementary  sequence  lengths  are  permissible  only  in 

lengths  which  are  even  numbers  and  formed  by  the  sum  of  two  squared 

integers. 

*  *  * 

A  list  of  all  possible  code  lengths  up  to     200    with  the  number  of 
ones  allowed  in  the     A    and    B     codes  is  given  in  Table     2.  lu 
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*  Indicates  possible  kernel 


Code 

Lengths 

*2 

4 

8 
*10 

16 
*18 

20 
*26 

32 

*34 

36 
40 

*50 

52 

*58 

64 
68 

72 

*74 

80 

*82 

*90 

*98 
100 


Unordered 
Weights  of(A,B) 

2,1)(1,0) 

3,3)(1,1)(3,1) 
6,4)(4,2) 

7,6)(4,3)(7,4)(6,3) 
10,  10)(6,6){10,6) 
12,9)(9,6) 

,  11)(9,7)(11,  7) 


16 
15 

20 

21 

21 

21 

24 
24 

30 
29 

31 
25 

34 
34 

36 

39 
39 

42 

43 
43 

46 
46 

46 
46 

51 
42 

56 

55 
57 
49 


15 
10 

16 

18 

16 

21 

22 
18 

25 
22 

27 
21 

31 

27 

28 

37 
31 

36 


36 

42 
38 

45 
37 

48 
39 

49 

55 
51 
43 


(11,10) 
(16,11) 

(16,12) 

(16,13) 
(18,13) 

(15,15)(21,15) 

(18,16) 
(22,16) 

(25,20)(29,28) 
(28,21)(22,21) 

(31,25) 
(27,21) 

(27,24) 
(31,24) 

(36,36)(28,28) 

(31,29) 
(37,29) 

(36,30) 


38)(36,31) 


(38,31) 

(38,34) 
(42,34) 

(37,36) 
(45,36) 

((51,42) 
(48,39) 

(49,42) 

(55,45)(45,45) 
(57,49)(51,43) 


Code 

Lengths 

104 
*106 

116 
*122 

128 
*130 

136 

144 
*146 

148 

160 
*162 

164 
*170 

*178 

180 

*194 

196 
200 


Unordered 
Weights  of  (A,  B) 

(58,56)(58,48)(48,46)(56,46) 

(60,55)(51,46)(60,51)(55,46) 

(65,61)(55,51)(65,55)(61,51) 

(67,66)(56,55)(67,56)(66,55) 

(72,64)(64,56) 

(72,69)(72,61)(61,58)(64,57) 
(73,66)(73,64)(64,57)(66,57) 

(76,70)(66,60)(76,66)(70,60) 

(78,78)(66,66)(78,66) 

(81,  76)(70,65)(81,70)(76,65) 

(81,79)(69,67)(81,69)(79,67) 

(88,84)(76,72)(88,76)(84,72) 

(90,81)(81,72) 

(91,83)(81,73)(91,81)(83,73) 

(94,87)(83,76)(94,83)(87,76) 
(92,91)(79,78)(92,79)(91,78) 

(97,94)(97,84)(84,81)(94,81) 
(99,93)(87,81)(99,87)(93,81) 

(106,101)(93,88) 

(106,93)(101,88) 

(105,  105)(91,91)(105,91) 

(110,  100)(108,  106)(100,90) 
(94,92)(94,108)(106,92) 


All  possible  complementary 
code  lengths  up  to  200  with 
unordered  possible  weights. 


TABLE     20  1 
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The  shortest  possible  complementary  pair  is    A=ll,    B=10„      This  pair  or 
any  of  its  transformations  is  called  a  kernel  of  length  two,   or  the  quad. 
A  kernel  is  a  basic  length  code  which  cannot  be  decomposed  into  shorter 
length  codes  by  an  inversion  of  the  standard  generating  methods  to  be 
explained  later.      Some  possible  kernel  lengths  are  2,    10,    18,    26,    34, 
50,    etc.     although  codes  for  all  of  these  do  not  exist.      Complementary 
sequences  which  are  not  kernels  are  called  composite  complementary 
sequences. 

The  above  list  of  kernel  lengths  did  not  include    n=4    or    n  =  8 
which  are  possible  complementary  sequence  lengths,,      These  if  they 
exist  must  therefore  be  composite.      All  possible  codes  of  length  four 
can  be  generated  from  the  exhaustive  list  of  possibilities  for  four  binary 
digits  given  in  Table     2„  2. 

0000  1000 

0001  1001 

0010  1010 

0011  1011 

0100  1100 

0101  1101 

0110  1110 

0111  1111 

Table     2.2.     All  Possible  Binary  Numbers  of  Length  Four. 

The  unordered  possible  pairs  of  ones  in  the  codes  are  (1,1),    (3,3) 

and  (3,  1).      Thus  all  possible  codes  of  length  four  have  either  one     1     or 

three     l's    because  of  this  limitation.      This  reduces  the  table  to  the 

eight  following  numbers  of  length  four: 

0001  0111 

0010  1101 

0100  1011 

1000  1110 

The  possible  unordered  (1,  1)  code  pairs  are  the  following  four: 

1000  1000  0001  0001 

0100  0010  0100  0010 
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The  possible  unordered     (3,  1)    pairs  of  codes  are  the  following  eight: 

0111  0111  1110  1110  1011  1011  1101  1101 

0010  0100  0100  0010  0001  1000  0001  1000 

The  possible  unordered     (3,3)  code  pairs  are 

0111  0111  1110  1110 

1011  1101  1011  1101 

Comparing  the  last  code  to  the  kernel  of  length  two,   A=ll,   B=10, 
shows  that  it  might  have  been  constructed  by  writing  in  time  sequence 

S1=AB=a1a2b]b2     =     1110, 

S2  =  AB=a1a2b162    =     1101. 
The  next  to  last  sequence  pair  might  have  been  formed  by  inter- 
lacing   A    and    B     in  the  following  manner: 

Ti  =  a1b1a2b2    =     1 110, 

Ts  =  a1^1a2^2    =     1011. 
AH  the  rest  of  the     16    pairs  could  be  considered  either  as  transforma- 
tions of  the  last  pair,   or  as  being  formed  by  the   same  two  composite 
generating  operations  used  above  on  the     8    transformations  of  the  ker- 
nel of  length  two.      A  general  proof  will  now  be  given  to  show  that  both 
the  time  sequence   scheme     (Si  »    S2)     and  the  interlace   scheme  (Ti,   T2) 
will  always  form  complementary  sequence  pairs  providing  that    A    and 
B    are  a  complementary  sequence  pair. 
Theorem    2„  7 

11    J\~  a-.a-ia-o.^.^eo-a      -.a 
1    Z    3  n- 1    n 

B=b, b0boe ...... b       ,b 

12    3  n- 1    n 

are  a  complementary  sequence  pair,   then 

C  =  a,  a-,.  .  .  a   b.b,,.,  ....    b 
1    Z  n    1    Z  n 

D  =  a,  a-,a~.  „  a  b,b0.  ,0.u.   b„ 
1    Z    3        n    1    Z  n 

are  a  complementary  sequence  pair. 
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Theorem    2„  8 

If    (A,   B)    are  a  complementary  sequence  pair,   then 

C  =  a,  b,  a0b->0  .  .  „  .  „  „  a   b 
1    1    2    2  n   n 

D  =  a1b1a2b2.e..o..anbn 

form  a  complementary  sequence  pair. 

The  following  notation  will  be  used  in  the  proof: 
i=n-j 
U       =    E       a.0a.,.  for  all    jl<  j<n-l 

L       =    1^,n_Ja.«a..  •  ©  1  for  all    ia<  j<n-l  8 

a  Si      i+j  J      —  j  - 

i=l 

It  is  understood  that  the  two  CD  sequences  in  the  two  theorems 

above  are  not  the  same  but  because  the  proofs  are  identical  two  sets  of 

symbols  will  not  be  used.     For.each  spacing     j    there  are  three  possible 

ways  for  bits  to  match,   within    A,   within    B,    or  from  a  bit  of    A    to  a 

bit  in    B,      These  possibilities  are  denoted  by    L       JL,      and     L.  ■,   or    U    , 

U,  ,    U    ,      respectively. 

The  necessary  and  sufficient  condition  for    C    and    D    to  be  a 

complementary  sequence  pair  is  given  by     L.     =  L.  +  L,  +L.  ,      and 

UJ  =  U    +UC  +  U   r     where     L     =  UJO 
d  a      b  ab  c  d 

Now     L.     =  U,      and     L,    -   U       ,  r.  f  ,  <. 

a  b  b  a    by  definition  ot  complementary. 

Adding 

a  o  a       b 

but  ,  t 

i=n-l  i=n-l_ 

Uu  =  Ur       since        S      b.«b.,        =        T.     b=  *  b  i+j 
b  b  i=1       i      i+J  i=1      ' 

so  substituting 

L     +   L,     =  U     +  Ur 
a  d  a  b 

All  bit  matches  from    A    to    B     in  the    C     sequence  are  paired  in  the     D 
sequence  with  bit  matches  from    A    to    B     therefore, 
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ab  ab 

and  adding 

L     +  L.+  L,   ,    =  U      +Uc  +  Ur  or     JL     =  U  , 

a  d  ab  a  b  ab  c  d 

*  *  * 

Golay  describes  two  other  methods  of  generating  composite  sequen- 
ces from  shorter  complementary  sequences.      These  methods  will  not  be 
used  too  extensively  in  the   rest  of  the  paper  and  will  therefore  just  be 
mentioned  rather  than  proved.     Given  two  complementary  sequences 
pairs  (A,  B)  (C,D),     A    of  length    m    and    C     of  length    n. 
Let 

Ux  =  ACxAC?.  ........  Acn    Bdl  Bd?.  ........  Bdn 


and 


i],  -    Adn  Adl       RCn  R^i 


V,  =  ACl  Bdl. .  .  .  ........  Acn  Bdn 

V^  -   Adn  u^n  a^i  RC! 


(2.6) 


(2.7) 


It  can  be   shown  that  the  pairs  of    U    are  complementary  as  are  the  pairs 
of    V.      Where,    if  an  exponent  is  a  one  the    A    or    B    code  is  left  un- 
changed,  and  if  the  exponent  is  a  zero  the    A    or    B     code  is  comple- 
mented.     The  lengths  of  codes  which  can  be  generated  from  these  two 
methods  are     2mn    where     m    is  the  length  of  the  (A,B)  sequence  pair 
and    n    is  the  length  of  the     (C  ,  D)  sequence  pair.      There  is  also  another 
special  method  of  generation  which  applies  only  to  codes  of  length    2    . 
This  method  is  explained  in  Chapter     5    where  it  is  necessary  for  the 
completeness  of  a  proof. 

The  kernel  of  length  2,   or  the  quad,    has  been  discussed  in  the 
previous  pages.      There  are  two  kernels  of  length     10.      These  are 
1001010001     and     1000000110,    and  a   second  kernel     0101000011     and 
0000100110.      The  next  possible  kernel  size  is     18    and  Golay  has  proved 
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by  exhaustive  search  that  no  kernel  of  this  length  exists.     Kruskal  has 

1  3 

since  then  completed  the  proof  analytically.  One  of  the  goals  of  the 

writer  in  his  study  was  to  make  an  exhaustive  search  for  codes  of  kernel 
length    26.      This  search  disclosed  there  was  only  one  kernel  of  length 
26.     Chapter     7    of  this  paper  describes  this  search  in  detail  and 
Chapter     8    describes  attempts  to  find  kernels  of  length    34. 
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CHAPTER    III 
THE    GROUP    OF    OPERATIONS 

A  large  portion  of  this  investigation  is  devoted  to  the  set  of   allow- 
able transformations  on  complementary  sequences.      These  allowable 
transformations  possess  the  characteristic  that  the  length  of  the  sequen- 
ces and  the  complementary  property  are  invariant  under  transformation. 
To  be  useful  in  searching  for  new  kernels  it  is  also  desirable     that  these 
transformations  have  the  following  features: 

1.  That  they  form  a  group  for  all  general  transformations. 
Where  a  general  transformation  is  one  which  can  be 
applied  to  any  complementary  code  pair. 

2.  The  transformations  as  represented  be  simple  manipula- 
tions on  the  elements  of  the  group. 

3.  The  supplementary  characteristic   (to  be  explained  in 
Chapter  6)  be  a  property  of  the  representation. 

4.  The  parity  check  be  easily  made. 

Some  of  these  transformation  operations  were  discussed,    and  the 
proof  of  invariance  given  for  all  group  generators,    in  the  background 
chapter,    Chapter  2,   on  complementary  sequences.     As  a  continuation, 
a  brief  review  of  the  transformation  operations  and  a  tagging  of  symbols 
to  these  operations  will  prove  useful.     Each  basic  operation  symbol  and 
its  definition  is  given  in  Table     3.  1„     More  definitions  are  given  than  are 
necessary  to  generate  the  group;     however,    the  redundance   seems  to  the 
writer  to  be  an  aid  to  understanding. 

A  check  through  Table     3.  1     shows  that  the  transformations  affect 

the  complementary  sequence  pair  in  four  different  ways.     Ci  ,    C2  , 

Tx  ,   T3    affect  all  the  bits  of  one  of  the  code  pairc     E,T,C    affect  all  bits 

of  both  codes,   while    Ai  and    A2  affect  half  of  the  bits  of  both  codes.      The 

identity,      1,    of  course  has  no  affect  on  either  codec     A  reasonable  way 

to  form  the  elements  of  the  group  to  account  for  the  above  affects  would 

be  to  divide  each  code  of  the  pair  into  two  units.     These  would  be     A 

(even  bits)  and    A     (odd  bits),    also    B     (even  bits)  and    B     (odd  bits), 
o  &  o 
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The  Identity 


a.  -»  a. 
1         1 


b.  ->b. 
1         i 


A  code  time   reversal 


ai   "an+l-i 


b.    -b. 
i         i 

B    code  time  reversal 


a.  -»   a. 
i  i 

b.  -  b    ,  .    . 
i  n+ 1  - 1 

Time  reversal  of  both  codes 


a.    -»  a    ,  , 
i  n+ 1  - i 


b.    -  b    ,,     ; 
i  n+ 1-1 


Ai  Alter  odd  bits 


a2i-l"a2i-i  a2i  "*   a2i 

2i-  1         2i-  1  2i  2i 


A2  Alter  even  bits 


a2i-T 

a2i-l 

a2i^ 

a->. 

2i 

b2i-l 

2i-  1 

b,. 

2i 

b,. 

2i 

Complement  the     A    code 


a.      - 
i 

a. 

i 

b.      - 
i 

b. 

i 

c 

omple- 

ment  the 

B    code 

a.      -» 

i 

a. 
i 

b.      - 

b. 

1  1 

Complement  both  codes 


a.      -»     a. 
l  l 


b.      -     b. 
i  i 


Some     operations. 


Exchange  the     A    code  with  the    B    code 

ai      -     bi  TABLE    3.  1 

b1      -     a1  32 

i  i 


o  1    3    5  n- 1 

A.     =  a->a  .a/.  ........  a. 

e         L    4    b  n 

B     =  b^-b^.  ........ b      . 

o  1    3    5  n-  1 

e  Z    4    o  n 

However,  this  particular  division  for  representing  the  complemen- 
tary sequences  proves  to  be  very  awkward  for  the  time  reversal  opera- 
tions (Ti  ,T2  ,    T)  and  for  the  ease  of  the  parity  check  condition. 

Fortunately,   one  simple  change,    that  of  time   reversing  the  second  and 

14 
fourth  row  (even  bits  of    A    and    B),    clears  up  these  difficulties. 

The  method  for  converting  from  a  complementary  sequence  pair 

into  the  standard  group  form  is  therefore 


1    c    3  n- 1    n 


1— a, i,&r, . . . o  a      i 


11=  a   a      , a7  (3.1) 

n    n-2  L  w«  x/ 


III=b1b3b5.  ..  ..bn_j 

B=b1b?b,. . . . . b       .b / 

1    <<-    3  n- 1    n        \  T,,   ,     ,  , 

N IV=b   b      -,.....  b0 
n   n-Z  L 

If  the  matrix  (I    II    III    IV)  is  taken  as  the   symbol  for  the  above  grouping 
the  Identity  operation,    I,    should  yield  (I    II    III    IV).      Post  multiplying 
(I    II    III    IV) 

10      0      0 
by         o      0      ^       0      §ives  (T    H    m    IV>'    therefore  I    can  be  taken 
0      0      0       1 

as  the  identity  matrix  of  rank    4,    using  the  operation  matrix  multipli- 
cation on  the  right.     Similarly    Tx  has  the  effect  of  exchanging  the  first 
and  second  columns  yielding  (II  I    III    IV).     Now  post  multiplying 


(I    II    III    IV)    by 

0      0  . 

yields  (II    I     III    IV)  the  desired  result  so  the 


0       10      0. 
10      0      0 
0      0      10 
0      0      0       1 


matrix  just  used  could  represent    Ti  .      Similarly, 
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1000  0100       1000     1000 

T  ..  0   1   0   0  T_  1   0   0   Or.O   1   o   or_o   1   0   0 

2~  0   0   0   1  0001  ^1_  00   1   0   3  0   0   I   0 

0010  0010       0001     0001 

1  0   0   0       1000  0010 

01000100  0001 

1=    0  0   1   o  A2~  0   0   10  10   0   0 

0001       0001  0100 

Matricies  are  often  a  very  convenient  method  of  expressing  trans- 
formations.     However,   when  the  elements  are  few  and  the  numbers  in 
them  simple  it  is  often  more  convenient  to  express  the  operation  dif- 
ferently.    For  instance  (I    II    III    IV)  operated  upon  by    Ti  yields 

(II    I    III    IV).      This  could  be  just  as  well  expressed   as 

I    II    III    IV 
11       II    I    III    IV     ' 

This  indicates  that  the  top  row  is  transformed  into  the  bottom  row  by  the 

operation,   or  in  even  more   shortened  form,    Ti  =  (II    I    III    IV)    with  the 

top  row  understood  which  means,   as  is  common  practice,     Ti  (I    II    III  IV) 

=       (II    I    III    IV).     Another  example  is    Ax  =  (I   II  III  IV),    [which  means 

Ai  (I  II  III    III)  =  (I  II  III"  IV)]   where  the  bar  across     I    and    III    indicates 

the  negative  element  in  the  matrix  or  in  actuallity  the  complement  of  the 

portion  of  the  code  contained  in    I    and     III.      Table     3.  2    is  a  list  of  all 

32    transforms  generated  from  Table     3.  1     both  by  symbol  and  group 

element  representation.     There  are  actually     64    transform  operations 

in  the  entire  group.      However,   the  exchange  operation  (III    IV     I    II)     is 

not  considered  along  with  all  the     32    elements  it  would  generate.      In 

this  paper  the  pair    (I    II)    and    (III    IV)    are  generally  considered  as 

being  an  unordered  pair,   which  eliminates  the  exchange  operation. 

The  second  feature  listed  under  desirable  characteristics,   that 
of  simple  manipulation  of  the  elements  of  the  group  is   satisfied., 

The  supplementary  property  which  is  discussed  in  Chapter    6    is 

also  satisfied  by  the  formulation  of  groups  from  complementary 
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sequences  by  this  method.      The  parity  check,   a.§b.#a    .,      ®b    ,,     .  =  1, 

11      n+l-i      n+l-i 

is  also  satisfied  under  all  transforms  and  is  easily  checked  since  each 

column  of  equation  3.  1     is  one   solution  of  the  parity  check. 

i=  a1a~a[-  .....a... ...a      , 

lib  i  n- 1 

n   n-2  n+l-i     2  (3.1) 

111=  b,b_bc  o...0b.   „  .  „  .  b      , 
13o  i  n- 1 

IV  =  b    b  b b 

n   n-2  '"00"Dn+l-lD2 

I  =  (I  ii   in  iv)  R  =  (I  IT  iv  ITT) 

Ai=  (T  ii  Hi  iv)  q  =  (n  I  iv  in) 

a2=  (i  IT  in   iv)  p  =  (IT  I   IV  HI) 

Tx=  (ii   i  in   iv)  o  =  (TT  T  in  rv) 

ts=  (i  ii  iv   in)  n  =  (ii   i  In  Tv) 

t  =  (II    I   iv    in)  m  =  (IS    IV  III) 

Ci=  (T  fi   in   IV)  L  =  (IT  I  HI  W) 

C2=  (I   II  HI  W)  K  =  (I  Tl  Tv  ST) 

c=finH!Iv)  j=(rTTTvITi) 

z  =  (ii  TTTi  iv)  H=  (i  ii  w  my 

Y  =  (IT  I  T!i   iv)  G  =  (IT  T  iv  in) 

x  =  (IT  i  in  Tv)  f  =  (i  TT  in  iv) 

w  =  (ii  T  in  Tv)  d  =  (T  ii  in  iv) 

v  =  (T   ii   iv   TT!)  O  =  (ii    i  Tv  TTT) 

u  =  (I  ii  fv   in)  B  =  (IT  I   iv  Tfl) 

s  =  (i  IT  iv"  in)  it  =  (nTTv  in) 

E  =  (III    IV    I    II) 
Elements  of  the  unordered  operations  group. 

TABLE  3.2 
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Another  feature  inherent  in  this  formulation  of  the  group,    which  is 
very  useful  in  searching  for  new  kernels,    is  the  interlace  breakdown 
into  the  group.      This  is  useful  because  a  kernel  may  be  considered  as  the 
interlace  of  two  half  length  codes  where,   as  will  be  shown  in  Chapter     7, 
the  number  of  zeros  and  ones  in  each  half  length  code  are  easily  pre- 
dete  rmined. 

Although  this  group  of  operations  is  indeed  a  group,    it  has  never 

been  proved  formally.      In  order  for  a  collection  of  elements  to  be  a 

15 

group,   they  must  satisfy  the  following  four  conditions: 

1.  A  group  has  closure. 

2.  A  group  is  associative  under  the  group  operation. 

3.  A  group  has  an  identity  element. 

4.  Each  element  in  a  group  has  an  inverse. 

That  this  operation  group  has  closure  is  demonstrated  by    Table     3.  3 
which  is  the  full  multiplication  table  of  this  group  of  order    32.     The 
elements  under  the  operation  multiplication  are  associative   since  they 
are  expressible  as  matricies  and  matrix  multiplication  is  associative. 
The  element    I    is  the  identity,   and  as  seen  from  Table     3.  3    all  elements 
have  an  inverse.      Therefore  these     32    elements  have   satisfied  all  the 
requirements  for  a  group. 

There  are  three   special  subgroups  of  order    4    which  within  the 
subgroup  deal  with  only  one  type  of  operation,    Ti  T2  TI,    Ai  As  CI, 
Cx  C2  CI.     All  of  these  subgroups  have  the  unusual  property  that  each 
element  in  the  subgroup  is  its  own  inverse.     This  group  of  order    32 
has  been  identified  as     isomorphic  with  Senior's  group  number    44. 
This  will  be  discussed  in  some  detail  at  the  end  of  this  chapter. 

A  study  of  complementary  codes  in  (I    II    III    IV)  form,    or  as  we 

will  call  it  sequence  quadruple  form,    reveals  many  interesting  proper- 
ties.     Some  of  these  properties  and  their  formal  proofs  will  follow. 
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In  these  proofs  the  complementary  pair  (A3B)  will  be  assumed  to  be  of 

length    n.     The   symbol     L.     will  mean      S      (a /©a      .     .  •  1)    for    all 

_.   1=1  J 

i,    1  <  j  <  n- 1     and  similarly    UD  =     J7    ,    «, 

J         -  J  —  '        B         S     b.'Ob    . .     .  „ 

.    ,      l      n+i-  ] 

i=l  J 

Theorem    3.  1 

If  (A,B)  are  a  complementary  sequence  pair  and  if  either  (I,  II) 
or  (III,  IV)  in  sequence  quadruple  form  are  complementary,  then  the 
other  pair  must  also  be  complementary, 

1.  Assume     (I,  II)     are  comlementary, 

2.  L,    =  UTT    by  the  definition  of  complementarity. 

3.  UT  =  L.XJ    by  the  definition  of  complementarity,, 

4.  Adding  the  equations  in  steps     2    and    3     gives 
LI+Ln  =  UI+UIr 

5.  Since  only  even  spacing  s  are  concerned  in  the  sequence 
quadruple  form,   when  related  back  to  the  original     (A,B) 
pair,    the  likes  of    A    are  equal  to  the  unlikes  of    A    for  all 
even  spacings.      This  will  be  expressed  as     L.    =  U  .     for 
even  space. 

6.  L>A=UR    and    U=L,R    by  the  definition  of  complementarity. 

7.  Therefore     L    =U0     for  even  spacings. 

r>       a 

8.  L    T+UTTT=k.     since  the  total  number  of  matches  at  any  spacing 
is  equal  to  the  likes  plus  the  unlikes,   and  is  dependent  upon    j, 

9.  LIV+UJV  =kj. 

1  0.      Adding     8    and    9    gives     LjTT+UT,j+L)jV+Ujy  =  (L,j^+L,jy)  + 

(unl+uIV)  ,  2kj  . 

1  l'      LIII+LIV  =  kj  '  by  step    7o 

1  2.       LTV-UTTT=  0    by  subtracting  the  equation  in  step     8    from  that 
in  step     1  lc 

1  3.       Therefore     L^  =Uttt    and  the  (III,    IV)  pair  is  complementary. 
Since  the  proof  would  have  been  the   same  if    (IIIS   IV)     were 
chosen  complementary  rather  that    (I,    II),   the  general  state- 
ment is  true0 
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Theorem    3.  2 

If    (A,B)    are  a  complementary  sequence  pair  in  sequence  quadruple 
form  and  if  one  of  the  pair     (I,  II)    or    (III,  IV)    is  not  complementary 
neither  is  the  other. 

1.  Assume     (I,  II)  not  complementary. 

2.  There  are  two  possibilities  for     (III,  IV)    either  complementary 
or  not  complementary. 

3.  Assume  (III,  IV)    are  complementary,    then  by  the  Theorem  3.1 
the     (I,  II)    pair  must  also  be  complementary,   but  this  is  con- 
trary to  assumption. 

4.  Therefore     (III,  IV)     cannot  be  a  complementary  pair. 

>}c  >;«  >;« 

Theorem     3.  3 

A  necessary  and  sufficient  condition  for     (I,  II)    and    (III,  IV)    to  be 
two  complementary  pairs  when  a  complementary  sequence  pair  is 
written  in  the  sequence  quadruple  form  is  if  the  likes  for  even  spacing  of 
one  of  the  code  pair  is  equal  to  the  unlike  s  at  the  same  spacing  for  the 
same  code  of  the  pair.  First  the  sufficient  portion: 

1.  L, .    =  UA     for  even  spacings. 

A  A  r  & 

2.  L,,  +  L.      =  L.     for  even  spacings. 

3.  UT  +  Ujt  =  U.    =   La     for  even  spacings. 

4.  LT  +  LTT  +  Uy  +  UJT  =  ^L.     for  even  spacings  by  adding  the 
equation  in  step    2    and  step     3. 

5.  JLT  +  UT  =  k.     since  the  total  of  likes  and  unlikes  must  equal  the 
possibilities,, 

6.  Ln  +  Un  ,  k.  . 

7.  LT  +  UT  +  L.TT  +  UTT  =  2k.  =  2JL       by  adding  the  equations  in  steps 

1  1  11  11  J  a 

5    and     6     and  comparing  to  those  in  step    4„ 

ith    k.     substituted  for     L, .  . 


LT  +  LTT  =  k.     rewriting  step     2    with    k. 
I  II         J  J 


A' 


9.       LTT  -  UT  =  0     subtracting  the  equation  in  step     5    from  that  in 
step     8. 
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10.  Therefore     !_,„  =  U.    and    ( Ia  II)     are  a  complementary  pair  and  by- 
Theorem    3.1     (III,IV)    are  also  a  complementary  pair. 

11.  For  the  proof  of  the  necessary  condition  assume  that     (I,  II)     are 
complementary  and  that    U      <  L       for  some  even  spacing. 

12.  Therefore     UL  +  U,,  <  LT  +  LTT    for  some  even  spacing. 

13.  However,   by  assumption  (I,  II)    are  complementary  and 
L,  =  UTT    and    U,  =  LJT. 

14.  Adding  the  two  equations  in  step    13    gives     LT+LJT  =  UT+UTT  . 

15.  This  is  incompatible  with  step     12,   therefore  it  is  necessary  that 
U«    =   L.     or    UR  =  Lt-.     for  all  even     spacing s  in  order  to  have 

(I,  II)    and    (III,  IV)    complementary. 

#  #  * 

Theorem    3.4 

If  the  sequence  quadruple  form  of  a  complementary  pair  can  be 
transformed  into  the  pattern  (I    II    I    II),    either  the  original  pair  is  the 
quad  or     (I,  II)    and    (II,  IV)     are  complementary  pairs. 

1.  Given  that  (I    II    III    IV)    is  transformable  into     (I    II    I    II),     and 
is   so  transformed. 

2.  Lt+"^T+LTIT+LTV=^I+^II+^HI+^IV=  ^i    by  the  supplementary 
property,   Theorem    6.  1. 

3.  Ln=l3n     and    Uii=UTT      since    ai®ai+j"  ^i®^i+i    by  Theorem    2-1- 

4.  Therefore     2L  +2L    =  21^+21^=  2k.    by  substituting  into  step     2, 
and  using  the  condition  of  step     1. 

5.  LT+LTT=  k.  . 

6         LT+UT=  k.     since  the  total  of  likes  and  unlikes  must  equal  all 
I       I        J 
possibilities. 

7.  LTT-UT=  0     subtracting  the  equation  in  step    6    from  that  in  step    5. 

8.  Therefore     L    =  Uj    and     (I,  II)    are  complementary. 

9.  Since     (I,  II)     are  complementary,    so  are  (I,  II)    and  any  general 
transformations  of    (I,  U)     and    (111,1V). 

%  *  * 
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This  last  property  should  be  extremely  useful  in  decomposing  long 
codes,    since  it  avoids  the  tedious  task  of  checking  to  see  if     (I  s  II)    are 
a  complementary  pair.     Whenever  a  complementary  pair  expressed  in 
sequence  quadruple  form  has  the  pattern    (I  II    I  TT) ,     certain  additional 
transformations  are  allowed  which  are  not  true  in  general.      These 
allowable  operations  are  the  time  reversal  of    I    and     III,    II    and    IV, 
or  both.     Also  allowed  as  an  operation  is  the  altering  of  either  even  bits 
or  odd  bits  for  all  four  sequences.      The  proof  of  all  these  operations 
being  allowable  is  quite   simple,   but  does  require  an  addition  to  the  no- 
tation to  cover  these  particular  operations.     Double  or  quadruple   sub- 
scripts will  be  used  to  denote  these  operations,     T^-jT^,     and  T, -,-.,. 
represent  respectively  the  time  reversal  of    (ISIII),    (11,  IV)     and  the 
time  reversal  of  both  pairs.     A,  .     represents  the  altering  of  the  odd 
bits  in  all  four  sequences  while    A-,,     stands  for  the  altering  of  even  bits 
in  all  four  sequences.      Actually  complementing  any  two,    or  all  four  of 
the  four  sequences  will  also  leave  a  complementary  pair  when  reassem- 
bled,  but  each  of  these  operations  is  already  accomplished  by  the 
general  transformations. 

Theorem    3.  5 

Whenever  a  complementary  sequence  pair  is  put  into  sequence 
guadruple  and  then  transformed  so  that  the  resultant  group  pattern  is 
(I    II    I    II),    I    and    III,    or     II    and    IV,    or  both,    can  be  time  reversed 
and  the  sequence  pairs  formed  on  reassembly  from  thequadruple,  afterthe 
operations     (T,,     T 2       or    T1234)»     will  still  be  complementary  sequence 
pairs. 

1.  Given,    (A,B)     is  a  complementary  pair  of  length    n    transform- 
able in  quadruple  form  to     (I    II    I  TT),    and  is   so  transformed. 

2.  (I,  II)    is  a  complementary  sequence  pair  by  Theorem    3.4. 
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3.  Therefore     ( I_,  II) ,    (I,  II)    and    (1,11)    are  complementary  pairs 

by  Theorem    2„  3. 

4.  From  this  it  is   seen  that  the  original  sequence  quadruple  after 
the  operations     T13»T24    and    Ti?34    are     (I    n  J    *^' 

(I    U    I    U)     and    (I    II    I  TT).     When  interlaced  in  the  customary 
manner  to  form  pairs,    the  pairs  they  form  will  be  complementary 
sequence  pairs  by  Theorem    2.  80 

*  #  # 

Theorem    3„  6 

Whenever  a  complementary  sequence  pair  is  put  into  sequence  quad- 
ruple and  then  transformed  so  that  the   resultant  group  pattern  is     (I    II 
I    II),    the  operation    A,.,     the  altering  of  odd  bits,    or  the  operation 
A-,  .,    the  altering  of  even  bits,   will  allow  the  quadruple  upon  reassembly 
to  still  be  a  complementary  pair. 

1.  By  hypothesis    A    and    B    are  a  complementary  pair  which 
when  expressed  in  sequence  quadruple  form  are  transformable 
to     (I    II    I  TT),    and  is  so  transformed. 

2.  (I,  II)     is  a  complementary  sequence  pair  by  Theorem     3.  4. 

3.  (I,  II)    when  operated  upon  by  either    A,     or    A^     is   still  a 
complementary  pair  by  Theorem    2.4„ 

4.  If  the  symbol  (  1,  2)  is  used  for  the  pair  (I,  II)  altered  it  is 
seen  that  the  pair  formed  from  the  interlace  of  (12  12)  is 
a  complementary  pair  by  Theorem    2.  8. 

*  *  * 

These  last  two  theorems  should  prove  quite  useful  when  trying  to 
estimate  an  upper  limit  on  the  possible  number  of  composite  comple- 
mentary sequences  of  any  particular  length.      An  example  to  illustrate 
the  use  of  Theorem    3.5    might  prove  helpful.      Taking  pair  number  one 

from  Appendix     II, 
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A=l 101000111011110 

B=1000010010001011 

1=10001011=1 
11=01111011  =  11 
111=10001011  =  1 
IV=1 0000100 =ir 

Since  this  is  already  in    (I    II    I    II)    form  without  transforming  it, 

T,  ,    can  be  applied  directly  giving 

1=11010001 
11=01111011 
111=11010001 
IV  =  10000100. 

Reassembling    gives 

A^llllOOllOlOlOllO 
B^IOIOOIIOOOOOOOH 

which  is  identical  with  number  three  in  Appendix    II. 

As  an  example,   to  illustrate  the  use  of  Theorem    3.6,   applying    A?. 

to  the   same     (A,B)  pair  gives 

1. =11011110 
2. =00101110 
3. =11011110 
7. =11010001, 

reassembling  give 

A2=10110111101U000 
B2=1110001011 101101. 

This  is  the  same  as  number     6    in  Appendix    II    with  the  transformation 
TC2=£>   applied. 
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Unorde  red 

Unordered 

T  ransform 

Equivalent 

Transform 

Equivalent 

1. 

I  =  (I  II 

i  ny 

A2 

17. 

R  =  (I  Ti  Ti  1) 

0 

2. 

A!  =  (l     II 

T  ny 

C 

18. 

q  =  (ii  1   n  1) 

J 

3. 

A2=(I  TT 

i   n) 

I 

19. 

p  =  (IT  i   ii   i) 

T 

4. 

Ti=(II    I 

i  n) 

s 

20. 

o  =  (TT1   l  TI) 

R 

5. 

T3=(I    II 

"n  i) 

X 

21. 

N  =  (II    I    T   II) 

U 

6. 

T  =(II    I 

"TT  I) 

p 

22. 

M=  (T  n  Ti   I) 

Y 

7. 

Ci=(T  IT 

I  IT) 

F 

2  3. 

l  =  (Ti  1  I  ii) 

V 

8. 

C2=(I    II 

1    II) 

D 

24. 

K  =  (I  Ti   ii  1) 

Z 

9. 

c  =(iir 

1    II) 

Ax 

25. 

J  =  (ii  T  ii  1) 

Q 

10. 

z  =(II  T  TIT)) 

K 

26. 

H  =  (I    II    II  ~D 

W 

11. 

y  =(ir  i 

1  IT) 

M 

2  7. 

G  =  7JI  T  TT  I) 

B 

12. 

x  =(ir  i 

i   n) 

T2 

2  8. 

f  =  (i  Ti  1  !!) 

Ci 

13. 

W  =(II  I 

i  ID 

H 

29. 

D  =  (T    II    I    II) 

c2 

14. 

v  =(T  ii 

Tf  I) 

L 

30. 

TT    =    (II      I      II     T) 

0 

15. 

U  =(1    II 

ii  i) 

N 

31. 

B  =  (TT  I  "H   T) 

G 

16. 

s  =(i  TT 

ii  i) 

Ti 

32. 

0  -  (II  T  II    I) 

E  =  (I    Ti    I    II) 

■a 
A3, 1 

Redundance  due  to    (I    II    I    II). 


TABLE     3.4 
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The     (I    II    I    II)    form,    or  one  of  its  transforms,    appears  quite 
often  in  the  study  of  composite  complementary  sequences,,     Since  two 
of  the  components  are  the  same,    and  the  other  two  are  complements, 
it  is  reasonable  to  assume  that  some  of  the  transforms  are  redundant, 
that  is  yield  indistinct  results.      Table     3.4     lists  all  the  transforms 
with  the  elements     I,    II    and  their  complements,     If    (A,  B)    are  taken 
as  an  ordered  pair,    the  number  of  different  codes  formed  from  the 
transformation  group  are     32    instead  of    64    because     As(l    II    I    II) 
and    E(I    II    I    TTj    are  the  same.      If     (A,B)    are  taken  as  an  unordered 
pair  the  number  of  transforms  forming  different  codes  is  decreased 
from     32     to     16.      Table     3.4    lists  each  transform  with  its  unordered 
mate. 

The  only  possibility  of  having  less  than     16     unordered  pairs  of 
codes  which  are  distinct  would  be  if  the  group  form  could  be  transformed 
into  the  pattern    (I    I     I    I).      The  only  case  known  to  the  writer  where 
this  is  true  is  the  quad.      For  one  version  of  the  quad  this  is 

A=ll  B=10 

1=1 
11=1 
111=1 
IV=0. 

Table     3.  5     shows  the  eight  ordered  pairs,   or  if  the  cross  index  column 

is  used  the  four  unordered  pairs,  for  the  quad.      The  four  possible  code 

transformations  for  the  quad  are: 

12  3  4 

A=      11  11  00  00 

B=      10  01  10  01. 
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Pattern 

1.  (I    I 

2.  (I   I 

3.  (I  T 

4.  (I    I 

5.  (Tl 

6.  (T  T 

7.  (I  T 

8.  (I  I 


Transforms  with  pattern  Cross  index  number 

for  ordered  pairs  for  unordered  pairs 


I>  Ti ,  H,  tt  3. 

AlfY,V,B  6, 

A3,W,S,0,       E  1. 

T2,T,C2,N  8. 

Cx,0,K,J  7. 

C,M,L,G  2. 

Z,R,Q,F  f. 

X,U,PSD  4. 

Redundance  due  to    (I    I    I  "I) . 
TABLE     3,5 

This  exhausts  the  transformations  of  special  types  and  the   redun- 
dancies    caused  by    (I    II    I    II).      The  operations  group  will  now  be 
compared  to  Senior's  group  number    44    of  length    32    as  was  mentioned 
earlier  in  this  chapter. 

The  identification  of    any  group  with  a  group  in  Senior's  list,    is 

accomplished  through  establishing  a  one  to  one  correspondence  between 

16  As  will  be  demonstrated,    Senior's  group  num- 

group  generators.  ° 

ber    44    of  order     32     is  in  one  to  one  correspondence  with  the  unordered 
operations  group  of  complementary  sequences  expressed  in  sequence 
quadruple  form.      Therefore  the  two  groups  are  isomorphic. 

Senior's  generators  will  be  given  in  lower  case  letters,   while  those 
of  the  operations  group  will  be  in  their  standard  upper  case   symbols. 
Five  generators  are  required  for  Senior's  number    44;     let  these  be 
a,b,c,d,e.      The  following  relationships  between  generators  are  re- 
quired for  Senior's     44: 

1.  a2=b2=c2=d2=e2  =  l  where     1     is  the  identity. 

2.  b~ab=a,    c-1ac=a,   d_1ad=a        where    x   x  is  the  inverse  of    x. 
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3.  c_1bc=b,   d-1bd=b, 

4.  d    1cd=c, 

5.  e    ^esa,    e    1be=b, 

6.  e    1ce=ca,   e    1de=bd. 

A  check  of  the  above  requirements   shows  that    a,b,c,d    form  an 
abelian  group  of  order     16    while  the     "e"     generator  is  commutative 
only  with    a,b    and  is  no n- commutative  with    c,d. 

Let    a=Ci,    b=C2,    c  =  Ti,   d=T2,   e=Ax  .     From  Table     3.3     it  is  seen 
that  the  following  relations  are  true. 

i.       c?=  c|=  t!=t|=  a?=i 

2.  C2C1  C3  =  Ci,   Ti  Ci  Ti  =  Ci  ,   T3  Ci  T2  =  Ci  . 

3.  Ti  C2  Tx  =  C2  . 

4.  T2  T1   Ts=  Ti  . 

5.  A1C1  Ax  =  Ci  ,       Ax  C2  Ax  =  C2  . 

6.  AiTiAi=     Ti   d  =  0,       A1T2A1=    C3  T2  =    H0 

This  one  to  one  correspondence  between  the  generators  of  Senior's 
group  number    44    and  a  set  of  generators  from  the  unordered  operations 
group  on  complementary  sequences  is  necessary  and  sufficient  to  prove 
the  two  groups  are     isomorphic. 

The  diagrams  for  Senior's  group  number    44    of  order     32    and 
those  for  its   subgroups  were  used  as  a  guide  line  to  check  the  generation 
of  subgroups  from  the  operations  group  of  unordered  complementary 
sequence  pairs.      These   subgroups  of  order     1,2,4,8    and     16    are   shown 
in  Appendix     IV,     All  groups  that  were   shown  in  Senior'  s  diagrams  are 
included  in  these  charts;     however,    Senior's  diagrams  included  only 
normal  subgroups  and,   although  much  care  was  taken  in  the  generations 
of  the  subgroups  of  the  operations  group,   there  is  a  possibility  that 
Appendix    IV    is  not  exhaustive  due  to  some  undetected  non- normal  sub- 
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Senior  in  conjunction  with  others  has  rewritten  his  work  and 

17 
changed  his  designating  method.  Under  this  new  method  his 

number    44    of  length    32     is  now     r,a,.     The  ordered  operations  group 

of  length     64    is  also  identified  as      fl ,. a ,  „ 

An  identification  of  the  group  with  Miller's  number    23     of  order 

32    was  also  accomplished  for  those  who  prefer  the  substitution  group 

1  8 
method  of  identification. 
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CHAPTER    IV 
HAMMING    DISTANCES    OF    COMPLEMENTARY    SEQUENCES 

The  Hamming  distance  is  a  well  known  property  of  certain  classes 
of  codes.  This  very  useful  property  has  had  most  of  its  applica- 

tions in  the  fields  of  binary  error  detecting  and  error  correcting  codes. 
This  chapter  will  demonstrate  the  usefulness  of  Hamming  distances  in 
the  field  of  complementary  sequences,    as  both  the  code  pairs  and  their 
sequence  quadruple  formulation  have  certain  invariances  in  Hamming 
distance  under  various  conditions. 

The  Hamming  distance,   D(U,  V),   of  two  binary  vectors  or  sequences 
U    and    V    is  defined  as  the  number  of  positions  in  which  these'two  binary 

vectors  differ. 

i=n 
D(U,V)  =        S        u .  •    v.  (4.  1) 

i=l 

For  example,   the  Hamming  distance  of  the  two  binary  vectors 

U=1011001 

V=0111001 

7 

is    D(U,V)  =    S     u.     §    v.     =  (1  +  1+0+0+0  +  0+0)  =  2 
v  .    ,      i  i 

i=l 

The  Hamming  weight  or  just  weight  of  a  binary  vector  is  its  distance 
from  the  null  vector. 

Theorem    40  1 

All  complementary  sequence  pairs  of  length  n  have  a  Hamming 
distance  of    n/2. 

By  hypothesis 


and 


./\.— a-i  a -» a  0  0  .  ©  .  & .  a 

1^3  n 


B=blb2b3.„....bn 


A    and    B    form  a  complementary  sequence  pair, 
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1.  All  complementary  sequence  pairs   satisfy  the  parity  check 

ai«an-i+l     •Wi+l     ■     L  for  all 
i,    1  <  i  <  n/2  by  equation    2.  4„ 

2.  Since     mod  2    addition  is  commutative 

a.  «b.  ®a      ...    ©b      • .,   a  1  . 
i        l        n-i+1        n-i+1 

3.  Therefore  either    (a.  <&b.)  =   1     or  (a         .,   lb         ,i)»   but  not 

11  n-i+1        n-i+1' 

both. 

4.  For  each    i,      1  <  i  <  n/2     one  of  the  two  vector  positions  in 
step     3    will  have  a  one  and  the  other  vector  position  a  zero. 

Summing  up  over  all  possible     i,   gives    n/2    for  the 
Hamming  distance. 


When  a  complementary  sequence  pair  is  decomposed  into  a  sequence 
quadruple,    each  pair  of  sequences  in  the  quadruple,   which  form  one  of 
the  complementary  pairs,   has  a  Hamming  distance  which  will  be  shown 
to  be  characteristic  of  the  code  length,,      The   sum  of  the  two  pairs  of 
Hamming  distances  from  the  quadruple  is  an  invariant  and  is  equal  to 
n/2     the  Hamming  distance  of  the  complementary  sequence  pair. 

For  example, 

A=1001010001 
B=1000000110 
are  a  complementary  sequence  pair  of  length     10o 

Decomposing  them  into  the   standard  sequence  quadruple,   (I    II    111    IV), 

give  s 

I  =  10000 

II  =  10110 

III  =  10001 

IV  =  01000 
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I     $    II    =    00110  D(I,II)  =  2 

III     §    IV    =     11001  D(III,  IV)  =  3 

I     ©    II     §    111     9    IV  =11111     D(H>II,  III0IV)  =  5  . 

This  distance  is  half  the  length  of  the  complementary  sequence  pair, 

and  is   required  because  of  the  necessary  condition  of  the  parity  check. 

If  the  distances     (2,  3)    for  the  pairs  given  in  the  above  example  are 
taken  as  an  unordered  rather  than  an  ordered  pair ,   this    (2,3)  pair 
will  be  invariant  through  all  general  transformation  operations. 

All  of  the  general  transformations  on  the  sequence  quadruples 
have  one  or  more  of  the  following  properties: 

1.  Change  the  order  of    (I,  II)    or    (IIIjIV)     or  both, 
for  example,  T1  ,    T2>    T  . 

2.  Complement  the  pairs     (I,  II)     or    III,  IV)  or  both, 
for  example,    Ci   ,C2  ,C  . 

3.  Complement  one  from  the  pair    (1,11)     the  other  from 
the  pair    (III,  IV),   for  example,     Ax  ,    As  „ 

4.  Exchange     (I,  II)     for    (III,  IV),   for  example,     E0 

Since  all  operations  in  the  group  can  be  generated  by  multiplication 
from  the   set  listed  above,   an  examination  of  the  effect  on  Hamming 
distance  by  these  transformations  will  be  adequate  to  prove  the  invar- 
iance. 

Theorem    4.2 

A  change  in  order  of    (I,  II)    or    (III,    IV)     or  both  will  not  change 

the  Hamming  distances  of  the  pairs     (I,  II).    or    (111,1V),, 

n 

1.  By  definition    D(U,V)=      E     i^  •    v. 

i=  1 

2.  u.  §v.  =  v.  ®u.     modulo    2     addition  is  commutative 

1111 

i=n  i=n 

3.  D(U,V)  =      S     u    ®v.    =        T,    v    «u       =D(V,U) 

i=l       x       '  i=l 

*  *  * 
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Theorem    4.  3 

Complementing     (I,  II)     or    (III,  IV)     or  both  will  not  change  the 
Hamming  distances  of  the  pairs     (I,  II)     or    (III,  IV). 

1.  u.  ©  v.     =    u.       •  v.  by    Theorem    2.1 

1111  } 

2.  Therefore     £  u.  9    v.     =    £  u.     ©    v. 

11  11 

3.  and    D(U,V)  =  D(U,V) 

#  *  # 

If  it  is  assumed  that  the  length  of  the   sequence  quadruples  formed 
from  the  complementary  sequence  pairs  have  a  length    K    and  that  the 
distance  of  the  first  pair  is     R,   then  the  following  theorem  can  be  proved: 

Theorem    4.4 

Complementing  one  of  the  pair    (I,  II)     and  one  of  the  pair    (III,  IV) 
results  in  the  exchange  of  their  Hamming  distances,    or  expressed  dif- 
ferently,   gives  the  same  distance  pair  in  opposite  order. 

1.  Let  D(I,II)  =  R,      therefore 

D(III,  IV)  =  K-R 

as  the  parity  check  requires  D(I<HI,  III  ©IV)=K 

2.  Complement  one  of  the  first  pair  in  the     (I  II    III    IV)  form, 
and  also  complement  one  of  the   second  pair. 

Each  vector  position  that  was  formerly  alike  and  summed 
to  zero  is  now  different  and  sums  to  one. 

Similarly  each  vector  position  that  was  formerly  different 
and  summed  to  one  is  now  alike  and  sums  to  zero. 

3.  D(T,  II)  =    D(I,II)     =    K-R    from    2. 

4.  DflllTlV)     =    D(III,TV)     =    K-(K-R)     =R    from    2. 

#  *  * 

Theorem    4.  5 

The  exchange  operation  has  the  effect  of  changing  the  order  of 
Hamming  distances  of  the  pairs     (I,  II)     and     (III,  IV)    when  expressed 
in  sequence  quadruple  form. 
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1.  The  pairs     (I,  II)     and    (III,  IV)     are  not  changed  in  anyway 
except  by  position  in  the  exchange  operation. 

2.  The  only  way  to  effect  the  Hamming  distances  of  a  pair  of 
vectors  is  to  change  some  of  the  components  of  the  vectors. 

3.  Therefore  the  two  weights  are  the   same  but  reversed  in 
order. 


An  examination  of  the  Hamming  distances  for  the  known  kernels  of 
2,    10    and    26    is  now  in  order  as  these  distances  are  a  characteristic 
of  the  kernel,    as  are  the  vectors  formed  in  the  process  of  obtaining  the 
distances  by  the  modulo     2     sums. 

The  first  case  to  examine  is  the  trivial  case  for  the  kernel  of 

length    2. 

A=ll  B=10 

Broken  down  into  sequence  quadruple  representation 

1=1  D(I,H)=0 

11=1 

111=1  D(III,  IV)  =  1 

IV  =  0 

There  are  two  kernels  for  length     10    to  be  examined. 

A=1001010001 
B=1000000110 

1=10000  I§  11=00110  D(I,II)=2 

11=  1 0 1 1 0 

111=10001  III  0  IV=11001  D(III,IV)=3 

IV=01000 

A=0101000011 
B=0000100110 

1=00001  10    11=10010  D(I,  II)  =  2 

11=10011 

111=00101  III©     IV=01101  D(III,  IV)=  3 
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and 


Both  kernels  of  length     10    have  the  unordered  Hamming  distances  of 

(2,3)  when  expressed  in  sequence  quadruple  representation.     However, 

their  "Hamming  vectors",    the  unordered  pairs  (00110,11001)    and 

(10010,01101)    are  different  and  are  characteristic  of  the  kernel  of 

their  origin  under  all  possible  transformation  operations.     Hamming 

vectors  will  be  considered  at  length  in  the  next  chapter  and  proof  of  the 

above   statement  will  be  deferred  until  then. 

There  is  only  one  kernel  of  length     26  (Chapter  7)  and  it  is 

A=0100 1101 11 10 101 11100  11 10  10 
B= 10 11 001 00001 1111 1100111010 

1=0010111110111  I    $    11=0000001101010  D(I,II)=    4 

11=0010110011101 

111=1101001110111  III    ©    IV=1111110010101  D(III,IV)=    9 

IV=0010111100010 

This  is  a  complete  listing  of  the  Hamming  distances  for  all  known 
kernels  and  it  is  exhaustive  for    n=2,  10,  18    and    26c      It  is  worthwhile 
to  note  that  the   sum  of  the  unordered  pairs  of  distances  is  always  equal 
to  one  half  the  length  of  the  code  and  that  the  Hamming  vectors  of  the 
pairs  are  complements.      Both  of  these  conditions  are  due  to  the  parity 
check  being  a  necessary  condition,, 

A  natural  extension  of  the  Hamming  distances  of  kernels  is  the 
consideration  of  Hamming  distances  of  composite  codes  when  put  into 
sequence  quadruple  form.     Six  different  code  pairs  of  length  eight  are 
listed  in  Table     4„  1„      These  pairs  will  be  examined  for  a  possible  in- 
variance.     After  this  invariance  is  noted  a  general  proof  will  be  given 
for  the  characteristic  distance  of  composite  codes. 

The  first  four  examples  in  the  table  when  decomposed  into  sequence 

quadruple  representation   have     (I,  II)     and    (III,  IV)     as  code  pairs  which 

are   still  complementary.      These  codes  are  half  the  length  of  the  original 
pair,   and  being  complementary,   have  a  Hamming  distance  one  half  of 
their  length,   or  n/4  as  compared  to  the  original  codes  of  length  eight. 
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(a)  A=11111001 

B=11001010 


H=(00110011) 


1=1110  I  ©    11=0101  D(I,  II)  =  2 

11=1011 

111=1011  III  9-    IV=1010  D(III,IV)=2 

IV=0001 

(b)  A=11110110  H-(001100in 
B=11000101  H-(OOllOOll) 

1=1101  I    «    11=1010  D(I,  II)  =  2 

11=0111 

111=1000  III    <B    IV=0101  D(III,IV)=.2 

IV=1101 

(c)  A=11101101  H-moooiim 
B=11100010  H-(OOOOllll) 

1=1110  I    ®    II  =0011  D(I,  II)  =  2 

11=  1 1 0 1 

111=1101  III    ©    IV  =1100  D(III,IV)=2 

IV=0001 

(d)  A=11011110  H-fOOOOllin 
B=11010001                       H-(OOOOllll) 

1=1011  I    ®    11=1100  D(I,  II)  =  2 

11=0111 

111=1000  III    ®    IV  =  0011  D(III,IV)=2 

IV=1011 


(e)  A=11101011 

B=11100100 


H=(000011  11) 


1=1111  I    ©    II  =  0110  D(I,  II)  =  2 

11=1011 

111=1100  III    0    IV  =1001  D(III,IV)  =  2 

IV=0101 

(f)  A=10111110  H=(00001111) 

B=10110001  l  ' 

1=1111  I    ®    11=1001  D(I,  II)  =  2 

11=0110 

111=1100  III    ®    IV=0110  D(III,IV)  =  2 

IV=1010 

Some  Hamming  distances  and  vectors 
of  codes  of  length    8. 

TABLE    4,  1 
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An  examination  of  the  last  two  examples  in  the  table   shows  that 
although  the     (1,11)    and    (III,  IV)    pairs  are  not  complementary  their 
Hamming  distance  is   still    n/4    when  compared  to  the  original  codes 
of  length  eight.      That  this  feature  is  true  in  general  will  now  be  proven 
by  examining  and  exhausting  all  the  known  ways  of  generating  composite 
codes. 

Theorem    4.  6 

Every  composite  code  of  length    n     (generated  by  one  of  the  standard 
methods)  when  expressed  in  the  sequence  quadruple  representation  will 
have     D(I,  II)  =  D(III,  IV)  =  n/4. 

There  are  presently  known  four  general  ways  of  generating  com- 
posite codes  from  shorter  cqdes,   (there  is  also  the  method  applicable 
only  to  codes  of  length    2r    given  at  the  end  of  the  proof).     These  are: 

1.  Time  Sequence 
Sx  =  AB 

S2  =  AB 

2.  Interlace 

Ti   =  a,  b,  a0b-,.  .  .  a  b 
1  112    2  n   n 

Ts  =  a  i  b\  a^b-,.  .  .  a  b 
112    2  n   n 

3.  Time  Sequence  Exponential 

cl      c2  cm    dl     d2  dm 

Ux  =  A    i  A    \  .  .  A  mB    _B    *\  .  .  .  B   m_ 

d        d        i  d  i     c  _     c        i         c, 

TT  A    m.     m-1        a     1R    mu    m- 1      „     1 

U2  =  A       A        ....A      d        B  „..B 

4.  Interlace  Exponential 

c,     d,     c-,     d,  c         d 

V1  =  A    VA      B    2..-..A    mB    m 

d        c         d        .     c        .  d       c. 

V2  =  A   mB   mA  "^^B   m-1.....A    lB    l 

Throughout  this  proof  it  will  be  considered  that     (A,B)  are  a  comple- 
mentary pair  of  length    r    and    (C,D)    are  a  complementary  pair  of 
length    m. 
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The  parity  check    a.  §    a    ,  .         §  b.     $  b    ,  .     .  =   1     is  a  necessary 

1  n+l-i  1  n+l-i  7 

condition  for  a  sequence  pair  to  be  complementary.      One  of  the  simplest 
methods  for  making  the  parity  check  is  to  fold  both  codes  in  the  middle 
and  double  them  back  on  themselves.  This  procedure  is  shown  in 

Figure    4.1.     Each 

ala2a3 ai  an/2-ian/2 

a    a       i.  ......  3.    tl     t 

n   n-1  n+1-1    •  •  •  -  *n/2+l*n/2+1 

blb2b3 bi  •••"bn/2-lbn/2 

b   b    -l.......b    ,,     .     ....b    /->.->b    , 

n   n  n+l-i  n/2+Zn/2+l 

Foldover  Parity  Check  Method 
Figure     4."1 

column  in  Figure    4.  1     satisfies  the  parity  check.     If  code  pairs  are 

written  in  the  sequence  quadruple  form  such  as  in  Figure    4.  2   ,   it  is 

seen  that  the  columns  of  this  array  are  the  same  as  the  column  in  the 

first  array  with  rearrangements  in  order  provided  it  is  allowed  that 

the  columns 

a  a 

r  u 

a  and        a 

u  r 

b  b 

r  u 

b  b 

u  r 

are  the  same.     Since  modulo     2    addition  is  the  operation  used  in 
making  the  parity  check  and  in  determining  Hamming  distance,   and  the 
use  of  each  column  involves  only  modulo     2    addition  which  is  commu- 
tative,  the  two  columns  can  be  considered  the  same.      Therefore  the 
fold-over  method  of  grouping  and  the  sequence  quadruple  form  of  grouping 

have  invariant  columns,    order  being  ignored. 
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ala3a5" •"    an-3  an-l 

anan-2 a4       a2 

blb3      -"-bn-3bn-l 

bnbn-2    ••'."••^b4        ^Z 

Sequence  Quadruple  Form 
of  Parity  Check 

FIGURE    4.  2 

Now  consider  the  composite  generating  methods, 

1.  Time  Sequence 

Si  =  AB 

S2  =AB~         writing  in  the  fold- over  parity  check  form 

a,a-)a_.<,...   a      0a      1  a 

12    3  r-Z    r- 1    r  _,  .     -..  n 

b   b b-      b?      b,  D(AsB)  =  r/2 

r   r-  1  3        2        1 

12    3  r-2    r- 1    r  _...     5x           /9 

c    t  c         c         c  D  ASB     =  r/2 

b   b       .  b~      b,      b,               \     »     / 

r    r-  1  3         2         1 

The  top  pair  are  complementary  and  therefore  have  a  distance  of 
r/2,    similarly  the  bottom  two  rows  are  complementary  and  have  a  dis- 
tance of    r/2    by  Theorem    4„  1.     To  write  the  time   sequence  code  in 
sequence  quadruple  form  would  not  change  the  distance  since  it  would  be 
purely  a  reordering  of  the  parity  check  columns.      Therefore     D(I>II) 
would  equal     r/2    as  would    D(I1I,IV).      In  terms  of  the  composite  code 
of  length    n    D(I,  II)=D(III,  IV)=  r/2  =  n/4    whether  or  not     (I,  II)     and 
(III,  IV)     happen  to  form  complementary  pairs„ 

2.  For  the  interlace  scheme  the  (I,  II)  pair  and  the  (III,  IV)  pair 
are  obviously  complementary  sequences  by  construction  and  will  have  a 
distance  of    r/2  =  n/4. 

3.  Writing  the  time  sequence  exponential  pair  in  fold-over  form 


gives  1,2  *     m 


cl     c2  c 

U]_— A       .A.      •  «•,..   -A. 

d         d  j  where  the   symbol    B     means     B    time 

n    m-R    m~  *  R  1  reversed. 
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d         d  d, 

U^A^A.™:.1 A_l 

B    *  B    2 B    m 

The  B's     are  all  time   reversed  but  the  time   reversal  of  one  of  the 

complementary  pairs  does  not  change  the  complementary  property 

D(A,B)=D(A,  B)  =  r  /2.      For  each  bit  of  the    C     code  and    D    code  there 

is  a  distance     r/2     since  every    A    is  matched  with  some  transformation 

of    B    for  each  of  these  bits.      There  are  a  total  of     m    bits  in  the    C 

code  and  in  the     D    code,    therefore  the  total  distance  of    Ui    folded  back 

on  itself  is     m  x  r/2  =  mr/2.     Similar  reasoning  gives     U2    folded  back 

on  itself  a  distance  of    mr/2.      A  rearrangement  of  the  columns  will  not 

change  the  Hamming  distances       Therefore     D(I,Ji)=mr/2    and 

D(II.l,  IV)  =  mr/2 ,    since  the  total  length  of  the  code  is     n=2mr,    D(I,II)  = 

D(III,  IV)=n/4. 

The  interlace  exponential  scheme  is  the  last  general  method  of 

code  generation  to  be  considered.      Written  in  the  fold-over  form  it  gives 
c,     d,     c-,cU  c       /->d       ,-. 

V1=A    XB    lA   2B    2.....  A   m/2B    m/2 

BdmAmBdm-}....    Bdm/2+1ACm/2+1 

V2=AdmBCmA^m-!....  A^/2  +  lBCm/2+l 
Ci       d,       c-,  c      /  ->         d/9 

B    l  A    l   B    2  B    m/2      A   m/2 

Again  each    A    is  matched  with  some  transform  of     B,   D( A.  B)=D(A,  B)= 
D(A,B),   etc.    =   r/2.     As  before,    rearranging  columns  into  sequence 
quadruple  form  will  not  effect  the  Hamming  distances.      The  distance  of 
the  top  two  rows  is     mr/2    and  the  bottom  two  rows     (V2)    also  is  equal 
to    mr/2.      The  total  length  of  the  code  is     2mr,     therefore 
D(I,  II)=P(III,  rV)=mr/2=n/4. 

5.  When  composite  codes  are  of  length    2    ,    they  are   sometimes 

formed  in  a  special  way  from  the  interlaced  sections  of  a  complementary 
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pair.     These  sections  or  pieces  of  the  codes  must  be  of  length    p=2J 

where     j<r.      Using  the  symbols    A   =a,a,a,a/.....a    ,     A->    =a    ,,a    .->... 
to  '  p      1    2    3    4  p  2p      p+1    p+2 

a?       and  so,    the  following  pair  would  be  complementary  if    A    and    B 

formed  a  complementary  pair. 

S,    =A    B   A-,    B.      .....  A       ,    B       , 

lp       p    p    2p    2p  np/p    np/p 

S.    =A   B"   A,    B-,  A      /   B       . 

2p       p    p    2p    2p  np/p    np/p 

Rather  than  try  to  use  the  general  term  to  prove  that    D(I,  II)=D(III,  IV) 
=  n/4    as  was  done  in  the  previous  cases,   a  short  general  example  of 
length    8     will  be  used.     These  two  codes     (A,B)    of  length  eight  will  be 
interlaced  four  bits  at  a  time     (p=2)     to  form  a  code  of  length     16. 

S14=ala2a3a4blb2b3b4a5a6a7a8b5b6b7b8 
S24=a1a2a3a4b1b2b3B4a5a6a7a8b5B6B7b8 

Writing  this  composite  pair  in  fold-over  form  yields 

S14=ala2a3a4blb2b3b4 
b8b7b6b5aga7a6a5 

S24=ala2a3a46lb2b3b4 
b8b7B6b5a8a7a6a5 

D(A,B)=4    by  Theorem    4.  1.     S,.     folded  over  is  just  a  column 
rearrangement  of  the  complementary  pair    (A,B),    and  since  a  re- 
arrangement of  columns  does  not  effect  Hamming  distance,    the  Hamming 
distance  of    S,.     folded-over  is  also    4.      Similar  reasoning  holds  for  S^^ 

Since  the  standard  sequence  quadruple  form  is  again  just  a  re- 
ordering of  the  columns,   the  distance  of  the  (I,  TI)     pair  and  the 
(III,  IV)    pair  is    n/4.     For  this  example, 

I=a, a0b, b,ara7brb7 

13    13    5    7    5    7  D(IfII)=4 


II=b8b6a8a6b4b2a4a2 
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III=a   a  b1b3a5a  b5b? 

D(III,  IV)=4 

IV=bftbAaftaAb  Ab?aAa-, 

86864242  D(I®II,III0IV)=8. 

Therefore  any  composite  complementary  sequence  pair  when  written  in 

sequence  quadruple  form  has    D(I,  II)=D(III,  IV)=n/4. 

*  *  # 

Although  there  is  no  proof  that  the     4     standard  generating  methods, 
plus  the   special  one  applicable  only  to    2       codes,   are  the  only  way  to 
form  complementary  codes  which  are  composite,   an  exhaustive  search 
of  code  lengths     16    and     20     revealed  no  codes  which  were  not  formed 
by  these  standard  generating  methods. 

The  full  usefulness  of  the' Hamming  distance  property  is,   of  course, 
not  yet  known.      The  property  of  the  distance  of  the  complementary  pair 
being  equal  to  half  the  code  length  as  shown  in  Theorem    4.  1  ,     is  used 
in  many  of  the  proofs  in  the  chapter  on  Hamming  vectors.      This  property, 
D(A,B)=n/2,   indicates  that  the  code  pairs  are  orthogonal,,     Because  of 
this  orthogonality  property  the  leakage  of  an    A    code  carrier  into  the 
B    code  filter,    or  vice  versa,    will  have  no  effect  on  the  output  at  exact 
match. 

The  primary  use  the  writer  has  found  for  the  Hamming  distance  of 

the     (I,  II)    and  (III,  IV)     pairs,    has  been  the  error  check  that  this  invariant 

property  offered  in  decomposing  codes  by  hand.     One  possible  use  for  the 

invariance  of    D(I,II)  and    D(III,IV)     could  be  in  the   search  for  kernels  of 

new  lengths.  If  a  characteristic  Hamming  distance  for  unknown  kernels 
could  in  some  way  be  deduced,  based  only  on  length,  this  invariant  prop- 
erty of  distance  would  act  as  a  most  powerful  screen  in  the  search.  The 
(I,  II)  and  (III,  IV)  pairs  are  themselves  orthogonal,  unless  the  code  is  a 
kernel.  This  may  provide  other  possible  methods  for  searching  for  new 
codes.      It  is  difficult  to   say  in  advance  what  utilization  will  be  made  of  any 

invariant  property,   but  it  is  fairly  safe  to  say  that  the  property  will  find  use. 
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CHAPTER    V 
THE  HAMMING  VECTORS  OF  COMPLEMENTARY  SEQUENCES 

The  Hamming  distance  between  the  code  pairs  in  a  complementary 
sequence  was  shown  by  Theorem    4.  1     to  be  one  half  the  length  of  the 
code.      This  Hamming  distance  was  obtained  by  first  finding  what  I  shall 
call  the  Hamming  vector,    H(A,B)=  A  #  B    or  for  convenience  just 
H=A§B,    and  then  counting  the  number  of  ones    H    contains.      The  Ham- 
ming vector  is  itself  useful  in  recognizing  the  composition  of  a  composite 
code  when  it  is  desired  to  break  the  code  down  into  shorter  lengths.      This 
decomposition  might  be  used  to  find  previously  unknown  kernels,    or  to 
set  an  upper  bound  on  the  possible  number  of  pairs  of  codes  of  any  length. 

An  examination  of  each  of  the  methods  for  generating  longer  codes 
from  shorter  ones  as  explained  in  Chapter    2    will  show  that  each  method 
has  a  characteristic  Hamming  vector  form,   although  it  may  be  necessary 
to  time  reverse  one  of  the  codes  to  obtain  this  form.     The  method  of 
generating  codes  of  length    2    ,   which  does  not  apply  to  complementary 
codes  in  general,    will  also  be   studied  with  its  characteristic  Hamming 
vector. 

In  the  following  proofs     (A,B)    form  a  complementary  code  pair  of 
length    n    while     (C,D)    form  a  complementary  code  pair  of  length    m. 
The  symbol     (0)       or  just     0       will  mean  a  string  of    n    zeros,    similarly 
(01)       will  symbolize  a  string  of    n    zero  ones,   and  (a.#b.)       will  be  a 
string  pf    n    elements  equal  to    (a  §b.)  . 

Theorem    5.  1 

The  time  sequence  form  of  generating  a  composite  code  has  the 

characteristic  Hamming  vector    H=(0)    (1)    , 

Use  the  time  sequence  generating  method  on    (A,B)  to  form 

Si=  AB 

Sa=  AB  .  62 


Take  the  Hamming  vector  of    (S1  ,    S2),    H=SL  ®  S2=  0nln,     this  is 

(5.1) 
seen  to  be  true  by  inspection  and  definition  of  modulo     2     sum. 

*  #  * 
Theorem    5.  2 

The  interlace  form  of  generating  a  composite  code  has  the  charac- 
teristic Hamming  vector    H=(01)    . 

Use  the  interlace  generating  method  on    (A,B)    to  form 

Tx=a,  b,  a-,b-.a0bQ a   b 

l     l    c    c    5    5  nn 

T2=a1b1a2b2a3b3....e...0anSn. 
Take  the  Hamming  vector  of    (T^Ts),   H=TX  9  T2=(01)n,  (5.2) 

this  is   seen  to  be  true  by  inspection  and  definition  of  modulo  2   sum. 

*  *  * 

Theorem    5.  3 

The  time  sequence  exponential  form  of  generating  a  composite  code 

n  ,      m  n  .n 


has  the  characteristic  Hamming  vector    H=(c,§d     )     (c-,©d        ,) 

&  1       m  2      m-1 

(c     ®d,)n  (d,   Gc      )n.  ...d     fl>Cl)n  . 
m      1  1        m  m      1 

Use  the  time  sequence  exponential  method  on    (A,B)    and    (C,D) 

to  form  .        , 

c,     c,  c         d,     cU  d 

U1=A    lA   2. ...A   mB_    1B_2..........    B_m 

d         d        ,  d,     c         c        ,  c  i 

u3=a  mA  T:.1 a  xb  mB  m;A .......  B  l . 

Take  the  Hamming  vector  of    (Ui  ,    U2  ), 

H=U!*  U2  =(Cl»dm)n    (c2®dm_1)n.  .  ...(c^Ad^cj". 


(d     0c, )n. 
m        1 


(5.3) 


This  is   seen  to  be  true  by  inspection  and  definition  of  modulo  2  sum. 
*  *  >!« 

It  is  to  be  noted  that  each  of  theTheorems     5.1,5.  2,  5„  3    have 

Hamming  vectors  which  are  antisymmetric  about  their  center.      This  is 

required  to  satisfy  the  fold- over  method  of  making  the  parity  check. 
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Theorem    5.  3    also  shows  that  along  with  the  anti-symmetric  characteristic 

the  Hamming  vector  for  time  sequence  exponential  composition  has  unit 

clusters  of    n    zeros  or     n    ones  and  is  also    anti-symmetric  in  both  halves. 

This  is  easily  seen  because  the  first  sum  is    (c  ,  ®d     )    while  the  last  sum 

1      m 

to  operate  on  the    A    code  is     (c     <&d,  ),     but  it  is  known  through  the  parity 

check  that    c,©d     ©c      Qd.    =   1  .      Therefore  one  of  the  sums     (c ,  <&d     )     or 
1      m     m      1  1      m 

(c      ©d.  )    must  be  zero  and  the  other  must  be  one.      To  summarize,   for 
m      1 

the  time  sequence  exponential  form  of  composition  the  first  half  of  the 

Hamming  vector  is     ant i -symmetric  about  its  center  (the  quarter  length 

point)  and  the  first  half  and  the  last  half  of  the  vector  are  also  anti-sym- 
metric. 

Theorem    5.4    will  concern  the  interlace  exponential  method  of 
forming  composite  codes  and  it  will  be   seen  that  this  method  will  again 
form  a  Hamming  vector  anti-symmetric  about  its  center,   but  with  a  dif- 
ferent arrangement  of  the  clusters  of  zeros  and  ones. 

Theorem    5.  4 

The  interlace  exponential  method  of  generating  a  composite  code  has 

the  characteristic  Hamming  vector    H=(c,®d     )    (d.®c      )    (c-,©d        ,)    .... 

°  1       m  1      m  c      m- 1 

(c     »d,)n(d     »c,)n. 
m       1  m      1 

Use  the  interlace  exponential  method  on     (A,B)     and    (CSD)     to  form 

c,     d,     c->     d->  c         d 

VX-A    *B_A      B    ?....   A   mB    m 

d         c         d         ,  d,     c , 

V2=A        B        A       ......   A      B      . 

Take  the  Hamming  vector  of    (Vi,Vs) 

H=Vx«>  V2=  (c,»d     )n(d,§c      )n(c?®d        ^"....(c     ®d,)n(d     §c,)n      (5.4) 
x  v    1      m  1      m         2      m- 1  ml  ml 

This  is   seen  to  be  true  by  inspection  and  definition  of  modulo  2  sum. 
*  *  * 

The  interlace  exponential  method  again  gives  a  Hamming  vector 

which  is  a nti -symmetric  about  its  center  due  to  the  condition  of  satisfying 
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the  fold-over  form  of  the  parity  check,    but  its  clusters  of  zeros  or  ones 

are  now    Zn    in  length  rather  than  just    n    as  was  the  case  for  the  time 

sequence  exponential  method  of  generation,,      This  characteristic  of 

clusters  of    Zn     in  length  is  caused  by  the  necessary  condition  of  the 

parity  check  for  the     (C,D)  code  pair.     For  example,    the  first  two  clusters 

are     (c,#d     )       and    (d,®c      )       and  the  second  two  clusters  are  (c.-,®d        ,) 
1       m  1      m  Z      m-  1 

and    (d?®c         ,)    .     Applying  the  parity  check  to  the     i=l     and    i=Z    bits 

in  the     (C,D)    pair  gives  the  following  two  equations: 

c,®c     ©d,»d     =1  and         c.§c        ^d.&d        n=l   . 

1      m       1      m  c      m-1       l      m- 1 

Adding  one  modulo    Z    to  both  sides  of  both  equations  and  rearranging 
gives    (c1»dm)  ^(djGc^l)  =0    and    (c2Wm_1)  •(d2»cm_1»l)=0,'    It  is 

true  in  general  that    x  (ft  1  =  x;    therefore     (c ,  ®d     )  4Md,<&c      )=0    and 

°  1ml       m 

(c,<ftd        .  )  ©  (d,(ftc         ,  )  =  0  .      The  sums     (c .   %  d     )     and    (d,    9  c      )    must 
4    Z      m-1         v    Z      m-1  K    1  m1  l    1  m 

be  alike  to  satisfy  this  constraint  with  each  of  the  sums  giving     n    ones 

if  unlike,    or    n    zeros  if  alike,   and  since  the  pairs  are  in  time   sequence 

the  clusters  must  be  in  groups  of    Zn„     A  similar  reasoning  will  show 

that  the  parity  check  for    i  =  Z    gives  identical  results  for  the  third  and 

fourth  terms  in  the  Hamming  vector  for  the  interlaced  exponential  form 

of  generating  composite  codes.      Therefore  equation    5C  4    can  be  written 

H(V1,V2)=  (cj  (ft  dm)2n    (c2»dm_1)2n..„0(dm»c1)2n0  (5.5) 

Before  going  into  the  special  case  applicable  only  to  codes  of  length 

Z       it  might  prove  helpful  to  give  an  example  of  generating  a  code  by  one 

method  and  then  decomposing  it,    to  show  that  it  could  have  been  composed 

by  a  second  method„      The  time  sequences  exponential  scheme  will  be  the 

generating  method  used. 

A=10  C=1001010001 

B=ll  D=1000000110 

Ui=100101  1001 100101011011000000000000111 100 

U3=01 1010010101010101 10001111 11001100111100 

H=U!(BU2       =11111111001100000000111111  11001100000000 
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In  passing  it  should  be  noted  that  the  shortest  clusters  of  zeros  and  ones 
in    H    are  of  length    2    and  that    H    is  anti- symmetric  about  its  center. 

If    Ux   is  now  time  reversed,    indicated  by    U,    ,      the  complementary 
pair  is  now: 

U!=ooi  11 10000000000001 101 101 01  oo  noon  01 001 

U2=01 10100101 01010101 10001111110011001 11 100 
H=  U     9  Us =01 01010101 010101 01 0101010101 010101010101   =  (01)20 
This  form  of  the  Hamming  vector  indicates  that    H    is  composed  of  an 
interlace  pair  of  length    20.      Decomposing  both    U,      and    Us    by  assum- 
ing an  interlace  pair  gives 

Ux=01 100000010111010110 

0110000001 1000101001 
H  =0000000000111111111  I  =  o10i10 
U2= 01 1000000101 110101 10 

1001111  11001110101 10 
H  = 1111111 1110000000000  =  i 10o10 

The  Hamming  vector  of     U,     indicates  that  this  code  of  length    20 

was  composed  by  the  time  sequence  method  from  codes  of  length     10. 

If  the  form    AB    is  assumed  for  the  codes  of  length     10     then 

A=0110000001 
B=0111010110 

This     (A,  B)     pair  is  seen  to  be  the  transformation  ETC,=  EG      on  the 
original  kernel  of  length     10     which  was  used  to  generate  the  code  by  a 
different  method,,      Similarly  the  pair  from    U2  is  also  a  transformation 
of  the  original     (C  ,  D)  pair.      Therefore  the  given  code  of  length    40 
could  have  been  generated  by  two  different  method s0 

A  second  example  will  now  be  given  using  the  same     (A,B)    and 
(C,D)     kernels  as  were  used  in  the  first  example,   but  his  time  utilizing 
the  exponential  interlace  scheme  of  generation  rather  than  the  time  se- 
quence exponential  form. 

DO 


A=10  C=  100 101 0001 

B=ll  D=1000000110 

Vi= 10 110 10001 001 00001 001 0000 10001 1101  111000 
Va=01 001 01 110 1101 1101 0001 1101 0001 11 01 11 1000 

H=VX§V2       =111111111111111100001 . 1 110000000000000000 

Note  that  the  smallest  cluster  is  four,    as  would  be  expected  of  an  inter- 
lace exponential  scheme  with  the     (A,B)     pair  of  length  two.      The     H 
vector  is  again  anti-  symmetric  about  its  center.      Time   reversing     Vx 
gives 

V    =0001  111  01 110001 0000 1001 00001001 0001 01 101 
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V2 =0100101 1101 101 1101000 11 1010001 1101 11 1000 
H=Vlt®V2      =0101010101010101010101010101010101010101  =  (01)' 

As  before  this  form  of    H    indicates  an  interlace  code  pair.     Decomposing 

both    V,     and    V?     into  interlace  code  pairs  gives 

Vj =001 11 10 1000 1000 101 10 
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01101000010001000011 
H  =01010101010101010101  =  (01) 
V2=  00111101000100010110 

101101111011101 11100 
H  =10101010101010101010  =  (10)10 

V,     and    V2  both  break  down  into  interlaced  pairs  of  length     10.     The  top 
row  of    V       breaks  into    A=01  10000001  ,   B  =  01 1 1  010  1  10    which  is  the 
transformation    EG    operating  on  the  original    (CSD)    pair  as  in  the 
previous  example.      Each  of  the  codes  for    V2    would  break  down  in 
similar  fashion  to  some  transformation  of  the  original     (CSD)  code  pair. 

One  use  for  the  method  just  shown  could  be  to   search  for  larger 
kernels  than  those  currently  known.     An  example  might  be  kernels  of 
length  fifty  found  through  the  use  of  time   sequence  exponential  or  inter- 
lace exponential,   where  both    (ASB)    and    (CSD)    are  codes  of  length  10. 
These  would  form  codes  of  length    200,      From  this  as  a  starting  point, 

a  decomposition  to  codes  of  length     100    and  then  to  length     50    might 
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be  possible  through  this  technique.      This  technique  will  be  discussed 
much  more  thoroughly  towards  the  end  of  this  chapter,   but  it  is  neces- 
sary to  consider  the  special  case  of    2       code  composition  and  the 
Hamming  vectors  of  the   sequence  quadruple  form  before  continuing  with 
this  interesting  possibility. 

The  special  generating  method  which  is  applicable  only  to  codes  of 
length    2       is  actually  a  generalization  of  both  the  time  sequence  and  the 
interlace  methods.      As  mentioned  in  the  previous  chapter,   this  method 

takes  a  code  pair  of  length     2       and  from  them  forms  a  code  pair  of 

r+1 
length    2  .      The  allowable  combinations  are  formed  by  interlacing 

pieces  of  the    A    code  and  the     B     code  of  length    2      ,    where  m=,  0,  1  ,  2,  .  ,r. 

As  an  example  a  code  of  .length     8    can  be  used  to  generate  codes  of 

length     16.      In  this  case     r=3    and    m    has  the  possible  values     0,  1,2,3. 

These  values  of    m    give  section  lengths  of     1,2,4,8     respectively. 

B=b1b2b3b4b5b6b7bg 
for    m=0, 

SiraTbla2b2a3b3a4b4a5b5a6b6a7b7a8b8 

S21=albla2B2a3b3a4b4a5b5a6b6a7b7a8b8 
H=S11§S21     =0101010101010101  =(°1)8  =  (0    I1)8 
This  is  exactly  the  same  as  the  interlace  procedure  which  is  a   standard 

generating  method. 

The  section  lengths  are     2    bits  long  for     m=l0 

S12=a1a2b1b2a3a4b3b4a5a6b5b6a7a8b7b8 

S22=ala2blb2a3a4b3B4a5a6b5b6a7a8b7B8 
h=s12®s22     =0011001100110011  =  (0212)4, 
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m=2  gives 

S14=ala2a3a4blb2b3b4a5a6a7a8b5b6b7b8 

S24=ala2a3a4blb2b3b4a5a6a7a8b5b6S7b8 
H=S140S24     =0000111100001111  =  (0414)2, 

For    m=3    the  lengths  are  eight,   which  is  the  entire  code. 

Slg=AB 

S28=AB, 
This  is  the  same  as  the  time  sequence  form  giving  a  Hamming  vector 

H=0000000011111111  =  (0818). 
The  Hamming  vectors  from  this  special  construction  are  altering 
clusters  of  zeros  and  ones,   each  cluster  being  the  length  of  the  segment 
used  in  the  generation. 

Theorem    5.  5 

r+1 

Complementary  sequences  of  length    2  formed  from  an    (A,B) 

pair  of  length    n=2     using  the  special  generating  method  with  segments 

Zn/m 
of  length     m  have  a  characteristic  Hamming  vector    H=(0      1      ) 

Use  the   special  generating  method  utilizing  segments  of  length    m 

of  the  pair    (A,  B)     to  form 

S,      =a1a-)a_0..a       b^b-j^.b       a        .a         ,        a7       b        ,b      ,->...  b 
lm       12    3  m      1    2  m     m+1     m+2...     2m     m  +  1    m+2  n 

i?T     =a1a-)a0..„a       b,  b,, ,  ,  u       a      .  i&      .->.»«  a  -,       d      ■  i  t)      ,->...  b     . 
2ml23  ml2  m     m+1    m+Z  Zm     m+1    m+Z  n 

Take  the  Hamming  vector  of    IS,      ,    S->      ) 

&  lm        2m 

H=Sl      *S7       =    (0mlm)2n/me 
lm       2m 

This  is  seen  to  be  true  by  inspection  and  definition  of  modulo  2    addition. 

>!<  'fi  * 

It  now  seems   suitable  for  completeness  to  include  a  theorem  which 
has  been  alluded  to  many  times  but  not  proven.      This  is  the  anti- sym- 
metric property  of  all  Hamming  vectors  of  complementary  sequences. 
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Theorem    5.  6 

The  Hamming  vector  of  a  complementary  sequence  pair  has  a  pat- 
tern which  is  anti- symmetric  about  the  center  of  the  Hamming  vector. 


12    3    4 
12    3    4 


a       ,a 
n-  1    n 


.......     b       .b  is  a  complementary  pair  of 

n"    !         n  1  4.1. 

length    n 

2.       H(A,  B)=  h.h^h^h..  .  „  h    /„,.h       by  definition  of  Hamming  vector 
1234  n/2n7  to 


3. 


a .9  a    ,  ,     .9  b.   «  b    ,  ,    .  =  1 
l  n+l-i        i  n+l-i 


4.       a.  ©  b.  =  h. 

1         »  1    and 


a n+l-i  §bn+l-i  =  hn+l-i 


by  equation    2.  4 


both  by  definition  of    h. 


5. 

6.       h.  §    h 


n+l-i 


.   =1 


7.       Add  (h       ,    .   ®  1)  to  both  sides  of  the  equation 


8. 


h    ,  .     .  =  h. 
n+ 1  - l  l 


Make  a  change  of  reference  to  the  center  of  the  vector,    which  lies 

between  m/2    and     m/2  +  1,    rather  than  the  end  by  having     i=n/2  -r    then 

n+l-i=n/2+l   +  r 

9.       h    ,-  =  h    ,,  +l+r    for  all    0  <  r  <n/2  -1  by  step    8. 

n/ 2  -rn/2  —      _  /  r 

1  0.      Therefore  the  Hamming  vector  is  anti- symmetric  about  its 
center  when  it  is  formed  from  a  complementary  pair. 

*  *  * 

The  Hamming  vectors  of    (I,  II)     and    (III,  IV)     in  the  sequence  quad- 
ruple form  are  invariant  under  transformation  if    H(I,II)     and    H(III,  IV) 
are  taken  as  an  unordered  pair.     The  following  four  theorems  which 
parallel  Theorems    4.2,4.3,4.4,4.5     on  Hamming  distances  will  dem- 
onstrate this  characteristic  of  unordered  Hamming  vector  pairs. 

Theorem    5.  7 

A  change  in  the  order  of    (I,  II)     or    (III,  IV)     or  both  will  not  change 
H(I,II)    or    H(III,IV).  7Q 


1.  H(U,V)=u.  *  v.  by  definition. 

2.  u.  ®  v.  =  v.  %  u.   since  modulo     2    addition  is  commutative, 

1111 

3.  therefore    H(U,V)  =  H(V,U) 

*  *  * 

Theorem    5.  8 

Complementing  the  pair  (I,  II)  or  (III, IV)  or  both  will  not  change 
H(I,II)    or    H(III,IV). 

1.  H(U,V)  =  u.  «  v.  by  definition, 

2.  u.  $  v.  =  u.  ©  v.  by  Theorem    2.  1, 

l  i  l  i  J 

3.  therefore     H(U,  V)  =  H(U,  V)„ 

#  #  # 

In  the  following  proof  it  is  assumed  that  each  vector  I,  II,  III,  IV  is 
of  length  m.  It  is  also  necessary  to  show  that  x  =  x.  This  is  easily 
seen  sincexl  l=x    and    x  #     l  =  x,     therefore    x  =  x®(l#    l)  =  x. 

Theorem    5.9 

Complementing  one  of  the  pair    (I,  II)     and  one  of  the  pair    (III,  IV) 

exchanges  the  Hamming  vectors    H(I,II)     and    H(III,  IV). 

I—  a1a0ar.......a      ■• 

13    5  n- 1 

11— a        a  -)oooeo.»ea.-4 

n   n-  2  2 

13    5  m-  1 

IV=b     b     -2......b, 

mm  c 


m 


1.  H(I,II)=h.h?h_h,  h 

2.  H(III,  rV)=h1h->h-h/1.  .  .  .  .  .  .  .h        ,  f  . .  ..        .       , 

y       '       '       1    2    3    4  m    because  of  the  parity  check. 


h     . 

m 


3.  H(I,  II)=H(I,II)    =h1h2h3h4.  . 

4.  H(!Ti,  rv)  =H(in,Tv)  =  K  ^R^R^.  . .  Km  =h1h2h3h4. . .  , hm 

5.  Therefore     H(I,  II)=H(lTl,  IV)  =H(III,TV), 

6.  and    H(III,  IV)  =  H(l,  II)  =H(I, TT). 

*  *  * 
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Theorem    5.10 

Exchanging  the     (I,  II)    pair  with  the     (III,  IV)    pair  exchanges  their 
Hamming  vectors. 

1.  A  change  in  the  Hamming  vector  can  only  be  caused  by  a  bit 
change. 

2.  Exchanging  the     (I,  II)    pair  with  the     (III,  IV)    pair  changes 
position  but  not  bits. 

>:=  *  * 

Earlier  in  this  chapter  codes  of  length    40    generated  from  the  quad 
were  decomposed  into  kernels  of  length     10    by  use  of  the  Hamming 
vector.      The  conjecture  was  made  at  that  point  that  perhaps  codes  of 
length    200     manufactured  from  kernels  of  length     10    could  be  decomposed 
into  kernels  of  length    50.      In  order  to  prove  this  decomposition  is  not 
possible  several  more  theorems  are  necessary.     These  theorems  will 
be  proved  before  proceding  to  the  decomposition  problem. 

Theorem    5.11 

Under  all  possible  general  transformations  of  a  complementary  pair 
where     n^  2    there  are  only  four  possible  Hamming  vectors,   two  of 
which  are  complements  of  the  other  two. 

1.  (A,B)     are  a  complementary  pair  of  length    n    where     n  >  2, 

2.  H(A,  B)  =  H       interlaced  with    H    ,    where     H       are  the  odd  bits 

1  '  o  e  o 

of  the  Hamming  vector  and    H       are  the  even  bits  of  the  Hamming  vector 

in  reverse  order.      The  symbol     H=H.    *    H.     will  be  used  for    H     inter- 
laced with    H       in  exactly  the  same  manner  that    (I,  II)    and    (III,  IV)    are 
e 

interlaced  to  form    A,B.     Also     h.  =  a.$b.     as  before. 

3.  Let  H1=(I®III)=a1®b1,a3®b3.  . .  an_]»bn_1  =  h^hg..  .  hR_  l , 

H2=(II§IV)=an®bn,an_2®bn_2o..a2®b2  =  h^...^, 
H3=(I*  rV)=ai0bn, a3«bn_2.  ....  an_ 1*b2=H1K3¥5. .  .  *n_ l , 
H4=(Il*ni)=an«b1 ,  an_  2®b3.  ....  a2»bn_ ,  =^n_  y  .  .  ^. 
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4. 
therefore 

5. 
therefore 

and 

6. 

7. 
8. 

9. 
10. 
11. 
12. 


U.®  V.    =  u.®  v. 
11         11 


by  Theorem    2.  1 


(I©III)=H1  (miV)=H2         (I«IV)  =  H3  (II«III)=H4. 

If    h.  =  u.©v.,    then    h  .®  1=  h.   =  uJv.   ®  1   =  u.®  v.  =  u.®  v. 
ill  i  ill  1111 


H(U,V)  =  H(U,V)=(U,V) 

H^IGUI)  =  ( I©  Hi)  H3=(T«IV)  =  (I«JV) 

h2=(TT®iv)  =  (ii®iv)         h4=(1T®iii)  =  (ii«ni). 

Hl   *  H2  =hlh2h3    •••••V 

H2*Hl=hnhn-l    ■••••hl. 

But    h     /-,      =h    /-,    .    ,  .        by  Theorem    5.  6. 
n/2-r      n/2+l+r       ' 


Hl   *H2=hl  h2 


h    /->  h    /-,...  h,h  ! 
n/2     n/2  2     1 


H,   *  H,    =  h,   h,.  .....  o  h    ,-,     h     ._.  .  .  h,h.  . 

2  1  12  n/Z       n/Z  2    1 

Therefore     Hj  *  H2  =  H^*lTj"  =     H£  *  H,. 

H3*H4  =  *1*2     •••••     V 

H4*H3=yn-l *!• 

Rn/2-r=5in/2  +  l+r. 


Therefore     H-  *  H.  =  H„  *  H,    =  H„   *  H_. 
3  4  4  3  4  3 


13. 

14.  Table     5.  1     is  an  exhaustive  list  of  transformations  on  a 

complementary  pair  in  the   sequence  quadruple  form  with  the  Hamming 
vector  for  each  transformation. 

H,    *  H-,,    H?   *  G,  ,   H    *H.,    H   *H.,     and  their  complements  are  the 
only  possibilities  which  exist.      Therefore  only  four  Hamming  vectors 
are  possible  and  by  steps     9     and     13    they  occur  in  pairs  which  are 
complements  of  one  another. 
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TRANSFORMATION 


I           ( 

I  II 

in 

IV) 

cl 

(T  ii 

in 

IV) 

C2 

(i  ii 

In 

W) 

C 

[T  Ii 

TTi 

Tv) 

Al 

a  n 

In 

IV) 

Az 

[I  u 

in 

iv) 

Tl 

in  i 

in 

IV) 

T2            ( 

[i  ii 

IV 

in) 

T 

[ii  i 

IV 

in) 

Z 

in  "i 

in 

IV) 

Y 

JT  i 

Hi 

IV) 

X              1 

[TT  i 

in 

Iv) 

W             1 

Ln  1 

in 

W) 

V 

[T  ii 

IV 

TTi) 

U              I 

J  ii 

Tv 

in) 

S              1 

j  IT 

Tv 

ill) 

R             1 

J  "n 

IV 

ill) 

Q             1 

Ln  r 

IV 

ITT) 

P             1 

[IF  i 

Tv 

Hi) 

O             1 

Ji  1 

in 

IV) 

N             1 

,ii  i 

TTi 

TV) 

M             1 

J  ii 

IV 

ill) 

L              1 

Ji  1 

Til 

Iv) 

K              I 

[I  Ti 

Tv 

Hi) 

J               < 

JTT 

Tv 

HI) 

H             1 

[I   ii 

Tv 

TTi) 

G              1 

[ITT 

IV 

ill) 

F              1 

[I  Ii 

TTi 

IV) 

D              { 

[I  ii 

in 

iv") 

TT                       ' 

,n   i 

Tv 

III) 

B               I 

[TT  i 

IV 

HI) 

0              1 

Ji  1 

Tv 

in) 

HAMMING  VECTOR 
H(I)=H1*H2 

H<C2=Hl*H\, 
H(C)  =H1*H2 
H(A1)=H1*H2 
H(A2)=Hl*H2 

H(T2)=H3*H4 
H(T)  =H2*H1 
H(Z)  =H4*H3 
H(Y)  =H4*H3 
H(X)  =H4*H3 
H(W)  =4I4-H3 
H(V)  =H3*H4 
H(U)  =H3*H4 
H(S)  =H3*H4 
H(R)  =H3*H4 
H(Q)  =H2*H1 
H(P)  =H2*HX 
H(O)  =H4*H3 
H(N)  =H4*H3 
H(M)=H3*H4 
H(L)  =H4*H3 
H(K)  =H3*H4 
H(J)  =H2*H1 
H(H)  =H3*H4 
H(G)  =H2*H1 
H(F)  aB1*f!2 
H(D)  =^^2 
H(it)  =H2*H1 
H(B)=li2*H1 


COMMON  VECTOR 
H(D 
H(I) 
H(I) 
H(I) 
H(I) 
H(I) 
H(TX) 
H(TX) 

H(I) 

HlTj) 

H(TX) 

HtTj) 

H(TX) 

H(TX) 

H(TX) 

H(TX) 

H(TX) 

H(I) 

H(I) 

HfT^) 

H(TX) 

H(TX) 

H(TX) 

S(TX) 

H(I) 

H(TX) 

H(I) 

H(I) 

H(I) 

H(I) 

H(I) 

H(I) 


The    E     transformation  just  permutes  the  first  two  and  last  two  pairs 
of    (I    II    III     IV)     and  since  modulo     2    addition  is  commutative  no  change 
would  take  place  in  the     H    vector.      The  last  column  identifies  all 
Hamming  vectors  in  terms  of  two  vectors  and  their  complements. 

TABLE     5.  1 
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Theorem    5.12 

If  a  complementary  pair    (A,  B)    of  length    n    is  operated  upon  by  a 
complementary  pair    (C ,  D)    of  length    m    to  form    codes     U1}    U2    by  the 
time  sequence  exponential  method,   or  to  form  codes     Vi  ,    V2  by  the  inter- 
lace exponential  method  and  either    Ux  or    U2  is  time   reversed,    symbo- 
lized by    U    or  if    Vi  or    V2  is  time   reversed  symbolized  by    V    then 

H(U)    =  H(V)  =  H(A1,    B1  )2m    where  the   (A1  ,  B1)    pair  indicates   some 
general  transformation  of  the  original     (A,B)    pair. 

1.  A    ©    B     =     H(M)  .  where     B    indicates  the  time   reverse  of 

the  original    B 

2.  A    ©     B     =    H(H). 


3. 

B    ©     A    =    A    ©    B  =H(TX). 

4. 

B     ©    A    =    A    ©    B~H(L). 

5. 

Referring  to  Table     5.  1, 

HOxti-Hcr!). 

H(H)  =  H(T1), 
H(L)  =  H(T1). 

6. 

c,       c9                      c          d, 
Ui=A        A      A         B 

Bdm 

d        d       ,               d,       c 
U2=A   mA   T'1            A    l    B    m...m. 

^1 
..   B 

and 

c,       d,       c->                  c_      d 

V1=A    '   B    A  A   " .AmBm 

d  c  d        ,        d,       c ! 

.,,       .     m  „    m    .     m-  1     A     1    _     1 
V2=A  B  A  .....  A        B      . 

Time  reversing  one  of  the  pair  in  each,    for  example,   the  use 

of    T2    yields 

c ,         c0  c         d,  d 

tt  a     1      a     <-  a     rrin     1  t-,    m 

Uj.  =  A         A.........     A       B      ...c,,eB 

c  1         c->  c         d,  d 

U2=  B    l     B    2  B    mA    X......A   m 


and 


c,  d,      c,                          c        d 

Vi  =  A  B        A       ........   A        B 

c  1        d  1      c  9  c         d 

V2=B     l  A    l  B    2........      B    mA   m 
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c.  c. 

7.  H(U)  and    H(V)    consist  of  clusters  of    (A   *  9  B    1)    and 

d.  d. 

clusters  of    (B   1  9  A     1). 

8.  The  Hamming  vectors  for  these  clusters  have  the  following 
four  possible  forms: 

A©    B   ,   A  <B  B  ,     B  ©  A,      and  B    ®  A. 
These  four  possible  vectors  were  seen  in  steps     1,2,3,4,5    to  all  be 
equal  to    H(Ti).     If    Ui    and    Vi    had  been  time  reversed,    using  the 
operation    T1}   the   result  would  have  been     H  (Tx)a 

1  0.  Since  all  the  clusters  are  identical  and  there  is  a  cluster  for 

each  of  the     m    bits  in    C     and  each  of  the  m    bits  in    D, 
H(U)  =  HV)  =    [HfTi)]2,™    for    U2    or    V2    time   reversed  or 

in  general 

?  l       i 

H(U)  =  H(V)  =  HIA1^1)  where     (A  ,B  )  signifies  some  transfor- 

mation of  (A  ,B). 

>;c  >!<  Jfi 

Theorem    5.  13 

The  kernel  of  length  2  when  operated  on  by  a  complementary  pair 
by  either  the  time  sequence  exponential  or  interlace  exponential  method 
can  always  be  decomposed  at  least  once  by  the  interlace  method. 

1.  If  one  of  the  code  pair  formed  by  either  exponential  method 

is  time   reversed,     H(U)    =  H(V)  =  HfA^B1)  2m  by  Theorem  5.  14. 

2.  There  are  only  two  possible  Hamming  vectors  for  the 
quad,    01     and     10. 

3.  Therefore     H(U-)    =     H(V^  )  =  ( 01  )2m  or  ( 10)2m    both  of  which 
are  decomposable  by  the  interlace  method. 
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Theorem    5.14 

A  complementary  pair    (U1}    U2)    formed  by  the  time   sequence 
exponential  method  cannot  be  decomposed  by  the  time   sequence  method. 

1.  Let    (Ux  ,    U2  )    be  a  complementary  pair  of  length    Zmn  formed 
from  the  complementary  pair     (A,  B)     of  length    n    operated 
upon  exponentially  by  the  complementary  pair     (C,D)     of 
length    m,    formed  in  the  same  manner  as  in  Theorem     5.  3. 

2.  H(U)=(Cl   Odmftc.Od        ,)n...(c     ©dj^d^c      )n.  .  .  (d     lc,)n. 

1  Z      m- 1  m      11       m'  m      1 

3.  The  first    mn    bits  in    H(U)    are    m    clusters  composed  of    n 

zeros  or    n    ones  dependent  upon  the  bits  of    H(C,D). 

a  o      tu  c    i     tj    nmnimn  ,  mn  r.mn      ,         , 

4.  By  Theorem     5.  1     H=0        1  or     1  0        ,   for  decomposition 

of  a  code  of  length.  2mn    by  the  time  sequence  method. 

5.  This  requires  the  distance  of    (C,D)    be  zero  or    m,     but  by 
Theorem    4.  1     D(C,D)  =  m/2,    therefore     H(U)     is  not  decom- 
posable by  the  time   sequence  method. 

6.  If    Ui    or    U2    is  time   reversed     H(U)  =  H(A\  B1)2™    by 
Theorem    5.  1  2. 

7.  This  is  incompatible  with  the  form    H=0         1  since  it  would 
require  the  first    m/2     clusters  of    H(A1,B1)    to  have  a  distance 
of  zero  and  the  second    m/2     clusters  of    HfA^B1)    to  have  a 
distance  of    n.      Both  are  impossible  by  Theorem    4.  1     which 
requires  a  distance  of    n/2. 

#  *  * 

Theorem    5.15 

A  complementary  code  pair    (Vi,V2)    formed  by  the  interlace  ex- 
potential  method  (using  the   same  definitions  and  symbols  as  Theorem  5.  4) 

can  be  decomposed  by  the  time  sequence  method  only  if    H(C,D)= 

nm/2    m/2  ,1x1/2  _m/2. 

0  1  or     1  0 
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1.  H(V)  =  (c    «d     )2n(c.©d        ,)2n (c      ©d,)2n    by  Theorem  5.  4 

y     '        x    i      m  2      m- 1  m       1 

and  equation    5.  5. 

2.  The  clusters  of  zeros  and  ones  in  this  Hamming  vector  are  in 
lengths  of     2n    and  are  dependent  in  value  upon  the  bits  of  H(C,D), 

3.  H=0        1  is  the  requirement  for  decomposition  of  this  code 
by  the  time  sequence  method. 

I  7  I  7 

4.  If    H(C  ,D)  =  0m       1  or  its  complement,   it  will  satisfy  step  3 

/flZn,m/2   .  ,Zn.\xr\IZ    ^mn,  mn    TT/-,r\  j   .  .      ->  £. 

since     (0      )  (1       )  =0        1        =H(V),   and  it  also  satisfies 

Theorem    4.  1 . 

5.  If    Vi    or    V2    is  time  reversed     H(V)  =H(A1,  B1)2"1    by 
Theorem     5.12. 

6.  This  is  incompatible  with  the  form    H=0        1  since  it  would 
require  the  first    m/2     clusters  of    H(AX,  B1  )     to  have  a  distance 
of  zero,   and  the   second  m/2     clusters  to  have  a  distance  of    n. 
Both  are  impossible  by  Theorem    4.  1     since  it  requires  a 
distance  of    n/2. 

7.  Step    4     showed  that    H(C ,  D)  =  0m  /2lmZ.      This   satisfies  the  re- 
quirement for  the  time   sequence  decomposition  by  inspection. 
Any  other  arrangement  of    H(Ci,D)  would  fail  to  allow  decom- 
position. 

#  >!<  >'f 

Theorem    5„  1  6 

Any  complementary  sequence  pair  (Vi,  Vs)  formed  by  the  interlace 
exponential  method  as  in  Theorem  504  (using  the  same  definitions  and 
symbols),    the  quad  as  the     (C,D)  pair,    is  time   sequence  decomposable. 

1.  The  quad  has  only  two  possible  Hamming  vectors,      01  and     10. 

I  7  /  7 

2.  This   satisfies  the  form    H=0  1  where     m=2,    the  quad 
length. 

3.  Vi    and    V2  are  therefore  time  sequence  decomposable  by 
Theorem    5.  1  5. 

*  *       7R  * 


Theorem    5.17 

A  complementary  pair  formed  by  either  exponential  method  is  de- 
n/2 
composable  by  the  interlace  method  if  and  only  if    H(A1,  B1)  =  (01. )  or 

n  /  Z 
(10)         .    All  definitions  and  symbols  are  assumed  to  be  the   same  as 

those  in  Theorems     5.  3,  5.  4    and    5.  12. 

1.  By  Theorem     5.  2    these  codes  are  decomposable  by  the  inter- 
lace method  only  if  they  are  of  the  form    H=(01)         or  (10) 

2.  H(U)    and    H(V)    are  clusters  of  at  least    n    ones  or    n    zeros 
and  therefore  cannot  have  the  form    (01) 

3.  H(U)  and    H(V)  =  H(A\  B1)2™    by  Theorem    5.12. 

4.  If  and  only  if    H^.B1)    has  the  form     (01)n/2    or     (10)n^2can 
H(U)    or    H(V)=    [(01)n/2]    2m  =  (01)nm,   or  using     10    instead 
of    01,    (10)        .      This  satisfies  the  requirement  of  Theorem 
5.2    that  the  Hamming  vector  be  of  the  form  (01)       or     (10)    . 


A  brief  description  of  the  method  used  to  prove  that  codes  of  length 
200    formed  from  kernels  of  length     10,    cannot  be  decomposed  into  ker- 
nels of  length    50,    will  now  be  given  as  a  guide  to  the  actual  proof  of  the 
theorem.      The  method  used  is      to      exhaust  all  possible  methods  of  de- 
composition by  examining  all  possible  Hamming  vectors  for  the     (ASB) 
and    (C,D)  pairs.      Table     5.2     lists  all  possible  Hamming  vectors  for 
kernels  of  length     10. 

Kernel     1  Kernel    2 

A=1000000110  A=1001101111 

B=1001010001  B=1100001010 
The    4  possible  Hamming  vectors.        The  4  possible  Hamming  vectors, 

0001010111  0101100101 

1110101000  1010011010 

0000101111  1100101100 

1111010000  0011010011 

All  possible  Hamming  vectors  of  kernels  of  length     10„ 

TABLE     5.2 
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The  two  methods  of  generating  codes  of  length    200    from  kernels 
of  length     10    are  the  time  sequence  exponential  method  and  the  interlace 
exponential  method  given  respectively  by 

U1=AClAC2....ACl0Bdl......Bd10 

d10Ad9         ,dljl0  Jl  (5°6) 


U2=A        A    '.  ..A  B        .....  .B 

cl        dl        c2  d2                      c10    d10 

V1=A         B_      A  L  B      .......A        B    iU 

d10      c10    d9  c9                     d10    cl 

V2=A          B      UA   y  B      .......  A    1UB 


(5.7) 


where     (A,B)    and     (C,D)    are  kernels  of  length     10. 

The  Hamming  vectors  of  these  codes  of  length    200  are  anti- sym- 
metric about  their  centers  and  are  in  clusters  of     10    zeros  or     10    ones 
for    H(U),   and  in  clusters  of    20     zeros  or    20     ones  for  the     H(V)     vector. 
Time  reversing  either    Ui    or    U2    and  similarly  either    Vx    or    V2    , 
gives     20    clusters  of  Hamming  vectors  of  the  forms  listed  in  Table  5.2, 
since  each  of  the     20     A    codes  is  now  matched  with  some  transform  of  B. 

The  known  possible  ways  for  kernels  of  length    50    to  generate  com- 
posite codes  of  length    20    are: 

1.  Either  interlace  or  time  sequence  kernels  of  length    50    to 
form  composite  codes  of  length     100    and  then  either  inter- 
lace or  time   sequence  these     n=100    codes  to  form  composite 
codes  of  length    200. 

2.  Take  kernels  of  length    50    and  either  time  sequence  expo- 
nential or  interlace  exponential  with  the  quad  to  form  codes 
of  length    200. 

3.  Take  kernels  of  length    2     and  either  time  sequence  exponen- 
tial or  interlace  exponential  with  kernels  of  length     50. 

An  examination  of  all  possible  Hamming  vectors  of  length    200    from 

each  of  these  methods  will  show  the  incompatability  of  these  Hamming 

vectors  with  those  generated  by  the  kernels  of  length     10. 

Theorem    5.18 

No  codes  of  length    200    formed  from  kernels  of  length     10    can  be 

decomposed  into  kernels  of  length     50    by  standard  methods  of  decomposi- 
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1.  Given    Ux    U2    and    Vi     V2    both  complementary  pairs  of  length 
formed  by  the  time   sequence  exponential  and  interlace  expo- 
nential methods  as  given  in  equations     5.  6    and    5.  7  respectively. 

2.  The  time  sequence  method  of  decomposition  is  impossible  for 
each  of  the  following  Hamming  vectors  for  the  following  reasons. 

a.  H(U)     can  never  be  decomposed  by  time  sequence. 
Theorem    5.  14. 

5    5 

b.  H(V)     requires  a  form    0    1       by  Theorem     5.15. 

Table     5.2    contains  no  vector  of  this  form. 

c.  H(V  )     and    H(U)      can  never  be  decomposed  by  time 
sequence.      Theorems     5.14,15. 

3.  The  interlace   sequence  method  of  decomposition  is  impossible 
for  each  of  the  vector  forms  for  the  following  reasons. 

a.  H(U)    and    HfV)    can  never  be  decomposed  by  interlace. 
Theorem    5.  1  7. 

b.  H(V)      and    H(U  )     are  decomposable  by  interlace  only  if 
H(A1,B1)is  of  theform    (01)5    or    (10)5    by  Theorem    5.  1 7. 
Table  5.  2    contains  no  such  vectors. 

4.  Assume  that  one  of  the  Hamming  vectors    H(U)    or    H(V)    has  a 
second  exponential  form,   first  with  the  kernel  of  length    2 
operated  on  exponentially  by  the  kernel  of  length    50,   and 
secondly  by  the  kernel  of  length    50     operated  on  exponentially 
by  the  quad.      Let    H1  be  the  designator  for  the  Hamming 
vectors  of  the     2    and     50     combinations. 

a.  Consider  first  the  quad  acted  on  exponentially  by  the 
kernel  of  length    50o      H^U)=H1(V)=  H^A^B1).     But  since 
A^B1  is  the  quad,   by  Theorem    5.  13,    it  is  decomposable 
by  the  interlace  scheme.      This  is  not  possible   since 
step     3    above  exhausted  all  possible  interlace  decompo- 
sitions. 

b.  Next  consider  the  kernel  of  length     50    acted  on  exponen- 
tially by  the  quad.      Considering  first  the  exponential 
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interlace  scheme,   by  Theorem     5.  15     if    Hi(C,D)=01  or     10, 
which  it  does  for  the  quad.     Therefore  the  code  is  time   se- 
quence decomposable.     However,   in  step    2    all  of  these  time 
sequence  possibilities  were  exhausted  and  therefore  this  de- 
composition is  not  possible.     H1(U)    would  have  a  form  con- 
taining either     50    zeros  or     50     ones  in  clusters„      For  the 

original     10     length  codes  to  form  these  would  require    a 

5   5 
0    1       Hamming  vector  which  is  not  listed  in  Table     5.  2. 

This  exhausts  all  possible  Hamming  vectors  that  might  lead  to  de- 
composition,   therefore  codes  of  length    200    formed  from  kernels  of 
length     10    cannot  be  decomposed  into  kernels  of  length    50    by  standard 

methods. 

*  *  # 

Table     5.  3     shows  the  possible  Hamming  vectors  for  codes  of 

length    26. 

A=01 00 1101 11 10 101 11 1001 110 10 

B= 10 11 001 00001  111  111001 11010 

The  4  possible  Hamming  vectors. 

11111111111101000000000000 

00000000000010111111111111 

11101110111010100010  001000 

00010001000101011101110111 

All  possible  Hamming  vectors  of  the  kernel 
of  length  26. 

Table  5.  3 

An  examination  of  Table     5.  3    which  is  exhaustive  for  the  Hamming 

vectors  of  length    26  and  Table     5.2    for     10,    shows  that  the  next  two 

theorems  can  be  proved  by  identical  methods  to  those  used  in  Theorem 

5.  18.      Since  these  proofs  exactly  parallel  that  of  Theorem     5.  18    they 

will  not  be  given. 

Theorem    5.  19 

Complementary  pairs  of  length     520    formed  from  the  time  sequence 
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exponential  or  interlace  exponential  from  kernels  of  length    26    and     10 
cannot  be  decomposed  by  standard  methods  into  kernels  of  length     130. 

$  >jc  3jt 

Theorem    5.  20 

Complementary  pairs  of  length     1.352     formed  from  the  time   sequence 
exponential  or  interlace  exponential  method  from  kernels  of  length    26 
cannot  be  decomposed  by  standard  methods  into  kernels  of  length     338. 

*  *  # 

The  study  of  Hamming  vectors  thus  far  has  disclosed  that  a  pattern 

/.._1.n/2  _n/21n/2  ..  ,  .        .  ,  ,         .       _ . 

oi     (01)  or    0  1  or  their  complements  in  either  the     (A,B) 

pair  or     (C,D)  pair,   is  required  in  order  that  a  code  formed  from  one 
kernel  length  by  an  exponential  method  be  decomposed  into  a  pair  from 
a  different  kernel  length.      The  quad  is  only  kernel  which  thus  far  satis- 
fies this  criteria.     However,    the  other  kernels  come  close,    as  a  check 

i* 

4         4  2  2 

shows  that    0    101       and     (01)    10(01)       are  among  the  Hamming  vectors  of 

12         12 

length     10,    and  0       101        is  a  Hamming  vector  of  length    2  6. 

If  the   quad  is  involved  in  the  formation  of  composite  codes  it  makes 
their  decomposition  more  likely.      This  is   shown  in  Theorem     5.  12     and 
Theorem    5.  16.      Appendix     II    demonstrates  this  fact  quite  clearly  since 
it  contains  all  codes  of  length     16,    less   some  operational  redundancies. 
Each  of  the  code  pairs  and  its  time  inverse  is  decomposable  by  some 
method  into   shorter  length  codes.      Appendix     III    lists  all  the  pairs  of 
length    20,    less  some  operational  redundancies.     Each  code  pair  or  its 
time  inverse  can  be  broken  into  its  originating  pair  through  the  use  of 
Hamming  vectors,   but  no  code  and  its  time  inverse  can  both  be  decom- 
posed by  the   standard  Hamming  vector  methods. 

The  writer  would  like  to  end  this  chapter  with  a  conjecture  that 
kernels  of  length     50     do  not  exist.      This  conjecture  is  based  upon  three 
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bits  of  evidence,    the  proof  of  Theorem    5„  18    being  the  first.      The  second 

2      2      2      2 
piece  of  evidence  is  based  upon  the  fact  that     50=7   +1    =5    +5    .      From  this 

one  might  assume  that  codes  of    200  formed  from     10's     might  cover  just 
one   set  of  the  ones  determined  by  equation    20  5.     This  proved  not  to  be 
true  however,    since  the  number  of  ones  in  both  the  exponential  formations 
of  codes  of  length     200    from  kernels  of  length     10,    if  they  were  decompo- 
sable into  kernels  of  length    50,   would  exhaust  all  possible  unordered 
pairs  of  ones  for    n=50    as  determined  by  equation     2„  5„      The  last  bit  of 
evidence  to  base  this  conjecture  upon  is  that  the  code  length  is  not  twice 
a  prime  number.      The  only  previous  kernel  which  might  have  existed,   but 
was  not  twice  a  prime  number,   was    n=18,   and  it  did  not  exist.     Similar 
reasoning  extends  this  conjecture  to  kernels  of  length     130    and    338. 
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CHAPTER    VI 
SUPPLEMENTARY  AND    CYCLIC     COMPLEMENTARY    CODES 

This  chapter  is  concerned  with  two  classes  of  codes.      The  first  of 
these  classes  consists  of  quadruples  of  sequences  with  the  property  that 

the  total  number  of  likes  at  each  spacing  equals  the  total  number  of  un- 

2 
likes  at  the  same  spacing.      These  are  called  supplementary  codes.        The 

second  class  of  codes,   which  will  be  defined  later,    consists  of  the  cyclic 
complementary  codes.     Both  of  these  code  types  have  complementary 
sequence  pairs  as  a  possible  subset;     this  allows  the  use  of  their  prop- 
erties as   screens  in  the   search  for  new  kernels. 

Supplementary  codes  will  be  discussed  first,    since  the     n=26     search 
involved  the  use  of  this  property  as  a  screen,   while  the  cyclic  comple- 
mentary property  was  not  used  until  the  search  for  kernels  of  length    34. 
Only  one  theorem  will  be  proved  for  each  of  these  types  of  codes,    since 
these  two  theorems  were  the  only  ones  applied  in  the  actual  search  for 
new  kernels. 

If  u,  U-.U- u 

12    3  n 

V  .  V-x  V  V 

1    c    5  n 

w    w.~w  w 

1^3  n 

x.  x-.x_.  ......  x    are  sequences  of  zeros  and  ones  and 

12    3  n  ^ 

satisfy  the  constraint  of  equation    60  1,    they  form  a  quadruple  which  is 
supplementary. 

£     (u.0u  .«1)  +  (v.0v  .01)+w.«w    ,       .91)  +  xJx  .§1)  = 

v    i      n+i-j      '       x    i      n+i-j      '  l       n+i-j  i      n+i-j 

i=l 

J 

E    (u.§u  .)  +  v.9v    ,.     .)  +  (wJw    ,.     .)  +  x.Qx    ,.     .)  =  2j  (6.  1) 

._    v    i      n+i-j'  i      n+l-j'       v     i       n+i-j'  i      n+i-j'      J  \    •     / 

for  all  j,    1  _<  j_<  n  -    1.     Note  that  there  are     2j    like  pairs  and    2j    unlike 

pairs  at  a  spacing  of    n-j. 
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For  simplicity  in  the  proof,    we  use  the  symbolism  that  was  used  in 

Chapter  2, 

J 

let  LTT  =    E    u.  9  u  .©1  for  all  j,    1  <  j  <  n- 1 

U  ,      i         n+i-i  J        —  j  — 

i=  1  J 

j 
UTT  =    £    u.  0  u    .  .    • 
i=l      X  J 

A  restatement  of  equation    6.  1     in  this   symbolism  is 

Lu  +  Lv  +  Lw  +  Lx  =  uu  +  uv  +  uw  +  ux  (6«  2) 

Theorem    6.  1 

Any  complementary  sequence  pair,    (A,B),    written  in  standard 
(I    II    III    IV)    form  has  the  property  that    I,  II,    III    and    IV    are  a 
supplementary  quadruple. 

1.  L    =UR  and  U.=.LR  by  the  definition  of  complementary. 

2.  Considering  just  the  even  values  of    j     rather  than  all  values, 

LA=LT  +  LII» 

uA=uI  +  un, 

LB=1jIII  +  LIV 

UB=UIII+UIV 

3.  Lj  +  Ln  =Um  +  Uw. 

4-       Lm+Liv=ui+Uir 

5.  Adding  the  equations  in  steps     3    and    4    gives 

Li  +  Ln  +  Lm  +  Liv  =  ui  +  Un  +  um  +  urv  • 

6.  This  is  the   same  as  equation    6.  2    and    I,    II,    III,    IV    formed 
in  the  sequence  quadruple  form  from  a  complementary  pair 
are  supplementary. 

Since  the  supplementary  property  applies  only  to  the  even  spacings 
in  a  code  pair,    the  converse,    that  all  interlaced  supplementary  quad- 
ruples are  complementary,   is  not  true. 

*  *  * 
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The  cyclic  complementary  property  is  useful  in  itself  for  a  communi- 
cations or  telemetering   system;     however,   at  this  time  our  concern  is 
for  its  property  as  a  necessary  condition  for  a  code  pair  to  be  complemen- 
tary.    A  cyclic   sequence  or  code,   as  the  name  suggests,    is  a  never 
ending   sequence  of  zeros  and  ones  which  has  a  period  of  n  bits.  A  cyclic 
complementary  sequence  pair,   is  a  pair  of  cyclic  codes,    each  of  period 
n,   where  the  number  of  likes  of  one   sequence  equals  the  number  of  un- 
likes  of  the  other  sequence  for    n    possible  matches  for  all  spacings 

from     1     to    n-1,    or  stated  in  terms  of  the  bits  of  the    A    and    B  codes 

n  i=n 

C.  =      S    a.®  a  =    S    b.®        L.®  1         1  <  i  <  n- 1 .  (6.3) 

J        i=1     i       n~J+1         i=1     i   n-j+i  -  J  - 

where    a     ,  .  =  a.     and    b     ..   =  b.     since  the  period  is    n. 

n+i  i  n+i     '    l  c 

Theorem    6.  2 

A  complementary  sequence  pair    (A,  B)    of  length    n,   if  written  in  a 
cyclic  fashion,    is  always  a  cyclic  complementary  pair  of  period    n. 

1.  Since     (A,  B)    are  a  complementary  pair  they  satisfy 
equation    2.  3. 

j  i=j 

f .  =    E   a.®  a        ,.      =    S      b.®  b       ._,_.©    1      fori,    l  =  j=n-l. 
J     i  =  1      i       n-j+i         i=1       i       n-j+i 

2.  Expanding  for    f,     gives 

a,®  a     =  b,®  b   ®  1     which  is  equal  to    a   ®  a,    =  b   ®  b,®  1. 
lnln  nlnl 

3.  Expanding  for    f      ,     gives 

(ax®  a2)  +  (a2®  a3)  +  a3®  a4)  + (an_  }®  an)=(b1®b2®l)  + 

(b2®  b3®  1 )  +  (b3®  b4®  1 )  +.  t +(bn_  !®  bn®  1 ). 

4.  Adding  the  equations  of  step     2     to  step     3     gives    C,„ 

5.  Similarly  expanding    f?     gives 

(aL®  an_1)  +  (a2®  an)  =  (b^  bn-1®  1)  +  (b2®  bn®  1) 
which  is  equal  to 

(anl®  ax)  +  (an®  a2)  =  (bn-1®  b^  1)  +  (bn®  b£  ®    1). 
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6.  Expanding     f      _,     gives 

(a1«a3)  +  (a2®a4)  +  (a3§a5)  +  .  .  .  (an_2»an)  =  (b^b^  1)  + 

(b.©b.®l)  +  (b-^b^l)  + +  (b      9®b    01) 

c      4  3       b  n-Z       n 

7.  Adding  the  equations  of  step     5    to  step    6     gives     C->. 

8.  Continuing  this  same  procedure  through  all  possible  values  of 
j    would  show  that    C.  =  f .  +  f       ..     Therefore  the  cyclic  com- 
plementary  constraint  equations  are  based  on  the   sum  of  two 
restricted  portions  of  the  complementary  constraint  equations. 
Therefore  all  complementary  pairs  have  to  be  cyclically 
complementary.      The  converse  is  not  nece  ssary  true  since 

C  .    is  a  sum,    and  in  a  sum  the  addend  and  augend  are  not 
unique. 

#  #  * 

Since    C.  =  f.  +  f       ■,    each  constraint  of  the  complementary  property 

J  J  J 

is  used  twice  as     j    varies  from     1     to    n-  1 ;     therefore     C.     is  symmetric 
and  is  centered  at    C     /?.      In  the  application  of  the  cyclic  complementary 
property  in  Chapter     8,    half  of  the  characteristic  cyclic  number  is  de- 
leted due  to  this  symmetry  property^ 

Both  of  the   screens  described  in  this  chapter  were  suggested  by 

14    19 
Dr.    Golay;  however,    to  the  best  of  my  knowledge  the    theorems 

and  their  proofs  have  not  appeared  in  print. 
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AN    EXHAUSTIVE    SEARCH    FOR    KERNELS    OF    LENGTH    26 

CHAPTER    VII 

In  his  paper  Golay  emphasized  the  need  for  an  exhaustive  search  for 
complementary  code  pairs  of  length    26.      This  chapter  describes   such  a 
search.      The  purpose  of  this  search  was  to  determine  if  there  were  any 
kernels  of  length    26     and  if  any  of  these  exist,   the  number  of  such 
kernels.      It  would  prove  very  difficult  to  surmise  the  number  of  kernels 
in  a  code  of  this  length  predicated  on  the  known  kernels  of  other  lengths, 
since    n=2    had  one  kernel,    n=10     had  two,   and    n=18    had  zero.     An  ex- 
haustive search  for  kernels  of  length     26     might  throw  some  light  on  a 
possible  general  method  for  finding  kernels  of  longer  lengths,    or  might 
disclose  some  sort  of  pattern  showing  the  distribution  of  kernels  among 
the  possible  code  lengths. 

The  only  feasible  method  to  accomplish  this  search  was  with  a  high 
speed  digital  computer.      This  chapter  is  concerned  with  the  computer 

program  and  the  results  of  the  exhaustive   search  for  kernels  of  length  26. 

52 
At  first  glance  it  would  appear  that  there  are     2         possibilities  to 

be  screened  for  an  exhaustive  search  for  all  possible  kernels  of  length 
2  6.  Although  this  is  true,  the  application  of  some  of  the  theorems  de- 
veloped in  the  earlier  chapters  immediately  eliminates  from  considera- 

52 
tion  large  blocks  of  the     2  possibilities  which  are  either  redundant  or 

impossible.      This  is  a  very  necessary  procedure  since  even  with  the 
highest  speed  computers  of  today  it  would  be  impossible  to  investigate 
this  number  of  possibilities  in  a  life  time. 

The  program  was  developed  for  the  CDC    1604;    this  computer  ac- 
complishes approximately    200,000     operations  per  second.      This   speed 
coupled  with  the  powerful  screens  used  in  programming  made  possible  a 
reduction  of  the  computer  run-time  to  approximately    75    hours.     A  brief 

outline  of  the  various   screens  used  and  their  reduction  factor  will  now  be 
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given  before  continuing  in  more  detail  with  the  actual  programming 
techniques  used.      The  word  code  or  code  pair  as  used  in  this  chapter 
and  also  in  the  next  chapter  will  be  understood  as  possible  complemen- 
tary code  pair. 

The  first  screen  utilized  was  the  number  of  ones  which  must  appear 

2  2 

in  each  of  the  code  pairs.      This  was  given  by    n=(n-p-q)      +  (p-q) 

as  derived  in  chapter    2,   where    n    is  the  length  of  the  code,     p    is  the 

weight  of  the    A    code  and    q    is  the  weight  of  the     B    code.     The  possible 

solutions  for    n=26    were  the  unordered  weight  pairs     (16,15),    (16,11), 

(11,  10),    (15,  10).      The  pair  (p=l6,   q=15)    was  arbitrarily  chosen  to  be 

used  in  the  program,      This   reduced  the  total  number  of  possible  code 

-  952    Q    in15     ,  .       26!  261  _    A    in13 

pairs  from    2     ^8x10         down  to    .  ,  ,  .  , , —  x     nrrxT1     s*    4x10       .     This 

was  a  reduction  by  a  factor  of    200,    which  brought  the  life  time  search 

down  to  a  little  less  than  a  year,   but  this  amount  of  time  was  still  not 

feasible  for  a  computer  search.     Golay  remarked  in  a  footnote  that  a 

complementary  pair  might  be  thought  of  as  being  composed  of  two  inter- 

9 
laced  half  length  codes,    in  this  case    n=13.    '       Since  it  was  obvious  the 

total  number  of  ones  in  both  of  the  interlaced  pair  must  equal  the  number 
of  ones  in  one  of  the  complementary  code  pair,    equation     2.  5    was 
applied  again,   this  time  for    n=130      The  unordered  weight  possibilities 
for    n=13  were     (9,  7),  (9,  6),   (6,4),    (7,4).      In  order  to  be  compatible  with 
the  (16,  15)  pair  for    n=26,    the   (9,  7)  pair  was   selected  for  the  code  con- 
taining    16    ones  and  the   (9,  6)  pair  was  chosen  for  the  other  code  which 
contained     15     ones. 

The  basic  problem  had  now  been  reduced  to  all  possible  combinations 

of  the  four  thirteen  bit  sequence  families.      This  was  equal  to 

13'  13'  13'  13' 

q,    .'      x     Q,   '  |      x    tt-tt — •  x    7»  1>"  =     715x715x1716x1716    or  approximately 

12 

1.5x10         possibilities.     Although  this   reduced  the  problem  by  one  more 
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order  of  magnitude  it  was  still  much  too  large  for  a  practical  search. 
Dr.   Golay  in  a  private  communication  with  the  author  pointed  out  the 

possibility  of  using  the  supplementary  characteristic  as  an  additional 

14 
screen.  Theorem     6.  1     gives  the  necessary  and  sufficient  condition 

that   complementary  codes   sequence  quadruple  form  be  supplementary 
sequence  quadruples.     The  codes  of  length     13    were  therefore  categor- 
ized according  to  their  first     3    and  last     3    bits.      The  number  of  like 
pairs  at  spacing  twelve,    at  spacing  eleven,   and  at  spacing  ten  was  com- 
puted for  each  of  the     13    bit  sequences.      For  example,    0111100111011 
has     0     likes  for  spacing  twelve,      1     like  for  spacing  eleven,   and    3     likes 
for  spacing  ten.      The  block  number  given  to    0111100111011     was  there- 
fore    013.      Out  of  sixty  four  possible  combinations  only  eighteen  different 
block  numbers  were  generated  for  each  of  the  two'  number  sequences  of 
all  possibilities  of  thirteen  bits  with  nine  ones  and  with  seven  ones.     An 
examination  of  Figure    7.  3  reveals  that  there  are  a  widely  varying  number 
of  members  within  each  of  the  blocks.      The   supplementary  characteristic 
states  that  the  total  number  of  like  pairs  for  each  spacing  must  equal  the 
total  number  of  unlike  pairs  for  the   same   spacing.     With  four  codes  at  a 
spacing  of  twelve  there  are    4    possibilities,    therefore  two  of  them  must 
be  likes;     at  a  spacing  of  eleven  there  are  eight  possibilities,    therefore 
four  of  them  must  be  likes;    and  likewise  at  a  spacing  of  ten  there  are 
twelve  possibilities  of  which  six  must  be  likes  to  satisfy  the   supplemen- 
tary property.      Block  numbers  were  added  for  each  possible  combination 
to  see  if  the  total  were     246,    and  in  cases  where  the  total  was     246    those 
particular  codes  within  these  blocks  were  then  sent  to  the  next  necessary 
but  not  sufficient  condition  for  the  codes  to  be  complementary.      If  the 
total  was  not     246,   the  blocks  of  codes  were  rejected.,      The  power  of  this 
screen  is   shown  in  the  example  where  say  two  of  the  blocks  had  code 
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numbers     000,    000.      The  only  possible  combination  to  add  to     246 
would  be     123,  123    for  the  other  two  codes.      This  cuts  the  total  number 

of  possibilities  in  this  particular  block  combination  down  from 

106 
1716x1716  =  2.  9x        to  only    64     possibilities.     Although  this  admittedly 

was  the  most  extreme  example,   this  screen  reduced  the  number  of 

q 
possible  codes  checked  to  about    2x10       or  further  reduction  of  three 

orders  of  magnitude. 

The  next  necessary  condition  to  be  checked  was  the  parity  test, 
equation    2.  4.     All  code  quadruples  which  satisfied  the  parity  test  were 
then  sent  to  the  necessary  and  sufficient  like  pair,    unlike  pair  check  for 
each  spacing.      The  like,   unlike  check  subroutine  was  quite  long  and  also 
involved  word  unpacking,,      This  subroutine  which  was  ignored  in  the 
calculation  of  run  time  increased  the  computer  search  time  from  a  calcu- 
lated    70     hours  to  an  actual  run  time  of  the  order  of    75     hours. 

The  half  length  codes  in  the  program  were  interlaced  in  the  same 
manner  as  was  done  in  the  operations  group  formulation,,     This  method 
of  combination  made  the  parity  check  easy  to  calculate  since  the  bit 
positions  in  the  computer  words  were  the   same  as  the   sequence  posi- 
tions of  the  codes.     Modulo     2    addition  was  directly  applied  to  the  code 
quadruple  to  form  a  vector  of  all  ones,    providing  the  parity  test  held. 
If  the  code  vector  was  not  all  ones  the  code  quadruple  was  rejected. 

It  is  to  be  noted  that  by  selecting  the  ordered    (16,  15)    pair  for  the 
number  of  ones  in  the    A    and    B    code,   the  operations  Ci,    C2    ,   and    C 
were  eliminated  from  possible  redundant  consideration.      These  were 
eliminated  by  their  requirement  for    (10,15),    (16,11),     and    (10,11) 
ordered  pairs  of  ones  respectively,   based  on  equation    2.  5.      The  inter- 
lacing of  the  half  length  codes  with  ordered  pairs  of     (9,  7)    ones  and 

(9,6)     ones  respectively  eliminates  from  consideration  the  operations 
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Ti  ,    T2  ,    T     and    Ai  since  these  would  have  required    (7,9,)  (9,6), 
(9,7)(6,9),    (7,9)(6,9)    and    (4,7)(4,6)     ones  respectively.      The  only  basic 
transformation  not  deleted  by  controlling  the  weight  in  each  of  the  code 
quadruples  was    As.      This  is  seen  to  be  true  because  if  the  pair 
(9,7)(9,6)    has  the  operation    A2    performed  on  it,    the   result  is  (9,6)(9,7). 
If  this  code  is  now  transformed  by  the  trivial  exchange  operation,    E,   a 
code  pair  of  the  original  form  is  generated.      The  scheme  which  was  used 
to  eliminate  the    A2    redundancy  will  be  discussed  late  in  this  chapter. 

Figure     7.  1     is  a  rough  flow  diagram  for  the  program  used  in  the 
exhaustive  search  for  kernels  of  length    26.     With  modifications  in  data 
as  necessary,    this  program  was  also  used  for  exhaustive  searches  of 
composite  codes  of     16    and    20. 

A  brief  description  of  the  blocks  in  Figure     7.  1  will  give  the  reader 
some  familiarity  with  the  programming  philosophy  used  in  this   search. 
The  actual  programs  as  written  for  the  CDC    1604  are  given  in  Appendix  V. 

The  number  generation  subroutine  is  able  to  generate  any  number  up 
to    48     bits  in  length  (word  size  of  the   1604)  with  any  number  of  ones  up 
to  a  maximum  of    48.      This  is  accomplished  in  the    1604  by  inserting  a 
number  of  the  form    00000.  ..  00001  1 1 1.  .  0  1 1 1     into  the  accumulator, 
where  the  total  length  of  the  number  is     48    bits  and  it  contains     p    ones 
and  has  a  code  length  of    nc      This  number  is  then  shifted  left    48-n  bits. 
This  shifted  result  is  then  checked  to  see  if  it  is  negative,   which  indicates 
a  one  in  the  most  significant  bit  position.      A  procedure  is  now  set  up  to 
count  a  one  or  zero  and  then  shift  left  one  bit  and  repeat  the  process  for 
the  entire    n    bits.      If    p  ones  are  counted  the  number  is  stored,    if  not  it 
is  rejected.     A  number  one  larger  than  its  predecessor  is  inserted  into 
the   system  and  the  process  is   repeated  until  the  ones  are  all  in  the 
uppermost  bit  positions.  q  . 


Blocks     2,3,4    are  used  to  compute  a   supplementary  characteristic 
number  for  each  of  the     715    thirteen  bit  numbers  with    9    ones  and  also 
for  each  of  the     1716     thirteen  bit  code  numbers  with     7    ones;     then  in 
accordance  with  their  supplementary  numbers  these  are   stored  into 
group  blocks.      The  supplementary  number  was  derived  by  using  a  mask 
to  expose  for  consideration  only  the  first     3    and  last     3    bits  of  each  of 
the     13    bit  code  numbers.      These  masked  numbers  were  then  compared 
with  an  exhaustive  list  of  all  possible  combinations  of  zeros  and  ones  in 
these     6    bits.      These     64    possible  combinations  were  in  an    M    to  one 
correspondence  with  a  list  of     18    possible  supplementary  characteristic 
numbers   (see  Figure     7.2).     This  correspondence  allowed  the  temporary 
attachment  of  this   supplementary  number  to  each  of  the  codes.      The  code 
numbers  were  rearranged  into  blocks  according  to  supplementary  number 
and  at  the  same  time  had  their  temporarily  attached  supplementary 
number  deleted.      This  part  of  the  program  is  purely  for  putting  the  data 
into  a  useable   state  and  is  used  only  once  per  run,      The  rest  of  the  pro- 
gram is  highly  iterative. 

Consideration  of  the     (I    II    III    IV)    form  of  the  possible     26    bit 
codes,    shows  that  each  of  the     18    blocks  of    I,   II,    III,  and    IV     must  be 
compared  against  each  other  indicating     1 8x1 8x1 8x1 8 ~ 1 05, 000  possible 
block  comparisons.      This  was  cut  down  somewhat  by  use  of  a  slightly 
different  technique.      This  was  accomplished  by  loading  the  computer 
accumulator  with    246    and  from  this  subtracting  the  first  block  number 
of    I,   then  subtracting  the  first  block  number  of    III,   and  then  subtrac- 
ting the  first  block  number  of     II.     Rather  than  subtracting  the  Block 
number  of    IV,    the  difference  obtained  was  instead  checked  against  all 
all  possible  block  numbers  of     IV    to  see  if  it  was  listed0     Since  if  the 

block  number  did  exist  only  one  could  exists      This  procedure  cut  the 
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possible  number  of  block  comparisons  by  a  factor  of  eighteen  and  left 
only    5800     combinations  to  examine   rather  than     105,000,      If  the  block 
did  exist  all  the  codes  within  this     (I    II    III    IV)     block  grouping  were 
then  given  the  parity  test.     If  the     IV    block  number  did  not  exist    II    was 
stepped  ahead  one  block  number  and  the  process  repeated.     After    II 
had  cycled  through  its     18    block  numbers,     III    was  then  stepped  ahead 
one  block  number  and  the  entire  process   repeated.     Similarly  when    III 
had  cycled  its     18    block  numbers,     I    was   stepped  ahead  and  the  process 
repeated.      To  avoid  the     As  operation  redundancy,    whenever     I    was 
stepped  up  one  block  number    III    was  started  from  this  same  block 
number  rather  than  at  the  first  block  number.      This  avoided  the  block 
number  combinations     (I    III) ,  and    (III    I)     both  appearing  except  when 
block    I    was  equal  to  block    III.      This  block  number  cycling  allowed  a 
vantage  breakpoint  for  partitioning  the  program  into  suitable  size  in- 
crements for  computer  run  times. 

The  parity  check  was  made  by  adding(Modulo     2)     a  code  from    I 
to  a  code  from     III    to  a  code  from    II.      This   resultant  was  then  checked 
to  see  if  it  existed  as  a  code  in  the  list  from  the     IV    block.      If  it  did 
exist  it  was  complemented  and  then  sent  to  be  unpacked  for  the  necessary 
and  sufficient  like,   unlike  test.      If  it  did  not  exist  another  set  of  codes 
was  sent  in  and  the  process  repeated  until  all  possible  code  combinations 
in    (I    II    III)    were  exhausted* 

The  methods  used  in  the  unpacking  and  like-unlike  subroutines  are 
quite  straight  forward  and  will  not  be  amplified  here,    although  it  is  to  be 
noted  that  the  like-unlike  subroutine  is  not  limited  in  the  lengths  of  the 
codes  it  can  test,    up  to  the  machine   storage  size,   whereas  all  the  other 
subroutines,    since  they  use  packed  words,   do  have  quite  restrictive 
code  length  limitations.  _, 


This  program  was  first  used  to  make  an  exhaustive  search  for 
codes  of  length     10    for  test  purposes  since  all  kernels  of  this  length 
were  known.   After  this  satisfactory  checkout  the  exhaustive   search  for 
codes  of  length    26    was  undertaken. 

During  the  first  part  of  the  computer  run  for    n=26    Dr.    Golay  in- 
formed the  writer  by  personal  correspondence  that  he  had  discovered  a 
kernel  of  length    26,    using  a    "by-hand"  technique.  The  exhaustive 

search  revealed  only  a  transformation  of  the  code  pair  that  Dr.    Golay 
had  discovered.      Therefore  it  was  proved  that  only  one  complementary 
kernel  of  length    26    existed  and  that  it  was  (ignoring  allowable  trans- 
formations) 

A=0 100  11 01  11 10  101,11 100  11 10  10 
B=101 100100001 11111100111010   . 
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Supplementary  Number  of     13  Number  of     13 

Characteristic  bit  codes  with  bit  codes  with 

Number  nine  ones  seven  ones 

0  14  70 

20  14  70 

110  28  140 

1  44  112 

11  44  112 

21  44  112 

101  22  56 

111  88  224 

121  22  56 

2  14  70 

12  98  196 

22  14  70 

102  28  140 

112  70  56 

122  70  56 

13  44  112 

103  22  56 

123  35  8 

Total  715  1716 

Number  of  codes  within  each 
supplementary  block  for    n=26 

Figure     7.  2 
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CHAPTER    VIII 
A    PARTIAL    SEARCH    FOR    KERNELS    OF    LENGTH    34 

The  next  possible  kernels  for  investigation  after  length    26    were 
those  of  length     34.     As  will  be  shown  in  this  chapter  an  exhaustive 
search  for  kernels  of  length     34    was  not  feasible.      However,    an  impor- 
tant subset  of  possible  kernels  was  exhaustively  searched  and  a  scheme 
which  can  be  adapted  to  the  general  search  was  programmed.     No  ker- 
nels were  found  in  this  partial  search.      To  use  the   same  method  of 
attack  as  that  used  on    n=26    did  not  seem  too  feasible,   both  because  of 
the  length  of  time  involved  in  the   search  and  because  there  would  be  no 
guide  to   searching,   beyond  random       picking  of  possibilities  until  all 
possibilities  were  exhausted.      The  length  of  time  involved  for  an  ex- 
haustive search  of    n=34    was  estimated  fairly  reasonably  by  checking 
the  number  of  possible  codes  of  length     26    against  the  number  of  pos- 
sible codes  of  length     34,    since  the  time  for  the     26     search  was  known 

with  reasonable  accuracy.      For    n=26    there  were     71  5x715x1  71  6x1  716^* 

12 
1.  5x10         code  possibilities  at  one  point  in  the   screening  process.      At 

the   same  point  in  the   screening  process  for    n=34    codes  there  were 

12376x12376x1  9448x1 9448^  5„  8x10         possibilities,   which  is  a  ratio  of 

4 
4x10       to  one  as  compared  to    n=26„     A  time  estimate  from  the     26 

4  6 

length  code  to  the     34    length  code  was  therefore     75x4x10    =3x10       hours. 

This  length  of  time  was  of  course  not  practical  for  a  search  of  an  ex- 
haustive nature. 

The  use  of  the  cyclic  complementary  property  as  a  screen  was 

pointed  out  by  Dr.    Golay  as  highly  practical  because  some   reasonable 

19 

guess  could  be  made  as  to  which  cyclic    sets  might  contain  a  kernel. 

This  would  lead  to  a  hybrid  "by- hand"  and  computer  search.      It  would 

also  be  easier  to  document  the  areas  searched  by  the  cyclic   method  as 

compared  to  the  supplementary  method  and  would  therefore  avoid  the 
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duplication  of  effort  by  future  n=34  length  searchers.      As  will  be  dis- 
cussed at  some  length,   there  is  in  addition  an  even  more  important 
reason  leading  to  the  decision  to  substitute  the  cyclic  complementary 
property  for  the  supplementary  property  as  a  screen. 

The  generation  of  cyclic   sets  which  are  to  be  utilized  in  the  search 
for  new  kernels,    is  accomplished  in  a  somewhat  unusual  fashion.      This 
method  of  classifying,    puts  within  the  same   set  all  possible  numbers 
which  are  cyclic  permutations  of  one  another,   as  is  to  be  expected; 
however,    there  are  also  contained  within  the   same      set      all  numbers 
which  are  formed  by  the   removal  of  every  other  bit  in  a  cyclic  fashion 
until  all  bits  are  used.     As  an  example  the  five  bit  number     11010    will 
be  used  to  show  all  the  numbers  contained  within  its      set.         First, 
selecting  every  other  bit,    starting  with  the  first  bit,    and  repeating  on 

the   result  gives:        1„  11010 

2.  10011 

3.  10101 

4.  11100 
1.  11010 

and  second,   by  permuting  the  bits  of  each  of  these  in  the  normal  cyclic 

permutation  manner  yields 

1.  11010,     01 101 ,     10110,     01011,     10101 

2.       loon,    nooi,    moo,    oino,    ooiii. 

It  is  not  necessary  to  generate  the  cyclic  permutations  from    3    and    4 

since  they  are  contained  within    1     and     2. 

Counting  cyclically     the  number  of  likes  for  each  spacing  gives  the 

code  number     1331    for     1     and  the  code  number     3113    for    20      These 

code  numbers   start  with  a  spacing  of     1     and  include  up  through  a  spacing 

of    n- 1     or  four  in  this  particular  example.      The  code  number     3113 

means  that  at  a  spacing  of  one  there  were     3     likes  in     10011    counted 

cyclically,     at  a  spacing  of  two  there  was     1     like,    at  a  spacing  of  three 
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there  was     1     like,   and  at  a  spacing  of  four  there  were     3     likes.      The 

example  just  given  is   somewhat  trivial  since  all  binary  numbers  of 

length     5    with    3     ones  fall  into  the  same  sete      However,    this  is  not  the 

case  for    n=26    as  there  are     7    sets  for  codes  of  length     13     with    9 

ones  and     76     sets  for  codes  of  length     13    with     7    ones. 

An  examination  of  the  kernel  of  length    26,    shows  that  its  sequence 

quadruple  of  codes  came  not  from  four  different  cyclic  sets  as  one 

would  expect  but  came  instead  from  just  two  cyclic   sets. 

Using  the  kernel  form 

A=01 00 1101 11 10 101 11 100  11 10 10 

B  =  101100100001  1111  1100111010      one  can  decompose  this  into 

1=0010111110111 
11=0010110011101 
111=  1  101001  1101  11 
IV=  0010111100010 

Using     I    to  form  its  cyclic   set  gives: 

1.  0010111110111 

2.  0111111001101 

3.  0111011111010 

4.  0101100111111 
1.  0010111110111 

Now  take     3    and  write  it  twice,   this  gives     01110111(1101001110111) 

11010.     Picking  off  from  the  ninth  to  the  twentieth  second  bit  gives  a 

code  which  is  the  same  as     III    above  0      Therefore     I    and    III    are 

members  of  the  same  cyclic  set.     Similarly  take     II    and  form  its 

cyclic   set;    this  gives: 

1.  0010110011101 

2.  0110111001010 

3.  0111000101011 

4.  0100001110111 

5.  0001111100101 
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6.  0011011011100 

7.  0101110011010 

8.  0010100111011 

9.  0110101000111 
10.  0111011100001 
1  1.  0101001111100 
12.  0001110110110 

1.  0010110011101 

To  make  a  comparison  between    II    and    IV    the  complement  of    IV    will 
have  to  be  used  in  order  that  it  will  have     7    ones  as  does     II. 
IV=   1101000011101 

Writing    4    from  this  list  of     12    twice  gives     01000011101(1101000011101)11. 
Picking  off  from  the  twelfth  to  the  twenty  fifth  bit  gives     IV,    which  in- 
dicates that     II    and    IV    are  from  the  same  cyclic  set. 

A  check  of  kernels  of  length     10    to  see  if  their  quadruples  came 
from  cyclic   sets  in  pairs,   as  the     2  6    kernel  did,   proved  that  both  ker- 
nels of  length     10    did  in  fact  come  in  pairs.      This  was  meaningless 
however,    since  as  was  seen  in  the  example  of  length  five  with  three 
ones,    there  was  only  one  cyclic  set  possible  for  each  pair. 

In  the  hope  that  if  the  kernel  existed  in  the    n=34    case,    it  would 
also  be  formed  with    I,    III    from  one  cyclic  set  and    II,    IV  from  another 
cyclic   set,    the  program  was  modified.      This  change  along  with  giving  an 
exhaustive  search  for  code  possibilities  taken  two  at  a  time  from  cyclic 
sets  would  also  give  an  estimate  of  the  computer  run  time  for  the  ex- 
haustive search  in  general  using  the  cyclic  complementary  property. 
This  was  the  primary  reason  for  the  change  in  the  program. 

Eight  octal  digits  were  used  to  designate  the  cyclic  complementary 
number  for  each  of  the     12,376    possible  half  codes  of  length     17    with 
11     ones  and  for  each  of  the     19*448    possible  half  codes  of  length     17 
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with     10    ones.      Only  eight  digits  were  necessary  as  each  cyclic  count 
of  spacing     i    actually  counts  both  the     i     spacing  and  the     n+l-i    spacing 
simultaneously.      This  was  demonstrated  by  step    8  in  the  proof  of 
Theorem    6.2.      The  example  of  the  code  of  length  five  given  earlier  in 
this  chapter  showed  this   symmetric  property,   as  the  code  numbers  were 

1331     and    3113.      The  first  two  digits  of  these  cyclic  complementary 
designators  contain  all  the  information  available. 

The  revision  made  to  the  computer  program  was  not  too  extensive 
as  only  blocks     2,3,4,5     in  Figure     7.1     were  modified  to  handle  the 
cyclic  complementary  property  rather  than  the  supplementary  property. 
Of  course  all  blocks  were  modified  to  handle  codes  of  length    34     rather 
than  codes  of  length    26. 

One  difficulty  encountered  was  the   lack  of  computer  memory  for  the 
large  blocks  of  data,   and  it  was  necessary  to  store  the  cyclic  sets  on 
magnetic  tape  rather  than  in  the  main  memory  core  as  the     26     case  was 
handled.      The  cyclic   sets  were  then  called  in  two  at  a  time  for  checking. 
This  slowed  the  search  down  a  little,   but  by  no  more  than     1%    of  the 
total  time.      Figure     8„  1   is  a  block  diagram  of  the  process  of  obtaining 
the  cyclic  complementary  numbers  and  it  will  be  worthwhile  to  remember 
the  memory  size  restriction  when  reviewing  the  procedure  used. 

A  general  description  of  Figure     8.  1     will  show  the  differences  in 
the  method  used  in  the  search  for    n=34    kernels  as  compared  to  that 
used  for    n=26    kernels.      The  procedure  in  the    n=34    data  formation  was 
identical  with  that  of    n=26     as  far  as  the  initial  possible  number  genera- 
tion was  concerned.      (Block     1     in  Figures     7.  1     and    8e  1).      There  was 
an  immediate  departure  from  the  old  method  in  blocks     2,3,   and    4 
where  one  of  the     17    bit  numbers  with     11     ones  was  shifted  left     17    bits 
and  then  added  to  itself*      This  formed  a     34    bit  number  which  was 
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actually  the     17    bit  number  written  twice  in  time   sequence.     A  mask 
was  then  put  over  the  last     17    bits  and  this  masked  number  was  used, 
with  an  equality  search  of  the     12,3  76    generated     17    bit  numbers,   to 
identify  this  particular  number  in  the  list.     After  the  number  equal 
to  the   reference  number  was  located  in  the  generated  list  it  was  tagged 
a  one,    indicating  cyclic   set  one.      The  whole     34    bit  number  was  then 
shifted  right  one  bit  and  a  mask  used  to  again  pick  off  the  last     17    bits. 
The  search  through  the  list  and  tagging  with  the  group  number  was  re- 
peated.    After  cycling  through  all     17    possible  cyclic  codes  and  tagging 
each,    the  original    34    bit  number  was  now  repeatedly  long  right  shifted 
one  bit  and  then  right  shifted  one  bit.     A  long  right  shift  saves  the  bit 
which  is  pushed  off  the    A    register  and  stores  it  in  the    Q     register, 
while  the   right  shift  just  pushes  the  extra  bit  off  the  end  of  the     A  reg- 
ister.     Therefore  every  other  bit  of  the     34    code  was   saved  in    Q    and 
this  formed  a  new     17    bit  code  within  the   same  complementary  cyclic 
set.      This  code  was  then  formed  into  a  double  length  code  and  the  list 
of  all  possible     17    bit  numbers  with     11     ones  was   searched  and  tagged 
with  a  one  and  the  process  which  was  first  described  to  generate  cyclic 
permutations  was   repeated.      This  process  of  first  forming     17    cyclic 
permutations  of  a  code  and  then  taking  every  other  bit,   was  continued 
until  the  possible  code  which  was  originally  operated  upon  reappeared. 
Then  cyclic   set    2    was  started  and  tagged  accordingly,    and  so  on 
through    98     sets.      The     19s448  codes  of  length     17    with     7    ones  were 
classified  in  exactly  the   same  manner  and  formed     150     cyclic  comple- 
mentary possibility  sets. 

After  each  series  of  sets  were  formed  a  mask  was  used  to  search 
the  set  tag  numbers  and  pick  off  each  set  for  transfer  to  magnetic  tape. 

To  insure  that  the  pick- off  was  exhaustive,   a  zero  was  inserted  on  the 
list  in  place  of  each  removed  tagged  number. 

106 


The  sets  were  now  ready  to  be  called  back  into  the  computer  for  use 
in  the   search.     Normally  only  one  set  was  called  in  at  a  time  because 
the  pair  from  the     I,    III     set  was  checked  against  all  the     II,    IV     sets 
before  a  new  I,    III    set  was  called  into  the  computer.      The     I,    III    sets 
contained  the     11     weight  codes  while  the     II,    IV     set  contained  the     7 
weight  codes.     When  a  set  was  called  into  the  computer  each  cyclic 
permutation      subset      of    17    possible  codes  was  given  a  cyclic  comple- 
mentary number  to  categorize  the  entire   subset.     Each  group  contained 
either    4     or     8    of  these  cyclic  subsets.      These   subset  numbers  were 
formed  by  counting  the  likes  for  the     I,    III    codes  and  the  unlike s  for  the 
II,    IV     codes  for  the  eight  unique   spacmgs.      The  codes  within  the   subsets 
were  checked  against  each  other  only  if  the  likes  equaled  the  unlikes  of 
the  cyclic  complementary  subset  numbers  for  all  spacing,    or     LT    + 

LIII  =  UII  +  uiv 

Figure     8.  1     is  joined  to  Figure     7.  1     just  beyond  block    5     and  takes 

the  place  of  blocks     1,2,3,4,5    in  Figure     7.1.     The  same  feedback  paths 
shown  in     Figure     7.  1     to  insure  the  checking  of  all  codes  within  a  set 
and  to  feed  in  a  new  set  after  a  check  has  been  completed,   are  still  in 
operation  for  the   same  purpose.      Appendix    VI    contains  the  computer 
program  for  Figure     8.  1. 

The  total  run  time  for  the    n=34     search,   taking     I,    III    from  the 
same  cyclic   set  and    II,    IV    from  another  cyclic   set  was  approximately 

15    hours.      No  codes  were  found.     An  extrapolation  from  this  run  time 

5 
to  the  exhaustive  run  time  for    n  =  34    gives     15x150x98^2.2x10       hours, 

which  although  one  order  of  magnitude  less  than  the  supplementary 

method,   is  still  not  reasonable. 

A  second  partial  search  for    n=34    kernels  was  attempted  after 

observing  that  the  kernel  of  length     10 
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A=0101000011 
B=00001001 10 

with  its  center  quad  removed  yields 

A^OlOlOOll 
B^OOOOOllO 

which  is  a  complementary  pair  of  length    8.     Breaking    A1,  B1   into  its 

(I    II    III    IV)    configuration     gives 

1=0001 

11=1011 

111=0001 

IV=0100 

where     (I,    II)     and    (IIIS    IV)    are  both  complementary  pairs. 

An  attempt  was  made  to  draw  an  exact  parallel  from  this  observa- 

tion  to  codes  of  length     34.     One  possible  solution  to  equation    2.5    which 

determines  the  number  of  ones  necessary  in  each  of  the  kernels  of 

length    34    is     (21,  18).      If  this   same  equation  is  applied  to  codes  of 

length     17,    two  of  the  possible   solutions  are     (11,10)    and    (11,7).      The 

(11,  10)    pair  when  interlaced  will  total     21     ones,   and  similarly  the 

(11,  7)    pair  when  interlaced  will  total     18     ones.     Whereas  the  kernel  of 

length     10    with  the  quad  removed  had  a  sequence  quadruple  form  with 

codes  of  length    4,   the     34    code  with  the  quad  removed  would  have  a 

sequence  quadruple  form  with  codes  of  length     1 60      Appendix    II    gives 

a  list  of  all  possible  codes  of  length     16    with  weight     10    in  the    A  code 

and  weight     6    in  the     B     code.      This  list  does  not  include  the  results  of 

the  time  reversal  transformations  which  have  to  be  included  since  weight 

is  invariant  under     Tx  ,    T2  and    T„      It  was  observed  that  if  these  codes 

of  length     16    were  interlaced  to  form  codes  of  length     32    and  then  sep- 

erated  in  the  middle  and  the  quad    A=ll,    B  =  10    inserted,    all  conditions 

as  far  as  the  number  of  ones   required  in  both  interlace  codes  of     17    to 
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form  codes  of  length     34    would  be  satisfied.     An  example  might  serve 
to  clarify  the  preceding  statements. 

C=0101 1111 110001 10 

D=1001001100001010 

is  a  complementary  pair  of  length     1  60      Using  the  interlace  method  to 

form  a  pair  of  length     3  2    gives 

A=01 100011 10101111 1010000001101100 
B=001 101 101 11 110101 11 101010011 1001  . 

Splitting  this    A,B    pair  in  the  center  and  adding  the  quad  to  this  center 

gives 

A^OllOOOlllOlOllil     11     1010000001101100 
B^OOllOllOlllllOlO     10     1111010100111001    . 

Breaking  this  pair  into  standard    (I     II     III    IV)    form  gives 

1=01 011111111000110 

11=01010000111001001 

111=01011111111000110 
IV=10101111000110110  . 

I    has     11     ones.      Interlacing     I    with    II,    which  has     7    ones,    yields 
a  total  of     18    ones.     Similarly    III    has     11     ones,   and    IV    has     10    ones. 
When  these  are  interlaced  the  code  of  length     34    has     21     ones.      In  the 
actual  computer  run  the  codes  of  length     16    with     6    ones  were  converted, 
by  complementing,    into  codes  with     10     ones.     Therefore  both  the    A    and 
B    codes  contained     10    ones.     An  exhaustive  list  of  these  and  their  time 
reversals  was  generated.     Each  code  of  length     16    was  then  split  in  the 
center  and  a  zero  inserted  for  the  ninth  bit,   forming  an  over  all  code  of 
length     17    with     10     ones.     A  second  list  was  then  formed,   using  a  one 
rather  than  a  zero  for  the  inserted  ninth  bit,   forming  codes  of  length     17 
with     11     ones. 

These  two  lists  of  numbers  were  then  fed  into  the  program  of 

Figure     7.  1     at  block  two,   and  the  regular  supplementary  check  program 
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was  used.      The  codes  with     11     ones  were  used  for     (I,    III)    and  the 
codes  with     10    ones  were  used  for    (IT,    IV)0     When  these  codes  were 
combined  in  the  interlace  scheme,     I    interlaced  with     II    gave     21     ones 
and    III    interlaced  with     IV    gave     18    ones.     All  four  interlace  combina- 
tions    (I    II    III    TV),  (II    I    III  7V),(I    II    IV     III)  and    (II    I    IV     III)    were 
run  exhaustively.      The  run  time  was  approximately  three  hours  and  no 
codes  were  found.      These  partial  searches  for    n=34    are  documented 
here  because  they  represent  two  fairly  obvious  approaches  to  the  prob- 
lem,  and  it  would  be  very  wasteful  of  time  for  someone  to  duplicate  this 
effort.      The  program  as  given  in  the  appendix  can  with  a   slight  modifica- 
tion be  used  for  an  exhaustive   search  for  kernels  of  length    34    or  for  a 
partial  search  based  upon  other    "judicious  guesses". 
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CHAPTER    DC 
CONCLUSIONS 

The  study  which  has  been  presented  here  had  two  principal  objec- 
tives.     The  first  objective  was  to  make  an  exhaustive   search  for  kernels 
of  length     26    to  determine  whether  they  existed,   and  if  they  did  exist, 
to  determine  the  number  of  such  kernels.      The  second  objective  was  to 
develope  a  methodical  scheme  of  code  decomposition  from  composite 
codes  back  to  their  generating  kernels.     We  shall  summarize  in  this 
chapter  the  ways  in  which  each  of  the  major  portions  of  this  dissertation 
is   related  to  the  accomplishment  of  these  objectives.      Also  a  number  of 
suggestions  for  further  research  will  be  presented. 

The  operations  group  and  its  family  of  theorems  were  originated  to 
formalize  the  elimination  of  possible  redundancies  in  the   search  for  new 
kernels.      The  predetermination  of  the  number  of  ones  in  each  of  the 
sequence  quadruple  vectors  was  a  vital  screen  that  eliminated  many 
such  redundancies.     Also,    the   starting  of  both  the     I    and    III    vectors 
from  the   same  possible   supplementary  code  groups  eliminated  an  ad- 
ditional half  of  the  possible  codes.      The  theorems  in  the  operations 
group  chapter  also  were  useful  in  the  code  decomposition  part  of  the 
problem,    since  they  show  the  necessary  and  sufficient  conditions  for 
the     (I,    II)     and    (III,    IV)    pairs  to  be  complementary,      Theorem    3.4, 
which  states  that  if  the  sequence  quadruple  form  is     (I    II    I    TT7    then 
(I,    II)    must  be  a  complementary  pair,    is  most  helpful  since  the  process 
of  checking  for  complementarity  is  very  tedious  by  hand  methods. 

Hamming  distances  were  the  next  major  topic  considered.      They 
were  found  to  be  useful  in  a  number  of  the  proofs  in  the  Hamming  vector 
chapter.      They  were  also  a  great  aid  in  checking  hand  decomposition, 
since  it  was  easy  to  determine  whether  D(I,  II)=D(III9  IV)=n/4    for  compo- 
site codes  or  whether  they  had  the  characteristic  distance  in  the  case  of 
kernels.  Ill 


Hamming  vectors  were  the  most  powerful  tool  developed  in  this 
paper  for  code  decomposition.     The  theorems  cover  all  known  possible 
decompositions  and  also  give  strong   support  to  the  conjecture  that  ker- 
nels of  length    50    do  not  exist.      The  Hamming  vector  concept  should 
allow  a  fuller  understanding  of  complementary  sequences  in  general, 
due  to  its  characteristic  qualities  such  as  anti- symmetry,   and  also  be- 
cause there  are  only  four  possible  configurations  of  Hamming  vectors 
for  any  code  pair.     Similarly  in  sequence  quadruple  form  it  is  easy  to 
check  whether  a  decomposed  code  may  be  complementary  by  using  its 
Hamming  vector. 

Several  exhaustive  searches  for  codes  of  various  lengths  were  con- 
ducted during  this   research.      The  most  important  of  these  was  the 
search  for  kernels  of  length     26.      This   search  revealed  that  only  one 
kernel  existed.     An  exhaustive  search  for  codes  of  length     16    and    20 
revealed  that  all  codes  of  these  lengths  were  generated  from  shorter 
codes  by  standard  methods.      This  gives  considerable   strength  to  the 
conjecture  that  these  generating  methods  are  the  only  ones. 

The  possibilities  for  future  research  in  this  field  of  complementary 
sequences  are  many.     Some  of  these  follow  directly  from  this  paper 
while  others  lie  quite  far  afield  from  the  ideas  exposed  here.     An  inter- 
esting extension  of  the  present  work  would  be  an  exhaustive   search  for 
kernels  of  length    34.     It  is  felt  that  the  cyclic  complementary  screen 
which  was  used  in  the  partial  search  is  not  powerful  enough,   and  a 
better  screen  is  needed  before  this  search  can  be  conducted  on  an  ex- 
haustive basis  using  today's  computers.,      The  writer  had  considered  the 
combination  of  both  the  cyclic  and  supplementary  properties  as  screens 
in  the   same   search.     However,   a  full  study  was  not  made  of  the  gain 

which  might  be  obtained  by  this  screening  procedure.      The  programming 
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of  such  a  screen  would  not  prove  too  difficult  and,    since  the  cyclic 
property  would  logically  be  the  first  screen  of  the  two  to  be  used,    the 
areas  of  search  could  still  be  catalogued  quite  easily. 

If  it  would  be  possible  to  determine,   before  the  search,   the  charac- 
teristic Hamming  distances    D(I,    II)    and    D(II,    IV)    of  kernels,   these 
would  act  as  a  most  potent  screen  in  the  search  for  new  kernels.      It  is 
therefore  felt  that  further  research  on  characteristic  Hamming  dis- 
tances of  kernels  could  prove  quite  fruitful  for  future  applications. 

It  was  conjectured  that  kernels  of  length    50,130    and     338    do  not 

exist,    based  on  several  pieces  of  evidence,    which  although  strong  are 

not  conclusive.     Further  investigation  along  these  lines  might  offer  a 

conclusive  proof  without  the  necessity  for  an  exhaustive  search,    since 

an  exhaustive  search  for  these  length  codes  approaches  the  impossible. 

During  the  investigation  of  kernels  of  length    50    the  writer  noticed  an 

oddity  which  might  be  worthy  of  further  investigation. 

One  kernel  of  length     10    is 
A=1001010001 
B=1000000110. 

Divide  both  the    A    and    B    codes  into  segments  of  length    2,   then 
A=10  01  01  00  01 

B=10         00         00         01  10„ 

Now  let  the  symbols  10=A,01=A,      11  =  B,     and    00=Bo    Forming  a  code 

of  length     50    by  this  method  of  using  symbols  of  portions  of  the  quad 

to  represent  kernels  of  length     10    gives 

K^AAABA 


K2=ABBAA. 
This  code  pair  of  length    50     satisfies  the  parity  check,   and  also  equa- 
tion   2.  5     for  the  number  of  ones  in  both  full  codes  and  in  the  four  half 

codes.     Several  examples  were  checked  in  this  fashion  for  both  kernels 
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of  length     10,   and  all  satisfied  the  same  conditions  as  did  this  example, 
but  of  course  none  of  them  satisfied  the  necessary  condition  for  com- 
plementarity.    Along  with  helping  to  search  for  codes  of  length    50,   this 
offers  a  possibility  for  research  in  that  all  kernels   seem  to  have 
various  characteristics  in  common.     However,   the  quad  seems  more 
flexible.      Of  the  four  known  kernels,   the  quad  is  the  only  one  thus  far 
discovered  that  can  be  used  to  build  up  composite  codes  by  one  method 
of  generation  and  then  these  codes  can  be  decomposed  by  a  different 
method.      Therefore  a  study  of  kernel  characteristics  might  lead  to 
easier  ways  of  generating  new  kernels,   or  at  least  of  determining 
whether  they  exist  for  the  various  possible  lengths. 

Another  possible  field  of  research  is  cyclic  complementary  codes. 
A  communications   system  is  more  likely  to  use  a  cyclic  complementary 
carrier  modulation  instead  of  straight  complementary  modulation  due  to 
the  continuous  nature  of  the  carrier.      Therefore  an  investigation  of 
these  codes   seems  quite  in  order. 

Two  other  topics  worthy  of  research  in  the  complementary  sequence 
field  are  the  various  correlation  functions  of  complementary  sequences 
and  the  frequency  spectra  of  these  sequences.     For  secure  communica- 
tions purposes  it  is  important  to  have  a  code  that  is  noise -like.     A 
uniform  distribution  of  ones  and  zeros  gives  a  noise- like  appearance  to 
the  code.      For  example,   the  complementary  pair 

A=11010001 
B=11011110 

is  a  typical  complementary  code  pair  of  length    8.      If  the  bits  of  the    A 

code  are  considered  in  a  cyclic  fashion  one,    two,   and  three  bits  at  a 

time,    the   results  are   remarkably  uniform  as  is  seen  in  Table     9.  1. 

This  uniformity  in  output  would  make  it  extremely  difficult  for  someone 
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monitoring  the  code  to  decide  if  he  had  noise  or  an  actual  signal  on  his 
receiver.     However,    it  is   seen  that  the     B     code  is  not  nearly  so  uniform 
and  would  be  easier  to  detect  as  a  signal. 
A  code 


Bits 

Number 

Two  Bits 

Numb 

er 

Three  Bits        Nun 

1 

4 

11 

2 

111                         1 

0 

4 

10 

2 

110                       1 

01 

2 

101                        1 

00 

2 

011                        1 
100                       1 
010                        1 
001                        1 
000                       1 

B  code 

1 

6 

0 

2 

1  1 

4 

111 

2 

10 

2 

110 

2 

01 

2 

101 

2 

00 

0 

on 

2 

100 

0 

010 

0 

001 

0 

000 

0 

Noise  uniformity  of  a  complementary 
pair  of  length    80 

TABLE     9e  1 


An  investigation  of  a  large  number  of  codes  by  the  above  technique  would 
be  useful  in  determining  the  most  noise-like  code  pairs,   where  noise- 
like is  defined  as  a  uniform  output  with  bits  taken  one,   two,    three,   etc. 
at  a  time.      This  might  in  turn  lead  to  some  general  characteristic  of 
most-noise- like  complementary  sequence  pairs„  ** 
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Another  important  study  would  be  the  autocorrelation  function  of 
each  code  pair.     This  is   synonymous  with  the  output  of  a  code  detected 
by  its  matched  filter,,     For  example,   the  autocorrelation  of  the  code     A 
given  above  is 

10-30-10-18-10-10-301 
where  the  outputs  in  time  are   read  from  left  to  right,,      The  autocorrela- 
tion of  the     B     code  is 

-1030        10        18        101030-1. 
The  autocorrelation  functions  are  important  from  two  standpoints.      First, 
if  a  portion  of  the   system  should  need  to  be  shut  down  for  repairs  or 
routine  maintenance,    a  good  autocorrelation  function  would  allow  opera- 
tion on  just  one  carrier,,     Second,   the  frequency  spectrum  of  the  trans- 
mitted signal  is  the  Fourier  transform  of  this  autocorrelation  function. 
Therefore,   one  way  to  study  the  spectrum  is  first  to  autocorrelate  and 
then  to  take  the  Fourier  transform  of  the  transmitted  signals0      The  more 
widespread  the   spectrum,   the  more  difficult  is  a  jamming  procedure. 

A  study  of  the  cros  scor  relation  functions  of  complementary  sequence 
pairs  is  also  very  important  because  of  the  difficulty  of  isolating  the    A 
code     RF    carrier  from  the     B     code   receive  r0     It  was  pointed  out  in 
chapter  4  that  because  the    A    code  and  the     B    code  were   orthogonal,   at 
exact  match  cross  talk  would  be  no  problem„     However,    when  the   signals 
are  not  at  exact  match  cross  talk  might  be  a  large  source  of  noise.     An 
investigation  of  the  c  rosscorrelation  functions  of  the    A    and    B    code 
pairs  might  lead  to  a  class  of  complementary  sequences  where  the 
problem  of  crosstalk  is  minimized. 

Another  class  of  correlation  functions  which  are  important  are 

those  where  cyclic  errors  are  created  in  the   received  signal  detected  by 

the  matched  filter.      These  errors  could  be  caused  by  a  linear  shift  in 

116 


the  phase  of  the  received  signaL      There  may  be  certain  classes  of 
complementary  sequences  where  the  loss  in  coding   signal-to-noise  ratio 
is  minimized  when  noise  is  caused  by  an  error  of  this  type.      This  par- 
ticular class  of  signals  would  be  extremely  important  in  an  air  search 
radar  application. 

A  study  of  the  maximum  number  of  different  composite  codes  for 
any  length  would  be  an  important  contribution  to  the  field,   because  for 
anti-jam  reasons  it  is  very  important  to  be  able  to  change  codes,    and 
the  more  codes  with  good  characteristics  that  are  available  the  better. 

The  last  study  to  be   recommended  is  an  analysis  of  the  autocorrela- 
tion functions,   frequency  spectra,  cross  correlation  functions  ,   and 
linear  phase   shift  correlatipn  functions  for  codes  of  lengths  from     100 
to     200.      The  purpose  of  this  investigation  would  be  to  determine  if 
certain  kernels  are  capable  of  forming  composite  codes  of  better  char- 
acteristics than  other  kernels  for  radar  applications. 

In  conclusion,   the  writer  feels  that  complementary  sequence  pairs 
will  prove  very  important  in  future  applications  to  both  radar  and  com- 
munications  schemes.      However,   much  research  must  be  done  to  pick 
optimal  classes  of  these  codes  for  such  applications.      This  paper  has 
offered  many  tools  to  be  used  by  future  investigators  of  complementary 
sequences.     These  tools  have  found  important  applications  in  their 
present  state,   but  additional  investigations  along  the  lines  of  operations 
groups,    Hamming  weights  and  Hamming  vectors     should  further  in- 
crease their  usefulness. 
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GLOSSARY    OF    SYMBOLS    AND    TERMS 


(A,B) 

(I    II    III    IV) 

D(U,V) 

H(U,V) 


A_ 
Alter 


Complementary 
code 


Complement  a 
code 

Composite  code 


Cyclic   comple- 
mentary code 

Hamming  Distance 


Hamming  vector 


Hamming  weight 


Kernel 


Quad 

Sequence 

Quadruple 

Supplementary 
code 

Time  inverse 


Weight 


A  complementary  pair  of  binary  sequences, 
see  equation    2.  Z. 

A  complementary  pair  of  binary  sequences  in 
sequence  quadruple  form,    see  equation    3.  1. 

The  Hamming  distance  of  two  binary  vectors, 
see  equation    4„  1. 

The  Hamming  vector  of  two  binary  vectors, 
see  the  first  page  of  Chapter     5. 

Means  interlace  two  Hamming  vectors  in  a  pre- 
scribed manner,    see   Theorem    5.  11,    step    2. 

Is  transformed  into. 

The  time  inverse  of  a  binary  sequence  or 
vector    A. 

Means  to  complement  every  other  bit  of  both 
codes  of  a  complementary  code  pair. 

A  pair  of  binary  sequences  which  satisfy 
equation    2.  2,   also  called  a  complementary 
sequence  pair,    a  complementary  sequence 
and  sometime  just  a  pair  or  code. 

Means  change  the   sign  or  take  the  complement 
of  a  binary  sequence. 

A  complementary  code  which  is  reducible  to  a 
shorter  code        length  by  standard  methods. 

A  pair  of  binary  sequence  which  satisfy 
equation    6.  3. 

The  number  of  bits  in  which  two  binary  sequence 
differ,   notation  is     D(A,B),    see  equation    4.  1. 

A  binary  vector  or  sequence  formed  by  the 
modulo     2     sum  of  binary  sequences. 

Hamming  distance  of  a  binary  sequence  or  vector 
from  the  null  vector.     It  is  the  number  of  ones 
in  a  binary  vector. 

A  set  of  complementary  codes  which  is  irreducible. 
The  term  is   sometimes  used  for  a  member  of  the 
set. 

A  kernel  of  length    2. 

A  complementary  code  pair  expressed  as  in 
equation    3.  1 0 

A  quadruple  of  binary  sequences  which  satisfy 
equation    6.  1 . 

Means  to  make  the  first  bit  of  a  binary  sequence  the 
last  bit,  the  second  bit  next  to  last,  and  so  on;  also 
called  time  reverse. 

See  Hamming  weight. 
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EXPLANATION    OF    FORMAT     OF     THE    APPENDICIES 

Appendix    I    on  matched  filters  is  self-explanatory.,      Appendix    II 
is  the   result  of  an  exhaustive  search  for  codes  of  lengh     1  6„      There  were 
96    code  pairs  found  by  the  computer  after  the  program  screens  had 
eliminated  most  of  the  redundancies.     The  format  of  Appendix    II  gives 
the    A    code,   the     B     code,    H(A,B),    and  H[TX  (A,  B)]  which  signifies 
time  reversing  the    A    code.     It  is  worthwhile  to  note  that  each  code  pair 
and  the  pair  formed  by  the  time   reverse  of    A    were  decomposable. 
From  these     192    possible  decompositions     48    were  interlace,     48    were 
time  sequence,    and    96    were     2       special  method,   with    48     each  of  two 
different  kinds. 

Appendix    III    is  the   re_sult  of  an  exhaustive   search  for  codes  of 
length    20.      The  format  of  Appendix    III    gives  the    A    code,    the     B  code, 
if  it  is   required,    the  time   reversal  operation  used  to  put    H(A,  B)  in 
standard  form,    H(A,B)  and  from  which  of  the     2    kernels  of  length     10 
the  generating  codes  came.     All  of  the     24    code  pairs  or  their  time 
inverse  were  decomposable.      There  are     12    from  each  kernel,      8    of 
each  of  these  interlace  and  the  other    4    are  time  sequence. 

Appendix    IV     is  a  diagrm  of  all  the  subgroups  of  the  operations 
group.      It     starts       with  the  subgroup  of  order     1     and  works  up  to 
those  of  order     16.     An    x    in  the   row  indicates  that  the  operation  at 
the  top  of  the  column  is  a  member  of  the   subgroup. 

Appendix    V    is  the  CDC    1604  computer  program  in    AR    format 

and  also  in  machine  language  for  the     26     length  search.      There  are 

slight  modifications  in  the  various  sections  of  the  program  required 

for  other  length  codes.     Appendix    VI    is  also  in  both    AR    and  machine 

language,   and  is  for  codes  of  length    34.      This  program  parallels 

Figure     8.  1,   whereas  the  program  in  Appendix    V    parallels  Fig.    7.  1. 
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APPENDIX     I 
MATCHED    FILTERS 
A  Matched  Filter  is  by  definition  a  filter  which  maximizes  the  peak 
signal-to-noise   ratio.      Deterioration  of  the  signal  wave  form  is  accepted 
in  order  to  obtain  the  desired  maximum  ratio. 

One  tool  needed  for  the  derivation  of  the  matched  filter  character- 
istics is  the  Schwarz  inequality  which  is  a  special  case  of  the  Holder 

22 

inequality.  One  representation  of  the  Schwarz  inequality  is 


r  co  2foo  ?  f  °°  ? 

J        x(w)  y  (w)  dw    |   <  J      |x(w)  |    dw  J        |   y(w)    J     dw  (A  1 .  1 ) 


It  is  to  be  noted  that  the  equality  holds  when    y(w)    is  the  complex 
conjugate  of    x(w).     Rewriting  equation    A  1.  1     gives 


Poo  2 

J       x(w)  y(w)  dw    | 

~" <     1  (A  1.2) 


J        x(w)      dw  J      |y(w)   d 


2 
w 


The  derivation  for  the  matched  filter  given  here  closely  parallels 

23 

that  of  Turin. 

Let  f(t)  be  a  signal  impressed  across  a  filter  whose  frequency 
characteristic  is  H(jw)„  F(jw),  the  frequency  spectrum  of  f(t),  is 
given  by  the  Fourier  transform  of    f(t). 


L    *<*>  e   "JWtdt. 


G(jw)  is  the  frequency  spectrum  of  the  output  after  passing  through 

the  filter 

G(jw)  =  H(jw)  F(jw),  122 


or  transforming  back  to  the  time  domain  the  output  signal  voltage  is 

P  00 

g(t)  =    J^HOw)  F(jw)  e  JWtdf.  (Al.  3) 

Since  at  some  time,     t  =  A,     g(t)  must  be  a  maximum, 

max  g(t)  =  g(A)  =    J^HtJw)  F(jw)  e  JwAdf.  (A  1.4) 

To  complete  the  derivation  it  will  be  necessary  to  obtain  the  total 
noise  power,   and  the  total  power  in  the   signal.     If  the  noise  is  assumed 
to  be  white  noise  of    N    watts  per  cycle,   the  output  noise  power  density  is 

N     =  N    I   H(jw)    |Z, 
and  the  output  noise  power  is  therefore 

00  p       00 

J^     N  df  =  J^   N|H(jw)|    df„  (A  1.5) 

The  total  energy  contained  in  the   signal  is 


E  =  J 


which  is  also 


00     2 
«    f  (Odt, 


E  =  J        |F(jw)  I    df,  (A  1.6) 

_00  I  «•>         '    I 

by  Plancherel' s  Theorem.      It  is  to  be  noted  that    E     is  a  constant  not 

dependent  upon  the  filter  used. 

The  output  power  signal- to- noise  ratio  for  the  maximum  output  is 

therefore 
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(g(A)|2       =       1   L     F(jw)  H(jw)  eJwAdf 
No  nJ        |H(jw)|2df' 


Dividing  this  equation  by  a  constant  will  not  effect  the  time  of  maximum 
signal  output;     therefore  dividing  by    E,    equation    A  1.  6  ,    gives 


N 


g(A)    |2    .       |  J        F(jw)  H(jw)  e  jwA  df   |    2 


E  N 
o 


H(jw)|2df     f        |   F(jw)|2  df  (A  1.7) 


If  the  right  hand  side  of  equation     (A  1.  7)     is  compared  to  the 
Schwarz  inequality    (A  1„2)    with    x(w)  =  F(jw)  e  J  y(w)  =  H(jw), 

the  two  expressions  are  the  same.      But  for  the  expression  to  be  max- 
imum   y(w)     must  equal  the  complex  conjugate  of    x(w),    therefore 
H(jw)  =  F(-jw)  e  -jwA. 

A  matched  filter  is  therefore  a  filter  whose  frequency  character- 
istic is  the  complex  conjugate  of  the   signal  spectrum  to  which  is  matched. 
The  output  signal  for  a  given  signal  applied  to  a  matched  filter  is 
the  refore 

oo 

/^\  I  I  TT/   ■         \    |2  jw(t-  A) 

g(t)  =  JM     jF(jw)  I    eJ  'df  ^ 

at    t  =   A     the  output  is 

00 

g(A)  =     ]m      |F(jw|2  df  . 

This  is  the   same  as  equation    A  1„  6    which  is  all  the  energy  con- 
tained in  the   signal.      Therefore  at  exact  match  all  the  energy  in  the 
pulse  is  in  the  output  signal. 
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APPENDIX    II 

Exhaustive  Search  for  Codes 
of  Length     1  6„ 
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APPENDIX    III 

Exhaustive  Search  for  Codes 
of  Length     20. 
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APPENDIX   IV 

Subgroups  of  the  Operations  Group 
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APPENDIX    V 
REMARKS    ON    COMPUTER    PROGRAM 

In  Appendix    V    the  addresses  from  60000     to     60035    are  block     1 
of  Figure     7.  1.     All  numbers  used  in  the  program  are  octal.      The  program 
as  written  is  for  codes  of  length    26,    so  most  of  the  program  deals  with 
half  length  codes  of  length     13.      Words  to  be  changed  for  different  length 
codes  are     60000  (mask  of  ones  for  half  length  code  as  determined  by- 
weight),     60001     (length  of  half  code  minus  one   [upper  half],   and  length 
of  memory  word  minus  length  of  code   [lower  half]),    60004  (length  of 
half  code),    60034  (last  sequence  in  code  list).      For  exactly  the   same 
reasons  respectively  change     60015,    60016,    60021     and    60032. 

Addresses     17000     to     17045     separate  the  sequences  into  blocks  by 
supplementary  number.      Ir\dex  register    2     must  be  entered  with     100 
before   starting.     Addresses  to  be  changed  here  are     17000  (mask  of  first 
3    and  last     3    bits),    17005  (contains  the  total  number  of  sequences  for 
both  half  length  codes),    17007  (contains  the  number  of  sequences  with 
7    ones),    17015  (contains  a  mask  for  half  length  codes),      17026    (contains 
the  number  of  sequences  with     9    ones),    17032,  17033     (have  for  their 
instruction  address  the  base  word     50000     plus  the  number  of  half  codes 
with     7    ones),    17034    (contains  a  mask  for  code  length). 

Addresses     10000     to     10106     are  the  repetitive  part  of  the  program. 
The  only  memory  address  needing  modification  here  for  different  length 
codes  is     10106    which  contains     SVN    a  mask  equal  to  the  one  half  code 
length. 

Addresses     7200    to     7252    contain  the  unpacking  subroutine  and  the 
like-unlike   subroutine.      7201,7207,7215,     and     7223    all  contain  the 
half  code  length.     Similarly    7204,7207,7220,7223    and    7230    all  contain 
the  full  code  length  or  the  length  modified  by  one  or  two  to  fit  the  program 
needs.  135 


The  15000     to     15100     series  is  all  possible  bit  arrangements  for 
the  first     3    and  the  last     3     bits  of  codes  of  length     13,     while     16000    to 
16100     is  a  list  of  supplementary  block  numbers  in  one  to  one  correspon- 
dence with  this   list. 

The     10300    and     10400     series  respectively  contain  the  number  of 
codes  in  each  supplementary  block  for  the  half  length  codes  of  weight 
9     and  weight     7. 

16000     to     16020     is  the  final  list  of  the     18     supplementary  block 
numbers  and  is  in  one  to  one  correspondence  with  the     10300    and     10400 
series. 
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